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NAVFAC  Instruction  11012.138,  24  October  1975,  has  specified  that, 
for  buildings  designated  as  "critical  buildings,"  a site  seismicity 
study  be  performed,  in  lieu  of  the  lateral  load  coefficients  provided  in 
NAVFAC  P355,  Seismic  Design  of  Buildings,  April  1973.  Acceleration 
response  spectra  developed  from  the  study  shall  be  used  to  determine  the 
magnitude  of  the  lateral  loads.  This  report  is  intended  to  provide  a 
methodology  and  automated  procedure  for  computing  site  microzonation 
studies  and  recurrence  data  needed  for  design. 

This  report  cites  the  works  of  many  authors.  Much  of  the  material 
is  included  here  to  give  the  reader  adequate  background  to  understand 
the  developments  and  assumptions  made.  In  particular,  the  work  of 
Krinitzsky  (1974)  was  cited  verbatim  in  many  places. 

The  report  is  organized  into  nine  chapters.  Chapter  1 gives  the 
reader  an  introduction  into  earthquake  causes,  types  of  faults  and 
associated  behavior.  Chapter  2 gives  dimensional  relationships  between 
faults  and  magnitude.  Such  information  is  used  in  later  sections. 
Chapter  discusses  fault  movement  and  geologic  slip  rates  as  a proce- 
dure for  estimating  earthquake  recurrence.  Chapter  4 discusses  the 
historic  earthquake  data  base  and  its  limitations.  Chapter  5 presents 
stochastic  procedures  for  treating  earthquake  data.  Chapter  6 presents 
earthquake  site-induced  ground  motion  as  a function  of  epicentral  dis- 
tance and  magnitude.  Chapter  7 draws  upon  the  previous  sections  as 
background  and  presents  an  automated  procedure  for  determining  a prob- 
abilistic estimation  of  site  acceleration  based  on  the  historic  data 
base  and  adjusted  to  consider  geologic  data.  The  confidence  limits  of 
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the  estimate  are  given*  Such  a process  is  the  first  step  in  specifica- 
tion of  a design  level  earthquake  loading.  Chapter  8 provides  the 
second  step  - the  preparation  of  response  spectra.  Chapter  9 discusses 
risk  analysis  and  an  approach  for  econosuc  selection  of  design  criteria. 
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CHAPTER  1 - EARTHQUAKES  AND  FAULTING 

Plate  Tectonics 

The  United  States  is  located  on  the  North  America  plate,  the  west- 
ern portion  of  which  mets  the  Pacific  plate.  The  interaction  of  these 
two  plates  is  responsible  for  the  high  seismic  activity  which  has  in  the 
past  and  continues  now  to  take  place  in  the  Western  United  States. 
Plate  tectonic  theory  has  explained  much  of  the  geologic  activity.  Also 
explained  in  terms  of  plate  tectonic  theory  is  the  seismic  activity 
experienced  in  the  Central  and  Eastern  parts  of  the  United  States.  This 
midplate  activity  can  be  very  destructive. 

Figure  1-1  shows  a cross  section  of  the  earth.  The  lithosphere, 
composed  largely  of  basalt,  extends  to  an  average  depth  of  about  100  km. 
Below  the  lithosphere  is  the  asthenosphere  which  extends  to  a depth  of 
400  km.  Because  its  upper  portion  is  partially  molten,  seismic  velocity 
in  that  region  is  decreased.  The  lithospheric  plates  are  able  to 
"float"  on  this  plastic  layer.  Not  all  of  the  asthenosphere  is  molten; 
however,  there  is  a ligid  portion. 

Most  seismic  activity  is  located  at  plate  boundaries  and,  there- 
fore, boundaries  are  of  considerable  interest.  Walper  (1976)  is  quoted 
in  the  following  paragraph: 

Three  kinds  of  plate  boundaries  can  be  distinguished 
by  the  three  principal  modes  of  interaction  between  adjacent 
plates,  and  distinctive  types  of  plate  margins  are  identified 
with  each:  (a)  divergent  or  spreading  boundaries  with  con- 
structive or  accreting  plate  margins;  (b)  convergent  or  colli- 
sion boundaries  with  destructive  or  consuming  margins;  and 
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crust  [Figure  1-2 1 . 
Geologic  Faults  and  Earthquakes 


Since  the  San  Francisco  earthquake  of  1906  and  the  subsequent  work 
on  the  elastic  rebound  theory  of  earthquakes,  general  agreement  has  been 
reached  on  the  close  relationship  between  earthquakes  and  geologic 
faults.  Most  tectonic  earthquakes  of  the  type  that  cause  major  struc- 
tural damage  are  associated  with  fracture  on  a fault.  Plate  motion 
causes  stress  in  the  earth's  rock  crust.  Earthquakes  occur  when  the 
strength  of  the  fault  can  no  longer  withstand  the  stress  that  has  built 
up.  Fault  plane  solutions  and  earthquake  mechanism  studies  have  contri- 
buted to  a consistent  picture  of  the  earthquake  generation  process  which 
satisfactorily  explains  most  of  the  observed  facts. 

Krinitzsky  (1974)  states: 

A fault  is  a shear  rupture  in  the  earth  along  which 
opposite  faces  have  been  displaced  relative  to  each  other. 
Relative  motion  is  parallel  to  the  fault  itself.  ***  t A 
rupture  without  relative  displacement  is  termed  a joint . 

The  basic  kinds  of  faults  discussed  above  are  illus- 
trated schematically  in  ***  (Figure  1-3  showing  appropriate 
terminology).  They  are  defined  in  the  following  paragraphs. 

a.  Strike  slip.  Strike  is  the  line  of  direction  in  a 
horizontal  plane  produced  by  intersection  by  a fault  plane. 
Movement  along  this  line  ***  (Figure  l-3a]  is  strike  slip. 

Such  faults  are  usually  high  angle,  that  is  near  to  vertical. 

***  The  direction  of  movement  is  designated  as  right  lateral 


^Asterisks  (***)  indicate  deleted  material  in  quoted  passages. 


4 


or  left  lateral  depending  on  the  direction  (right  or  left)  the 
opposite  side  hss  moved  when  one  faces  the  fault.  ***  The 
soveaent  shown  in  ***  [Figure  1 -3a ) is  left  lateral. 

b.  No  real . A normal  fault  ***  [Figure  l-3b)  is  one  in 
which  the  hanging  wall  (the  overlying  side  of  the  fault)  has 
moved  downward  relative  to  the  footwall  (.the  underlying  side 
of  the  fault).  For  hanging  wall  and  footwall  ***  [see  Figure 
1-4).  This  movement  results  from  horizontal  tension.  ***  It 
is  also  called  a gravity  fault.  Note  that  AB  in  ***  (Figure 

1 - 3b ) is  slip,  AC  is  throw,  and  BC  is  heave. 

c.  Thrust  or  reverse.  This  is  a fault  ***  [Figure  l-3c) 
in  which  the  hanging  wall  has  moved  upward  relative  to  the 
footwall . It  is  caused  by  horizontal  compression.  ***  Its 
movement  AB  is  slip,  AC  is  throw,  and  BC  is  heave. 

Host  faults  are  combinations  of  strike  slip  and  one 
or  the  other  of  the  above.  Also,  the  fault  plane  itself  is 
likely  to  be  curved.  Blocks  may  be  rotated  relative  to  each 
other. 

Components  of  movement  ***  [Figure  1 - 3d ) are  staged 
as  AB,  net  slip;  AE  or  CB,  strike  slip;  AC  or  EB,  dip  slip; 
AD,  throw;  and  DC,  heave. 

The  fault  trace  ***  [Figure  1-4)  is  the  line  of  the 
fault  at  the  surface.  The  strike  of  the  fault  is  measured  in 
a horizontal  plane  as  degrees  of  deviation  of  its  direction 
from  north.  The  dip  is  measured  in  degrees  from  the  horizon- 
tal plane  to  the  fault  plane  and  is  given  a directional  compo- 
nent that  is  perpendicular  to  the  fault  strike. 

Features  commonly  associated  with  normal  faults  are 
shown  schematically  in  ***  (Figure  1-5).  Step  faults  are 
blocks  that  are  sheared  like  slices.  The  lowermost  block  is  a 
graben.  An  isolated  high  block  is  a horst. 


The  movements  illustrated  in  ***  [Figure  1-3)  are 
highly  siaplified.  In  nature  they  may  be  very  complex  with 
secondary  and  tertiary  slippages  related  to  primary  fault 
movements . 

Some  misunderstandings  have  occurred  - and  perhaps  some  significant 
differences  of  opinion  - about  the  direct  relationship  between  geologic 
faults  and  the  earthquake  hazard.  There  is  much  evidence  dating  from 
1906  to  suggest  that  destructive  ground  shaking  is  not  necessarily  at  a 
maximum  in  the  immediate  vicinity  of  the  causative  fault  (Hudson,  1972). 
More  often  than  not,  the  maximum  destructive  ground  shaking  is  observed 
to  be  some  miles  from  the  fault,  as  can  be  explained  by  a number  of  the 
features  of  the  generation  and  propagation  of  seismic  waves.  Classical 
photographs  of  the  1906  major  movements  along  the  San  Andreas  fault,  for 
example,  show  horizontal  surface  displacements  of  as  much  as  IS  feet 
passing  several  feet  from  a small  wood-frame  house  that  received  no 
significant  damage.  Similarly,  during  the  San  Fernando  earthquake  of 
9 February  1971,  a 5-foot  vertical  fault  scarp  passed  directly  through  a 
wooden  barn  just  a few  hundred  feet  from  a single-story  residence.  The 
barn  was  severely  damaged,  but  no  significant  structural  damage  to  the 
house  was  noted  (Hudson,  1972).  The  San  Fernando  earthquake  also  fur- 
nished numerous  examples  of  surface  faulting  passing  through  heavily 
populated  areas.  Although  severe  structural  deformation,  with  a result- 
ing economic  loss,  occurred  in  numerous  cases,  catastrophic  collapses 
leading  to  loss  of  life  and  serious  injury  were  not  directly  associated 
with  these  surface  breaks.  Hazardous  collapses  were  in  all  cases  the 
consequence  of  severe  ground  shaking,  which  is  pervasive  over  a large 
area  and  is  not  limited  to  the  vicinity  of  faults. 

The  focus  or  hypocenter  is  the  point  within  the  earth's  crust  where 
the  initial  rupture  occurs  and  from  which  the  first  waves  are  released. 
The  projection  of  this  point  to  the  ground  surface  is  the  epicenter. 
The  epicenter  and  hypocenter  do  not  necessarily  indicate  the  center  of 
total  energy  release  of  the  earthquake  but  rather  the  point  where  the 


seismic  energy  waves  were  first  created.  For  small  earthquakes  the 
center  of  total  energy  release  and  the  epicenter  are  not  far  apart 
because  the  fault  break  length  is  short;  however,  this  is  not  the  case 
for  large  earthquakes.  The  majority  of  earthquakes  in  the  United  States 
have  had  relatively  shallow  focal  depths  (0  to  40  km).  In  California, 
earthquakes  have  occurred  in  regions  where  surface  fault  patterns  were 
clearly  visible.  In  the  Puget  Sound  area,  earthquakes  are  focused  at 
deeper  locations  within  the  earth's  crust  so  that  a surface  rupture  is 
not  observable.  In  the  eastern  United  States,  the  relationship  to 
surface  rupture  has  not  been  identified  (Bolt,  1970;  Newinark  and 
Rosenblueth,  1971). 

A fault  undergoing  tectonic  creep,  or  one  with  abrupt  displacement, 
causes  changes  in  the  terrain  it  crosses.  Very  distinctive  patterns  are 
produced  where  active  faults  cross  streams,  such  as  landslides.  The 
ongoing  geologic  process  causes  scarps,  trenches,  sag  ponds  and  stream 
offsets  Figure  1-6  shows  a landform  with  an  active  fault  (Wesson  et 
al.,  1975) . 

Estimates  of  the  maximum  sire  and  frequency  of  earthquakes  on  a 
fault  are  based  on  the  geologically  determined  slip  rate  and  the  histor- 
ic record  of  ground  deformation  (where  available),  the  seismic  history 
of  the  fault  and  surrounding  tectonic  region,  a geological  evaluation  of 
the  tectonic  setting,  and  empirically  derived  relationships  between 
earthquake  magnitude  and  fault  length. 


Surface  Effects  of  Fault  Movements 


When  faults  are  considered,  the  assumption  is  commonly  made  that 
the  creation  of  entirely  new  faults  by  an  earthquake  is  unlikely 
(Krinitzsky,  1974).  Significant  surface  faults  and  their  activity  may 
be  found  by  proper  geologic  investigations.  Cluff,  Slemmons,  and 
Waggoner  (1970)  have  illustrated  schematically  the  character  of  typical 
surface  effects  by  faults.  These  are  shown  in  Figures  1-7  to  1-9. 


Strike-Slip  Faults.  Strike-slip  faults  are  usually  steep,  nearly 
vertical.  Figure  1-7  shows  the  damage  that  is  associated  with  strike- 
slip  movement.  Some  fault  displacement  pre-exists,  and  movement  occurs 
along  the  pre-existing  plane.  There  is  very  little  vertical  component 
of  movement;  hence,  there  are  negligible  topographic  effects  other  than 
the  slip  itself.  Krinitzsky  (1974)  describes  types  of  faults  below: 

The  San  Andreas  is  a nearly  vertical  strike-slip 
fault  that  in  recent  times  has  had  a cumulative  displacement 
of  many  kilometers.  It  has  manifested  this  activity  by  pro- 
nounced effects  on  the  (regional]  drainage.  Where  it  passes 
through  mountainous  country,  for  example  in  the  Santa  Cruz 
Mountains,  the  drainage  has  cut  a valley  along  the  active 
trace  of  the  fault.  In  the  Carrizo  Plain  of  southern 
California,  ***  the  fault  has  had  a dramatic  effect  in  turning 
drainage  that  crosses  the  fault.  The  displaced  drainage  now 
is  joined  in  a dogleg  along  the  fault  trace.  This  is  a mea- 
sure of  the  displacement  that  has  occurred  since  the  drainage 
was  developed.  Though  drainage  history  has  never  been  dated 
precisely,  interpretations  are  possible  of  the  manner  of 
movement  in  recent  times. 

Strike  slip  faults  [usually]  tend  to  have  greater 
continuity  along  their  fault  traces  than  do  other  faults. 

Also,  strike-slip  faults  will  occur  within  narrower  fault 
zones.  Along  the  San  Andreas  the  full  zone  in  some  places  is 
as  much  as  1.6  to  3.3  km  wide.  However,  the  most  recent  and 
most  pronounced  fault  movement  is  likely  to  be  on  single, 
identifiable  planes.  A recurrence  of  movement  has  a high 
probability  of  being  along  the  same  plane.  Thus,  the  trace  of 
movement  along  a strike-slip  fault  is  predictable  to  a far 
larger  extent  than  for  other  types  of  faults. 
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Vartoua  portion*  of  Krtnitxaky  (IU74)  are  quoted  below: 

land* I idea.  Wham  fault  movement*  are  poatnlated 
nr  at  an  eng  I itrr  r tug  aitr,  thr  Uacartl  from  landalide*  ahoold  hr 
evaluated  aa  a apeciat  category.  Thay  taay  ha  produced  in 
aaaoriat  ion  with  fault*  or,  atora  commonly,  thay  may  hr  trig- 
garad  elaewhere  by  earthquake*.  l.onaa,  aaturatad  dahria,  or 
jotntad  rook  with  weakened  plana* , on  ataap  alopaa  are  aua- 
ceptihle  and  havr  contributed  to  apart  ecu 1 ar  laudalidr*  in  I hr 
paat , 

Normal  Fault*.  Noratal  fault*  ***  (Figure  l-8|  Itavr 
more  pronounced  topographic  effect*.  Thr  downdroppe*.  block 
hraaka  along  thr  dragged  lip  and  fort**  aacondary  diaplacnarnt  a 
in  addition  to  thoaa  on  thr  major  plana.  Damage  by  fault 
dtaplacrmrnt  i*  concent  rated  in  thr  downdropped  block.  The 
opt  brown  block  may  remain  relatively  intact.  Tltia  aaHttatra 
that  movrmant  ia  confined  to  a aiuglr  major  fault  plana. 
Normal  faulting  ia  auacaptiblr  alao  to  mu 1 1 1 pi  e-plane  aleplike 
d tap  I at  amenta. 

Thruat  Fault*.  Thruat  fault*  ***  | Figure  1-U|  tend 
tc.  break  up  in  the  up  thrown  block.  Thr  downthrown  block 

remain*  intact,  except  that  near  the  face  of  the  major  fault 
plane  there  ia  landaliding.  Field  meaaurement  of  the  amount 
of  thru*ting  often  i*  difficult  bacauae  of  tltia  deal  ruction  of 
the  fault  plana.  The  breakage  in  the  upthrown  block  tend*  to 

be  highly  arcuate  and  irregular.  Throat  fault*  are  likely  to 

occur  along  multiple  plane*, 

Thruat  fault*  may  occur  along  a mountain  front  by 
activation  of  a acarp* forming  proceaa,  auch  aa  that  which 
occurred  during  the  San  Fernando  earthquake  of  1U7I,  Such 
movement  ia  to  aome  degree  predictable.  However,  thruat 

fault*  aiay  alao  occur  in  vary  complex  and  often  unpredictable 
aaaoctat ton* . A thruat  fault  occurred  at  Hanning  Hay  during 


the  Alaska  earthquake  of  1964.  It  raised  the  bay  floor  and 
cut  across  an  adjacent  ridge.  The  ridge  contains  a notch  that 
had  been  developed  during  earlier  movements. 

Geological,  Sei sinological , and  Soils  Investigations  and  Fault  Activity 

The  objective  of  the  geological  investigation  is  to  establish  the 
lithology,  stratigraphy,  structure,  and  history  of  the  general  region. 
Tectonic  structures  underlying  the  region  must  be  identified,  and  any 
evidence  of  fault  activity  determined.  The  seismologica 1 investigation 
compiles  a listing  of  all  earthquakes  of  record  which  may  have  affected 
the  general  area  of  the  site.  The  magnitude  of  the  earthquakes,  epi- 
center locations,  dynamic  characteristics,  and  durations  of  the  re- 
sulting ground  motion  are  determined.  Epicenters  within  about  200  miles 
of  the  proposed  site  are  of  particular  significance.  Geological  fault 
structures  within  this  approximate  radius  should  be  studied. 

The  type  of  information  required  includes  identification  and  delin- 
eation of  those  faults  that  are  active;  that  is,  capable  of  generating 
earthquakes.  Information  about  recency  of  faulting  is  useful  in  desig- 
nating particular  faults  as  active  for  specific  land  uses.  Designation 
of  a fault  as  active  can  be  based  on  the  assumption  that  the  more  recent 
the  faulting  (geologic  time)  the  more  likely  that  the  fault  will  undergo 
intermittent  displacement  in  the  geologically  near  future.  Movement  is 
presumed  less  likely  to  occur  in  the  immediate  future  (50  years)  along 
faults  that  have  progressively  longer  geologic  periods  of  demonstrated 
quiescence.  Because  of  the  apparent  great  range  in  frequency  of  move- 
ment, there  is  no  agreement  at  present  on  the  length  of  geologic  time 
pertinent  to  evaluation  of  the  near-future  behavior  of  faults.  Selec- 
tion of  the  timespan  used  to  designate  faults  as  active  from  age  of 
latest  movement  has  been  influenced  by  the  potential  consequence  of 
seismic  shaking  of  surface  faulting  on  the  specific  structure.  The 
greater  the  risk  to  be  incurred,  the  longer  the  timespan  that  must  be 
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considered.  Displacement  during  Holocene  time*  is  a generally  screpted 
criterion  of  activity  for  many  land  uses.  This  timespau  is  probably 
inadequate,  however,  to  assure  recognition  of  all  active  faults.  His- 
toric offsets  have  occurred  along  faults,  such  as  the  White  Wolf  fault 
in  southern  California,  that  had  no  previously  recognised  evidence  of 
Holocene  faulting.  Wentworth  and  Yerkes  (1971)  state  that  evidence  of 
displacement  during  late  Quaternary  time  (the  past  several  hundred 
thousand  years)  should  be  considered  evidence  that  a fault  is  probably 
active  for  low  seismicity  regions. 

To  quote  Krinitxsky  ( 1 9 7 A ) further: 


Fault  movements  have  occurred  throughout  geologic 
time.  The  oldest  geologic  strata,  it  may  be  supposed,  are  the 
most  faulted  but  this  is  not  a general  condition.  However, 
faults  which  moved  in  the  past  will  not  necessarily  move  again 
in  the  present. 

***  faults  l,for  an  engineering  analysis,!  are 
either  active  or  inactive.  Active  means  that  a fault  may  move 
at  some  time  in  the  near  future.  In  engineering  |terms|,  it 
means  it  may  occur  during  the  life  of  a structure.  Inactive 
means  that  it  will  not. 

***  The  Atomic  Energy  Commission  | 1973b)  has  devel- 
oped a set  of  guidelines.  ***: 

a.  Datable  movement  during  the  past  35,000  yr.  (The 
limit  of  reliable  radiocarbon  dating.) 

b.  Datable  movement  more  than  once  in  the  past  500,000 
yr.  (Narine  terraces  on  the  Weat  Coast.) 

c.  Structural  interrelation  whereby  a fault  can  he  shown 
to  move  if  movement  occurs  on  a different  fault  with 
proven  activity. 


The  most  recent  geologic  time  unit  - beginning  of  the  last  glacial 
retreat  12,000  to  20,000  years  ago  to  the  present. 


d.  Instrumentally  determined  macroseismic  activity 
relatable  to  a fault. 

e.  Projection  of  a proven  active  fault  through  or  into 
areas  where  all  evidence  of  the  fault  or  its  activity 
is  obscured,  as  by  thick  alluvium. 

The  International  Atomic  Energy  Agency  1 1972 ] adds 
two  additional  criteria: 

a.  Evidence  of  creep  movement  along  a fault  (creep  is 
slow  displacement  not  necessarily  accompanied  by 
macroearthquakes) . 

b.  Topographic  evidence  of  surface  rupture,  surface 
warping,  or  offset  of  geomorphic  features. 

They  would  further  classify  active  faults  on  a geomorphic 
basis  as  follows: 

Class  A - High  rate  of  movement,  greater  than  1 m per 
1000  yr 

Class  B - Topography  shows  clear  evidence  of  dislocation 

Class  C - Topography  shows  indistinct  evidence  of  dislo- 
cation 

Class  D - No  evidence  of  amount  or  rate  of  dislocation  on 
which  quantitative  assessment  can  be  based,  but 
fault  is  considered  capable  of  causing  surface 
faulting. 

In  general  engineering  practice  a fault  is  considered  to 
be  active  if  there  is  displacement  within  Holocene  ***  depos- 
its regardless  of  datable  evidence.  (Examples:  fault  dis- 
placements within  surface  gravels,  alluvium,  or  glacial  out- 
wash.)  ***  All  of  the  criteria  above  are  indicative  of  a 
present-day  hazard  from  renewed  movement  along  a fault,  and 
they  are  sufficient  to  identify  a fault  as  active. 
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Geophysical  and  Geological  Investigation 


1 


A general  description  of  the  types  of  tasks  performed  in  geophysi- 
cal and  geological  investigations  will  be  given  here  as  given  in 
Krinitzsky  (1974). 


Background  Information. 

Background  information,  ***  which  ***  [may 
assembled  and  studied  ***  includes: 


a.  Imagery.  Air  photos  and  other  [photos]. 

b.  General  geology.  Available  maps  and  reports  on  the 
rock  and  soil  types,  ages,  tectonic  history. 

c . Seismicity.  Historic  record. 

d.  Mapped  faults. 


Air  Inspection.  An  efficient  way  to  begin  is  to 
overfly  the  area  in  a helicopter  or  small  plane  and  inspect 
the  terrain  carefully  for  whatever  evidences  can  be  found. 
Early  morning  and  late  afternoon  hours  are  best,  since  long 
shadows  may  reveal  faults  that  might  otherwise  be  missed.  The 
faults  should  be  scanned  for  appearances  of  fresh  displace- 
ments and  a notation  made  for  checking  the  displacements  on 
the  ground.  The  prime  concern  in  checking  faults  for  earth- 
quake risk  is  to  determine  the  presence  or  absence  of  activity 
along  the  faults. 

Ground  Inspection.  Criteria  for  determining  recent 
fault  activity  are: 

a.  Scarps . Freshness  of  fault  breaks;  displacements  of 
surficial  gravels  or  soil  layers. 

b.  Turned  drainage. 

c.  Sag  ponds. 
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d.  Other  hydrologic  anomalies.  Springs;  lines  of 
greenery  in  arid  areas;  groundwater  barriers. 

e.  Land  displacement.  Anomalous  displaced  land  blocks; 
drowned  trees. 

f.  Tilted  trees.  Bends  in  trees  along  fault  zones. 
Bends  may  be  related  to  recent  movements.  Tree  rings 
may  help  to  date  fault  movements. 

g.  Sandblows . To  determine  inclusive  areas  of  large 
earthquake  disturbance. 

h.  Harrowed  ground.  Intricate  breakage  of  ground  by 
earth  shaking;  usually  on  local  topographic  highs. 

Associations . The  information  above  must  be  placed 
in  context  with: 

a.  Historic  seismic  events.  A caution  is  that  reported 
epicenters  are  sometimes  unreliable.  They  should  be 
rechecked  wherever  they  contribute  to  important 
judgments  on  location  of  fault  movements. 

b.  Structural  trends.  An  area  of  doubtful  activity  may 
be  determined  to  be  a continuation  of  a structural 
trend  in  which  elsewhere  there  is  proven  activity. 
Activity  may  be  inferred  by  association  and  through 
projection. 

Dating.  Dating  calls  for  imaginative  and  resource- 
ful approaches.  Some  of  the  more  common  are: 

a.  Stratigraphic  displacements.  Borings  may  show  net 
displacements  in  young  deposits  of  known  age.  Rates 
of  movement  may  be  postulated  for  such  faults. 

b.  Effects  of  datable  geomorphic  features.  Breaks  in 
marine  terraces  of  known  ages.  Displacements  in 
datable  alluvial  sediments  or  glacial  deposits. 
Similarly,  these  may  give  average  rates  of  fault 
movement . 

c.  Paleontology  and  palynology.  Age  dating  of  fossils 
in  displaced  layers. 
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d.  Radiocarbon  dating.  Age  dating  of  organic  matter 
(peat,  wood,  bones,  etc.)  in  sag  ponds,  scarp 
detritus,  and  displaced  layers.  Often  such  dating 
can  be  related  to  specific  earthquake  events. 

e.  Paleomagnetism.  Magnetic  dating  with  relation  to 
polar  wandering  during  recent  times. 

f.  Uranium  disequilibrium.  Disintegration  of  radio- 
active uranium  with  time. 

g.  Potassium-argon.  Particularly  useful  for  displaced 
volcanic  rocks. 

h.  Other.  Associations  with  young  cinder  cones, 
presence  of  fresh  scarps  without  datable  associa- 
tions, etc. 

Mapping.  Mapping  should  include  as  much  as  can  be 
determined  of  the  following  information: 

a.  Geometry  of  fault  movement. 

b.  Age  of  most  recent  movement. 

£ . Frequency  of  movement  or  rate  of  movement. 

d.  Reliability  of  interpretation. 

e.  Projected  activity  beneath  masking  deposits. 

f.  Association  with  landslides.  In  particular,  one  must 
be  certain  that  landslides  are  not  mistaken  for  fault 
movements . 


m 
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Strain  Measurement.  This  is  principally  a research  technique 
consisting  of  installation  of  strain  measuring  gages. 

Geodetic  Surveys.  Precision  surveying  of  a series  of  monuments  on 
opposite  sides  of  a fault  is  used  to  measure  fault  isovement. 

Geophysical  investigations  include  the  following: 

Bouguer  Gravity  Survey.  A regional  survey  of  the  earth's  gravity 
in  which  corrections  have  been  made  for  the  effects  of  topography. 
Anomalies  are  indicative  of  changes  in  the  mass  properties  of  rocks 
where  something  unusual  has  happened  in  the  earth's  crust.  This 
approach  does  not  identify  specific  faults,  but  may  help  delineate  areas 
of  faulting  in  the  subsurface. 

Magnetic  Survey.  A magnetic  survey  explores  local  changes  in  the 
earth's  magnetic  field.  The  results  are  affected  by  changes  of  rock 
types  (iron  content)  revealing  abrupt  changes  such  as  those  which  take 
place  where  there  is  movement. 

Seismographic  and  Electrical  Resistivity  Surveys.  These  are  the 
various  techniques  for  passing  energy  through  rocks  and  monitoring 
changes  in  rock  properties.  The  property  changes  occur  from  either 
stratification  or  inhomogeneities  in  the  rock  mass  or  from  displacements 
in  the  rock  mass. 
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Schematic  cross  section  of  earth  structure 
(from  Walper,  1976). 


a.  Divergent  or  spreading 
boundary 


TRtNCH 


b.  Convergent  or  collision 
boundary 


c.  Transform  or  transcursion 
boundary 


The  three  kinds  of  boundaries  between  plates 
i f rom  Walper,  1976) . 


Figure  1-2 
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Figure  1-3.  Fault  terminology  (from  Krinitzsky,  1974) 


HANGING  WALL 


Figure  1-4.  Fault  measurement  (from  Krinitzsky,  1974) 
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land forms  developed  along  recently  active  strike-slip 
faults  falter  U . S . Geological  Survey,  1975). 
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Damage  associated  with  movement  along  strike-slip 
fault  (after  Cluff,  Slemmons,  and  Waggoner,  1970) 


Figure  1-7 


BEFORE 


k.  DURING 


Damage  from  normal  fault;  displacements  induced  in 
the  downthrown  block  at  a distance  from  the  fault 
trace  (after  Cluff,  Slemmons,  and  Waggoner,  1970). 


b.  During 


C.  AFTER 


Damage  from  thrust  (reverse)  fault;  earth 
displacements  induced  on  the  upthrown  block 
at  a distance  from  the  fault  trace  (after 
Cluff,  Slemmons,  and  Waggoner,  1970). 


CHAPTER  2 - DIMENSIONAL  RELATIONSHIP  BETWEEN  FAULTS  AND  EARTHQUAKE 
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Earthquake  Magnitude 

Engineers  may  define  a design  earthquake  for  a site  in  terms  of  the 
earthquake  magnitude  M and  the  strength  of  ground  motion.  Factors  which 
influence  the  selection  of  a design  earthquake  are  the  length  of  geo- 
logic fault  structures,  the  relationship  of  the  fault  to  the  regional 
tectonic  structure,  the  geologic  history  of  displacement  along  the 
structure  and  the  seismic  history  of  the  region. 

The  design  earthquake  in  engineering  terms  is  a specification  of 
levels  of  ground  motion  that  the  project  is  required  to  survive  success- 
fully with  an  avoidance  of  loss  of  life  and  acceptable  damage  and  loss 
of  service.  A design  earthquake  on  a statistical  basis  considers  the 
probability  of  the  recurrence  of  a historical  event. 

Earthquake  magnitudes  may  be  specified  in  terms  of  a design  level 
earthquake  which  can  reasonably  be  expected  to  occur  during  the  life  of 
the  structure.  As  such,  this  represents  a service  load  which  the  struc- 
ture must  withstand  without  significant  structural  damage  or  interrup- 
tion of  required  operation.  A second  level  of  earthquake  magnitude  is  a 
maximum  credible  event  for  which  the  structure  must  not  collapse;  how- 
ever, significant  structural  damage  may  occur.  The  inelastic  behavior 
of  the  structure  must  be  limited  to  insure  the  prevention  of  collapse 
and  catastrophic  loss  of  life. 

The  selection  of  a magnitude  level  may  be  based  on: 

1.  Known  design-level  and  maximum-credible  earthquake  magnitudes 
associated  with  a fault  whose  seismicity  has  been  estimated 


2.  Specification  of  probability  of  occurrence  for  a given  life  of 
the  structure  (such  as  having  a 10%  chance  of  being  exceeded  in  25 
years,  NAVFAC  criteria) 

3.  Specification  of  required  level  of  ground  motion  as  in  a code 
provision 


4.  Fault  length  empirical  relationships 


Earthquake  magnitude  can  be  related  to  length  of  fault  for  shallow  depth 
earthquakes.  Data  have  been  plotted  by  Seed  et  al.  (1969),  Krinitzsky 
(1974),  Housner  (1965),  and  Tocher  (1958)  to  provide  the  curves  indi- 
cated (Figure  2-1).  It  is  important  to  note  that  in  some  regions, 
correlations  of  these  types  are  of  little  value  since  many  of  the  impor- 
tant geologic  features  ma  ’ be  deeply  buried  by  weathered  materials. 

Magnitude  as  measured  on  the  Richter  scale  is  calculated  from  a 
standard  earthquake,  one  which  provides  a maximum  trace  amplitude  of  1 
pm  on  a standard  Wood-Anderson  torsion  seismograph  at  a distance  of  100 
km.  Magnitude  is  the  log^  of  the  ratio  of  the  amplitude  of  any  earth- 
quake at  the  standard  distance  to  that  of  the  standard  earthquake.  Each 
full  numeral  step  in  the  scale  (two  to  three,  for  example)  represents  an 
energy  increase  of  about  32  times.  Experience  with  past  earthquakes  is 
presently  the  only  useful  basis  for  relating  fault  length  and  motion  to 
magnitudes  of  associated  earthquakes. 

Fault  Length  and  Earthquake  Magnitude 

A useful  insight  into  the  relationship  between  earthquake  magnitude 
and  length  of  observed  fault  slippage  was  presented  by  Iida  (1965).  He 
grouped  faults  of  all  types  and  obtained  the  plot  shown  in  Figure  2-2. 
He  presents  a spread  of  points  for  which  upper  and  lower  boundaries  are 
drawn.  iida's  spread  of  values  should  be  kept  in  mind  when  one  consid- 
ers the  linear  relationships  that  have  been  suggested  by  numerous 
authors. 
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Figure  2-1  shows  several  such  relationships  that  were  compiled  by 
the  Bureau  of  Reclamation  (1972)  and  reported  by  Krinitzsky  (1974). 
They  differ  from  each  other  only  to  a minor  extent.  The  curves  appear 
to  be  least  reliable  for  low  earthquake  magnitudes.  Krinitzsky  (1974) 
is  quoted  below: 

(Figure  2-3)  illustrates  the  concept  that  a small 
earthquake  may  radiate  displacements  around  a focus  (that 
point  which  is  the  center  of  the  earthquake  and  the  origin  of 
its  elastic  waves)  without  breaking  the  surface  in  any  consis- 
tent manner.  In  fact  it  may  not  break  the  surface  at  all. 

The  energy  release  can  be  accommodated  largely  or  totally  in 
the  subsurface. 

The  spread  of  related  data  is  explored  further 
in  ***  (Figure  2-4).  The  dashed  lines  were  drawn  to  encompass 
a spread  of  worldwide  earthquake  events  gathered  by  Ambraseys 
and  Tchalenko  11968]  ***.  In  addition,  a selection  of  large 
earthquakes  generated  specifically  by  thrust  faults  is  shown 
from  U.  S.  Geological  Survey  (USGS)  data.  From  the  Ambraseys 
and  Tchalenko  (1968)  data  there  is  a lower  limit,  or  cutoff, 
at  magnitude  5.4  below  which  surface  rupture  has  no  signifi- 
cance in  relation  to  magnitude.  Fault  movements  below  magni- 
tude 5.4,  with  very  few  exceptions,  are  contained  in  the 
subsurface  as  shown  in  ***  (Figure  2-3).  The  corollary  to 
this  is  that  if  any  active  surface  fault  movement  is  found,  no 
matter  how  short  it  is,  it  might  be  evidence  for  an  earthquake 
capability  of  greater  than  magnitude  5.4. 

Thrust  faults  ***  present  a more  difficult  problem. 

They  are  capable  of  huge  irregularities  between  magnitude  and 
surface  rupture  even  for  large  earthquakes  (magnitudes  greater 
than  7).  In  fact  their  greatest  irregularity  may  be  in  the 
largest  events.  This  may  result  from  energy  release  by  break- 
age along  multiple  parallel  planes,  much  more  so  than  is  the 
case  in  slippages  of  other  types  of  faults.  Thus,  where 
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magnitudes  are  to  be  estimated  for  thrust  faults,  one  must  be 
prepared  to  apply  values  that  can  represent  strong  earth- 
quakes, though  the  fault  lengths  may  not  be  exceptional. 

Fault  Length  and  Displacement 

Figure  2-5  shows  the  maximum  surface  rupture  on  main  faults  typical 
of  the  United  States  with  superimposed  lines  best  fit  of  world  and  North 
America  data  (Bonilla,  1967  and  1970).  Maximum  surface  displacement  is 
related  to  earthquake  magnitude  in  Figure  2-6  (Bonilla,  1970).  The 
displacements  are  total  displacement  measured  after  the  earthquake  and 
thus  may  have  a large  amount  of  creep  displacement  included.  This  data 
is  limited  and  should  be  used  only  as  a guide  to  minimums. 

Figure  2-7  shows  minor  branch  and  secondary  faults  associated  with 
a major  fault.  Bonilla  (1967)  relates  displacement  on  a minor  fault  to 
the  major  fault.  Figure  2-8  shows  an  envelope  of  the  data. 

The  broad  scatter  shown  by  the  data  is  partly  due  to  (1)  the  vari- 
ety of  field  conditions  in  the  areas  where  the  observations  of  faulting 
were  made,  (2)  the  inadequacy  of  the  length  of  surface  fault  rupture  as 
a measure  of  the  length  of  faulting  at  depth,  and  (3)  theoretical  limi- 
tations . 

Bonilla  and  Buchanan  (1970)  combined  the  data,  "the  study  of  the 
relationships  between  surface  fault  displacement,  surface  fault  rupture, 
and  earthquake  magnitude"  to  produce  Figure  2-9.  "Again,  the  data  have 
very  large  spreads."  Boundaries  were  added  by  Krinitzsky  (1974).  More 
recent  work  by  Mark  and  Bonilla  (1977)  reevaluated  the  data  taking  fault 
type  into  account.  The  results  are  shown  in  Figures  2-10,  2-11,  and 
2-12. 

It  is  important  to  note  the  spread  in  the  data.  Krinitzsky  (1974) 
concludes  "Any  fault  break  in  competent  rock,  no  matter  how  small, 
should  be  taken  as  indicative  of  the  capability  for  at  least  a 5.4 
earthquake."  It  is  important  that  consideration  be  given  to  the  local 
seismic  history  and  the  behavior  of  other  analogous  faults. 
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Absence  of  Surface  Faults 


Consideration  should  be  given  to  the  possibility  of  not  identifying 
all  the  faults  in  a region  which  may  be  active.  This  is  especially  true 
in  central  and  eastern  United  States.  To  account  for  this  a "floating 
earthquake"  (one  that  may  be  assumed  capable  of  occurring  anywhere  in 
the  r^feion)  should  be  considered  (Krinitzsky,  1974). 


LENGTH  OF  FAULT  SLIPPAGE 


SURFACE  RUPTURE  (KM) 


Figure  2-3.  Concept  of  a spreading  fault  novement  that  occurs 
during  an  earthquake  (from  Krinitzsky,  1974). 


Figure  2-4.  Spread  in  Magnitude  versus  surface  rupture  in 
data  assembled  by  Aabraseys  and  Tchalenko  plus 
data  on  selected  thrust  faults  from  the  U.S. 
Geological  Survey  (froa  Krinitzsky,  1974). 


EARTHQUAKE  MAGNITUDE 

Figure  2-5.  Length  of  surface  rupture  on  main  fault  as  related 
to  earthquake  magnitude  (from  Bonilla  and  Buchanan 
1970;  boundary  of  applicability  added). 


EARTHQUAKE  MAGNITUDE 


Figure  2-6.  Maximum  surface  displacement  on  main  fault  as 

related  to  earthquake  magnitude  (from  Bonilla  and 
Buchanan,  1970;  boundary  of  applicability  added). 
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Categories  of  faults 
(from  Bonilla,  1967). 


DISTANCE,  KILOMETERS 


0 I 2 3 4 S « 


DISTANCE  FROM  MAIN  FAULT,  MILES 

Figure  2-8.  Displacement  on  associated  faults  as  percentage  of 
displacement  on  main  fault  (from  Bonilla,  1967). 
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EARTHQUAKE  MAGNITUDE 

Figure  2-9.  Relation  of  earthquake  Magnitude  to  aurface 
•ngth  times  square  of  surface  displacement 
i rom  Bonilla  and  Buchanan,  1970;  boundary 
of  applicability  added). 
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EARTHQUAKE  MAGNITUDE 

Normal  oblique-alip  fault  data  (from  Mark 
and  Bonilla  1977). 
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CHAPTER  3 - FAULT  MOVEMENT 


Prediction  of  Fault  Movement 

Prediction  of  movement  along  specific  fault  zones  is  being  studied 
by  many  researchers  throughout  the  world.  Their  chief  approaches  thus 
far  are  given  in  the  following  paragraphs  taken  from  Krinitzsky  (1974) 
and  quoted  below.  11 

Field  Strains 

Where  faults  move,  they  are  the  visible  displace- 
ments that  occur  within  much  larger  areas  in  which  crustal 
strain  takes  place.  It  is  known  that  the  Niigata  earthquake 
was  preceded  by  vertical  movements  that  went  on  at  a steady 
rate  for  half  a century.  Fifteen  years  before  the  earthquake 
the  movement  became  more  rapid  for  about  10  yr.  Then  there 
were  5 yr  of  very  little  movement  that  ended  with  the  earth- 
quake. Studies  of  elastic  strain  accumulation  may  prove  to  be 
a guide  to  predicting  fault  movements  *** . 

Dil ation 

Fault  zones,  or  the  areas  closely  adjacent  to  major 
faults,  may  dilate  measurably  prior  to  relief  of  stresses  by 
slippage.  Monitoring  the  progress  of  dilation  by  [level 
surveys)  adjacent  to  specific  faults  may  identify  critical 
stages  that  preceded  movement. 

Seismologists  in  the  USSR  found  premonitory  changes 
in  the  ratio  V^/Vg  of  microearthquakes  prior  to  moderate  size 
earthquakes  in  Central  Asia  ***  (where)  ***  is  the  congres- 
sional velocity  and  is  the  shear  velocity.  ***  Corrobora- 
tive results  were  reported  in  the  Adirondack  Mountains  of  New 


York  and  similar  anomalies  were  found  to  have  preceded  the  San 
Fernando,  California,  earthquake  of  February  1971.  Tectonic 
strain  produces  a steady  increase  in  stress  until  dilatancy 
occurs.  ***  the  forming  of  many  small  cracks  ***  takes  place 
at  lesser  stress  levels  than  those  necessary  ***  (for  rup- 
ture]. Initially,  there  is  an  increase  in  pore  volume  without 
a corresponding  increase  in  pore  fluid  and  drops.  is 
not  affected  by  pore  fluids  so  that  Vp/Vg  also  decreases.  As 
water  flows  into  the  new  cracks,  the  ratio  will  return  to 
normal  and  pore  pressure  will  build  up.  With  continuing 
tectonic  stress  and  the  buildup  of  pore  pressure,  an  earth- 
quake will  result. 

Creep 


Similarly,  creep,  or  slow  movements  of  faults 
without  earthquakes,  may  indicate  the  gradual  relief  of 
stress  ***.  (It  may  indicate]  the  partial  relief  of  a locked 
condition  that  will  later  break  and  move  abruptly.  Clear 
identification  of  the  significance  of  creep  is  not  yet  possi- 
ble, but  creep  may  in  time  be  related  to  earthquake  potential. 
Microearthquakes 

Numerous  small  earthquakes  along  active  fault  zones 
may  indicate  the  approach  of  a large  event,  but  this  is  an 
(inconclusive]  relationship.  Of  use,  perhaps,  is  the  plotting 
of  microearthquakes.  They  may  show  a geographical  spread 
that,  correlated  with  historic  earthquakes,  can  suggest  the 
most  probable  locations  for  future  events. 


They  may,  however,  indicate  only  aftershocks  of  a previous  large  event . 


Stress-  and  Strain-related  Effects 

The  physical  properties  of  rocks  are  affected  by 
the  levels  of  stress  and  strain  they  are  undergoing.  Thus 
electrical  resistivity,  magnetism,  piezoelectricity,  the 
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transmission  of  seismic  waves,  the  emission  of  radon,  and 
possibly  other  factors  [may  show  changes]  as  rocks  approach 
the  stage  in  which  fault  slippage  is  about  to  happen. 

Geologically  Determined  Slip  Rates 

The  offset  of  distinctive  rock  units  establishes  the  rate  of  fault 
movement  within  fairly  wide  bounds.  Commonly  these  offsets  average  the 
rate  of  movement  over  millions  of  years,  and  sudden  slip  cannot  be 
distinguished  from  creep.  Data  for  the  San  Andreas  fault  suggest  an 
average  slip  rate  of  1-2  cm/yr  (0.4-0. 8 in./yr)  over  the  last  20  million 
years.  But  to  predict  movements  in  the  immediate  future,  the  most 
recent  few  hundreds  to  thousands  of  years  are  the  most  important  (Wesson 
et  al.,  1975). 

The  history  and  rate  of  fault  movement  have  been  obtained  within 
this  brief  time  period  in  a few  special  circumstances  in  southern 
California  using  absolute  age  dating  techniques. 

Figure  3-1  is  a simplified  sketch  of  a trench  wall  showing  vertical 
deformation  of  initially  flat-lying  sediments  and  sedimentary  contacts 
associated  with  predominantly  horizontal  movement  on  the  Coyote  Creek 
fault  of  southern  California  (from  Clark  et  al.,  1972).  The  trench,  dug 
shortly  after  the  1968  Borrego  Mountain  earthquake,  crosses  a branching 
break  of  the  fault  zone  along  which  about  50  mm  (2  in.)  of  vertical 
displacement  and  about  the  same  amount  of  horizontal  displacement  took 
place  during  the  earthquake.  Deposits  at  points  A,  B,  and  C were  dated 
radiometrically . The  vertical  displacement  of  the  sedimentary  contacts 
plotted  against  the  age  of  the  corresponding  deposits  yields  an  average 
rate  of  vertical  deformation  of  about  0.5  mm/yr  (0.02  in./yr)  for  the 
past  3,000  years.  This  suggests  a recurrence  interval  for  earthquakes 
the  size  of  the  1968  event  of  about  200  years  (Clark  et  al.,  1972) 
(Figure  3-2). 
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Recurrence  Intervals  From  Geologic  Slip  Rates 


II 

SI 


Wallace  (1970)  presented  an  approach  which  has  been  used  by  Lamar 
et  al.  (1973).  The  recurrence  interval  at  a point  can  be  estimated  by 


where  D 
S 
C 


displacement 
long-term  strain  rate 
tectonic  creep  rate 


(3-1) 


Lamar  et  al.  present  the  following  quoted  discussion. 

The  following  assumptions  are  made:  (1)  Slip  on  faults  occurs 
incrementally  during  earthquakes  and  will  continue  at  the  same 
rate  as  that  determined  from  geodetic  data  and  offset  of 
geologic  units.  (2)  Elastic  strain  accumulates  between  earth- 
quakes; the  displacement  during  an  earthquake  represents  the 
release  of  this  accumulated  elastic  strain.  (3)  Tectonic 
creep  is  aseismic  slip  which  reduces  the  accumulation  of 
elastic  strain  available  for  release  during  earthquakes. 

***  Recurrence  intervals  determined  by  Eq.  3-1  represent  a 
long-term  average;  there  is  however,  evidence  of  significant 
local  (Ambraseys,  1970)  and  worldwide  (Davies  and  Brune,  1971) 
time  variations  in  the  level  of  seismic  activity. 

For  most  faults,  creep  can  be  evaluated.  Therefore,  as  an  expedi- 
ent, Equation  3-1  is  simplified  as 


(3-2) 


40 


Equation  3-2  is  thought  appropriate  when  the  rupture  length  is  large 
coaipared  to  the  distance  of  the  site  to  the  fault.  When  this  is  not  the 
case.  Equation  3-2  is  multiplied  by  the  ratio  of  length  of  rupture  to 
total  fault  length  to  account  for  the  spatial  distribution. 

Application  of  Slip  Rate  to  Compute  Recurrence  Data 

Lamar  et  al.  (1973)  have  investigated  the  recurrence  data  for  six 
faults  in  the  southern  California  area.  Their  data  are  summarized  in 
Table  3-1.  The  displacement,  D,  used  in  Equation  3-2  to  calculate 
recurrence  intervals  at  a given  point  on  a fault  were  derived  from 
Bopilla  and  Buchanan  (1970),  except  in  the  case  of  the  White  Wolf  and 
Sierra  Madre  faults.  For  these  two  faults  historic  displacements  were 
available,  and  they  deviated  significantly  from  the  least-squares-fit 
curve  for  reverse  faults. 

Quoting  from  Lamar  et  al.  (1973): 

The  determination  of  long-term  slip  rates  and  recurrence 
intervals  provides  a new  approach  to  earthquake  hazard  assess- 
ment . The  results  can  be  strikingly  different  from  those 
based  on  the  historic  earthquake  record.  For  example,  the  *** 
Garlock  fault,  which  has  not  caused  damaging  historic  move- 
ment , may  have  accumulated  sufficient  elastic  strain  for  a 
destructive  earthquake.  On  the  other  hand,  historic  ruptures 
on  faults  such  as  the  ***  San  Fernando,  have  released  |some] 
accumulated  elastic  strain,  so  that  a destructive  earthquake 
(may  be  less  probable]  for  the  next  few  hundred  years  (depend- 
ing on  the  amount  of  strain  release  and  buildup].  The  recur- 
rence intervals  in  Table  ( 3- 1 ] must  be  considered  tentative 
and  are  subject  to  revision  as  new  information  becomes  avail- 
able. For  the  moat  part,  slip  rates  are  poorly  known,  and  the 


r 
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curves  relating  magnitude  to  surface  displacement  and  rupture 
length  are  based  on  meager  data  with  considerable  scatter. 
More  accurate  age-dates  of  offset  strata  on  faults  are  needed, 
and  additional  studies  following  earthquakes  throughout  the 
world  are  required  to  refine  the  empirical  relations  between 
magnitude,  surface  displacement  and  rupture  length.  This 
research  offers  the  prospect  of  more  quantitative  assessments 
of  earthquake  risk. 
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Table  3-1.  Recurrence  Intervals  on  Selected  Faults  in  Southern  California  [froa  Laaar  et  al., 
edited  by  Southern  California  Section  of  the  Association  of  Engineering  Geologists] 
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Most  likely,  based  on  historic  record. 


IN  METRES 


CHAPTER  4 - SE1SM0L0GICAL  STUDIES  OF  FAULTS 


Introduction 

The  National  Oceanic  and  Atmospheric  Administration's  (NOAA) 
National  Geophysical  and  Solar-Terrestial  Data  Center  in  Boulder, 
Colorado,  maintains  a file  of  earthquake  data  in  computer  format. 

This  earthquake  data  file  is  a magnetic  tape  file  of  earthquake 
origin  times,  locations,  magnitudes,  and  related  information  on  effects 
that  have  been  compiled  from  many  sources.  It  is  a useful  tool  for  pro- 
viding searches  of  the  seismicity  of  selected  regions,  for  preparing 
lists  of  earthquakes  of  given  characteristics  (large  magnitude,  inten- 
sity, tsunami,  etc.),  for  plotting  maps,  and  for  making  statistical 
studies.  It  is  updated  as  new  earthquake  data  become  available  in 
machine-readable  form. 

The  earthquake  data  file  (sometimes  called  hypocenter  data  file) 
contains  information  on  approximately  116,000  earthquakes  for  the  period 
June  11,  1638,  through  December  31,  1974.  The  file  includes  for  each 
event  the  date,  origin  time,  geographic  location,  focal  depth,  magni- 
tude, and  intensity  when  available.  It  contains  additional  information 
on  the  source  of  the  data,  quality  factors,  associated  phenomena 
(tsunami,  volcanism,  faulting,  etc.),  and  seismic  region  number.  In 
addition  to  earthquakes,  the  file  contains  about  500  phenomena  that  have 
surface  effects  similar  to  earthquakes,  including  known  or  suspected 
explosions  and  associated  collapse  phenomena,  coal  bumps,  rockbursts, 
quarry  blasts,  and  other  earth  disturbances.  At  present,  data  for  the 
1638-1899  period  are  practically  limited  to  earthquakes  in  the  United 
States,  a few  hundred  in  number,  which  had  a maximum  intensity  of  5 or 
greater.  A number  of  data  fields  for  some  events  are  unfilled  because 
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the  information  is  not  available.  Information  on  cultural  effects, 
intensity,  and  other  phenomena  associated  with  the  event  has  been  in- 
cluded for  earthquakes  in  the  United  States.  This  information  has 
sometimes  been  entered  for  non-United  States  earthquakes,  particularly 
since  May  1968,  although  significant  gaps  still  exist. 

The  quality  of  epicenter  determinations  varies  significantly  with 
the  time  period  studied.  Before  1900  locations  are  usually  noninstru- 
mentally  determined  and  are  given  as  the  center  of  the  macroseismic 
effects.  Most  instrumental  epicenters  prior  to  1961,  excluding  local 
earthquakes  in  California,  were  located  to  the  nearest  1/4  or  1/2  degree 
of  latitude  and  longitude.  Reliable  information  on  the  quality  of  many 
epicenter  determinations  is  lacking.  Beginning  in  1960,  epicenters  have 
been  determined  by  computer,  and  the  accuracy  is  generally  better. 
However,  although  stated  to  tenths  or  hundredths  of  a degree,  the  loca- 
tion accuracy  is  usually  a few  tenths  of  a degree.  Since  May  1968,  the 
latitude  and  longitude  values  for  most  events  have  been  listed  to  three 
decimal  places.  This  precision  is  not  intended  to  reflect  the  accuracy 
of  the  location  of  events  except  for  local  California  earthquakes  and 
special  epicenter  determinations.  Where  several  sources  have  determined 
an  epicenter  for  the  same  earthquake,  one  solution  has  been  designated 
as  the  most  reliable.  Usually  it  is  the  source  believed  to  contain  the 
best  data  set  for  the  earthquake.  In  come  cases,  data  from  two  sources 
were  combined  to  provide  a more  complete  record. 

Magnitudes  from  a number  of  different  sources  are  included  in  the 
earthquake  data  file.  Gutenberg  and  Richter  (1954)  and  Richter  (1958) 
discuss  the  development  of  the  magnitude  scale.  Many  magnitudes  pub- 
lished by  Gutenberg  and  Richter  (1954)  were  later  revised  by  Richter 
(1958).  The  revised  magnitudes  are  used  in  the  file  even  though  the 
source  is  identified  as  Gutenberg  and  Richter  (1954).  The  concept  of 
earthquake  magnitude  is  not  restricted  to  one  value.  Several  defini- 
tions are  possible,  depending  on  which  seismic  waves  are  measured. 
Three  different  magnitude  scales,  BODY  WAVE  (MB),  SURFACE  WAVE  (MS),  and 
LOCAL  (ML),  are  distinguished  in  this  file.  In  addition  another  data 
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field,  OTHER  MAGNITUDE,  has  been  included  when  it  was  unclear  which 
scale  was  used.  Richter  (1958)  and  other  modern  seismology  references 
provide  detailed  discussions  of  this  topic.  The  different  scales  do  not 
give  exactly  comparable  results,  and  different  values  frequently  are 
given  for  the  same  earthquake.  It  is  common  practice  to  average  the 
individual  magnitudes  from  different  stations  to  get  a more  uniform 
value  within  each  scale  (MB,  MS,  and  to  a lesser  extent  ML). 

In  general,  the  file  contains  earthquakes  of  magnitude  4.0  or  less 
only  for  the  United  States  region  and  for  areas  within  dense  seismic 
station  networks.  However,  no  claim  is  made  for  the  statistical  homo- 
geneity of  these  events.  Inclusion  of  earthquakes  of  magnitude  4.0  to 
5.0  also  is  influenced  by  the  proximity  of  seismic  stations  to  the 
source  or  epicenter. 

A maximum  intensity  is  listed  for  many  of  the  earthquakes.  Each  is 
assigned  according  to  the  Modified  Mercalli  Intensity  Scale  of  1931. 
Some  of  these  values  have  been  converted  from  reported  intensities  on 
other  scales. 

Correction  of  Epicenters 

Early  records  of  earthquakes  may  show  inaccurate  epicentral  loca- 
tions. Also  magnitudes  or  intensities  may  differ  from  values  that  would 
now  be  assigned.  As  an  example,  Krinitzsky  (1974)  reports  the 
following: 

The  location  of  an  NOAA  listed  earthquake  of 
15  May  1909  in  southern  Saskatchewan  is  shown  in  (Figure  4-1). 

It  has  a maximum  intensity  of  VIII  and  is  the  largest  recorded 
event  in  this  general  area.  Its  location  was  made  when  seis- 
mograph stations  were  few,  and  those  that  operated  were  far 
less  accurate  than  they  are  today.  The  intensity  of  VIII  was 
assigned  on  the  basis  of  the  large  felt  area. 
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The  event  was  checked  by  referring  to  the  newspaper 
accounts  of  this  tine.  This  was  not  too  formidable  a task. 

State  archives,  state  university  collections,  and  national 
archives  often  have  collections  of  local  newspapers  gathered 
in  central  depositories.  One  may  request  microfilm  copies  of 
these  papers  for  the  dates  of  interest,  review  them,  and 
assign  intensity  values  based  on  their  descriptions.  This 
exercise  provided  intensities  for  more  than  50  communities 
though  this  was,  and  still  is,  a thinly  populated  area.  The 
resulting  iso-intensity  map  shows  there  was  no  intensity  VIII. 

The  greatest  was  VI.  Also  the  region  of  VI  had  its  center  to 
the  east  by  about  1 deg  from  the  NOAA  location.  In  (Figure 
4-1 J the  revised  location  is  shown  in  relation  to  three  other 
earthquakes  of  1956,  1968,  and  1972.  Their  locations  were 
accurate  to  begin  with  because  of  better  instrumental  capabil- 
ity. One  is  associated  with  a seismically  interpreted  fault 
that  also  agrees  with  a geologically  mapped  fault.  Its  trend 
is  toward  the  three  other  events.  Thus  the  revised  location 
for  the  1909  event  helps  to  interpret  a fault  trend. 

Seismic  Arrays 

Seismometers  have  been  installed  near  known  active  faults  to  record 
microearthquakes.  Figures  4-2  and  4-3  show  the  location  of  recorded 
microearthquakes  for  a year  and  major  fault  zones  in  California.  The 
events  recorded  ranged  in  magnitude  from  0.5  to  1.5.  The  figures  show 
areas  where  the  microearthquakes  closely  trace  the  faults.  There  are 
also  regions  where  few  microearthquakes  occurred.  The  San  Andreas  fault 
north  of  the  Sargent  fault  exhibits  less  activity  than  southern  por- 
tions, perhaps  indicating  it  is  locked.  However,  sufficient  micro- 
earthquakes have  occurred  to  show  the  continuation  of  the  fault.  It  may 
be  reasoned  that  small  failures  might  occur  before  a major  rupture 
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occurs;  alternatively,  the  large  number  of  microearthquakes  demonstrate 
active  creep  which  may  be  sufficient  to  prevent  sizable  strain  accumula- 
tion and  preclude  a large  event. 

Since  earthquakes  are  associated  with  faults,  it  might  be  thought 
that  they  should  precisely  overlay  the  fault  location.  This  is  not  the 
case  because  the  distribution  of  seismometers  is  uneven  and  has  changed 
with  time.  There  are  limitations  in  the  accuracy  of  the  techniques  used 
to  locate  epicenters,  principally  from  variations  in  the  assumed  propa- 
gation velocities.  Further  explanation  for  the  location  of  epicenters 
being  off  their  associated  fault  comes  from  the  simplified  model  used  to 
locate  them.  This  is  explained  in  Figure  4-4.  Note  that  the  center  of 
earthquake  energy  is  located  at  the  focus-  For  an  inclined  fault,  the 
surface  location  (epicenter)  is  a distance  removed  from  the  surface 
fault  location.  It  is  only  in  vertical  faults  that  one  might  expect  the 
epicenter  to  lie  on  the  fault. 

Krinitzsky  (1974)  concludes  that  earthquakes  can  be  related  to 
existing  faults  and  that  the  possibility  of  formation  of  new  faults 
should  not  be  considered  in  design.  Large  earthquakes  require  fault 
breaks  of  considerable  distance.  The  uncertainties  that  occur  in  the 
association  of  earthquakes  with  faults  can  occur  only  for  small  events. 
Generally  in  the  western  United  States  the  extent  of  geologic  investiga- 
tion precludes  the  omission  of  a large  fault  remaining  unknown.  How- 
ever, there  are  uncertainties  associated  with  eastern  earthquakes.  For 
example,  causative  faults  responsible  for  the  New  Madrid  earthquakes  of 
1811  and  1812  have  not  yet  been  identified.  This  may  be  the  result  of 
insufficient  geologic  investigations.  The  importance  of  considering  the 
extent  and  quality  of  geologic  investigations  is  evident. 

Limitations  to  Historic  Data 

A period  of  demonstrated  quiescence  over  a geological  time  period 
indicates  inactivity  of  the  fault  and  probable  continued  inactivity. 
However,  inactivity  over  a period  of  historic  recording  (50  to  100 
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years)  does  not  isiply  future  inactivity.  Rather  it  may  point  to  a 

region  which  is  locked  and  through  which  a major  fault  rupture  may 
propagate.  The  last  two  earthquakes  producing  damage  in  southern 
California  (Arvin-Tehachapi , 1952,  and  San  Fernando,  1971)  occurred  on 
faults  lacking  historic  activity.  With  the  exception  of  the  San  Jacinto 
fault  system,  every  known  event  greater  than  magnitude  6 in  southern 
California  occurred  on  a fault  without  prior  historic  activity.  Caution 
must  be  exercised  to  recognize  that  the  accuracy  of  an  incomplete  data 
base  is  very  limited  when  extrapolated  for  return  periods  greatly  ex- 
ceeding the  length  of  the  period  of  recorded  data.  Furthermore,  after- 
shocks must  be  distinguished  from  main  shocks.  An  area  having  recently 
undergone  a large  event  releasing  strain  built  up  for  hundreds  or  thou- 
sands of  years  is  probably  safe  against  that  size  release  in  the  near 
future.  Thus  a recent  large  event  on  a fault  might  actually  indicate 
safety  in  the  immediate  future  rather  than  an  indication  of  activity.  A 
single  event  by  itself  cannot  give  an  accurate  measure  of  return  time. 
Krinitzsky  (1974)  is  quoted  below: 

In  the  United  States  the  history  of  earthquake 
activity  is  greatly  truncated.  At  best  the  record  covers  less 
than  350  years.  This  it  does  in  very  few  places.  For  most  of 
the  country  it  is  about  150  years.  This,  however,  is  general 
earthquake  history  and  not  the  record  of  movement  on  specific 
faults  or  in  fault  zones.  In  other  parts  of  the  world  the 
record  for  both  earthquake  history  and  interpreted  fault 
activity  is  better. 

N.  N.  Ambraseys  ***  (1971)  has  studied  the  record 
of  damaging  earthquakes  for  the  past  two  millennia  in  some 
portions  of  the  Near  East.  [Figure  4-5)  ***  shows  the  cumula- 
tive distribution  during  this  time  of  damaging  earthquakes 
along  (a)  the  Border  Zone,  which  is  a northeast-southwest 
trend  of  faults  that  extends  toward  the  Dead  Sea  along  the 
border  of  Israel;  and  (b)  the  Anatolian  Zone  that  is  roughly 
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Zone.  Ambraseys  points  out  that  (luring  the  first  five  centur- 
ies the  Border  Zone  was  quiescent  while  the  Anatolian  Zone  was 
active.  For  the  following  six  centuries  the  pattern  was 
reversed  only  to  he  reversed  again  during  the  eleventh  centu- 
ry. In  this  case  there  is  a factor  of  dependence  of  activity 
along  one  /.one  on  that  of  the  other.  The  two  join  to  encom- 
pass a miniplate,  and  movement  tends  to  shift  from  one  side  to 
another  on  a cyclical  basis.  This  tends  to  discredit  the 
validity  of  probabilistic  projections  based  only  on  a short 
history. 

Also,  as  noted  in  Krinitzsky  (1974)  there  are  likely  to  he: 

***  variations  in  rates  of  slippage  along  various  segments  ol 
any  one  long  active  fault  or  fault  zone.  This  is  known  to  he 
the  case  along  the  San  Andreas.  During  its  relatively  short 
historic  period,  it  was  noted  that  major  earthquakes  moved 
from  place  to  place  along  the  fault.  Portions  that  once 
moved,  notably  the  segment  that  slipped  during  the  San 
Francisco  earthquake  of  1906,  have  become  locked  while  slip- 
pages occurred  elsewhere.  |Figure  4-6|  is  a schematic  state- 
ment of  variations  in  rates  of  slippage  in  inches  per  year,  as 
given  by  Brown  and  Wallace  (1968).  The  slippage  is  that  which 
has  occurred  during  an  interval  of  60  years.  Though  this  is 
instructive  of  irregularities  in  the  rates  of  movement,  it  is 
not  intended  as  a guide  to  the  future.  If  anything,  the 
future  is  very  likely  to  he  different.  The  segment  between 
('ho  lame  and  Camp  Dix  appears  to  he  in  a locked  position. 
Stresses  are  building  up.  One  day  this  segment  will  rupture 
suddenly. 
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Figure  4-1 


Figure  4-4.  Epicenter  of  earthquakes  (from  U.S.  Geological 
Survey,  1975). 
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Figure  4-5.  Time  distribution  of  damaging  earthquakes  in 
Turkey  (from  Ambraseys  in  Nature,  1971). 
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Figure  4-6.  Variation  in  rate  of  slippage  along  the  San  Andreas 
fault  (from  Brown  and  Wallace,  1968). 


CHAPTER  5 - STOCHASTIC  APPROACH  TO  EARTHQUAKE  DATA 


Earthquake  Models 

Oliveira  (1975)  summarizes  as  follows: 

The  generation  of  earthquakes  as  related  to  the  geotec- 
tonic  mechanisms  of  faults  and  the  timing,  intensity  and 
location  of  future  earthquakes,  cannot  be  predicted  with 
certainty  in  a deterministic  way.  ***  The  generation  of  earth- 
quakes in  space  and  time  falls  under  the  general  category  of 
stochastic  processes  ***  [Cox  and  Miller,  1970;  Lomnitz, 
1974). 

[Time  Distribution.)  Both  the  Markov  and  Poisson  repre- 
sentations are  used  to  simulate  the  occurrence  of  earthquakes 
generated  in  time.  These  representations  are  stochastic  point 
processes.  The  Markov  model  differs  from  the  Poisson  model  in 
that  the  occurrences  of  new  events  depend  on  past  events 
whereas  in  the  latter  model  these  occurrences  are  independent 
of  past  events.  Results  obtained  using  these  two  models 
differ  somewhat.  The  Markov  model,  better  adjusted  to  the 
elastic  rebound  theory***  [Reid  (1911)1,  has  some  drawbacks 
due  to  difficulty  in  setting  the  initial  conditions  and  re- 
quires more  numerical  treatment.  The  Poisson  model,  on  the 
other  hand,  does  not  always  agree  with  experimental  data  for 
small  magnitude  earthquakes,  because  it  ignores  the  tendency 
of  earthquakes  to  cluster  in  space  and  time  ***  [Knopoff, 
1964).  Nevertheless,  the  Poisson  model  seems  to  give  adequate 
results  in  many  situations  except  for  cases  of  very  large 


return  periods  and  cases  when  analysing  structures  that  have 
just  been  damaged  [Rosenblueth,  1973).  Emphasis  is  given  to 
the  Poisson  model  in  the  present  treatment. 

******* 


P„(n,Xt)  = C X^t?n,  K5-D) 

where  p^  = the  probabilty  of  occurrence  of  n events  of  a given 
magnitude  range 

X = the  mean  rate  of  occurences  per  unit  of  time  t 

The  probability  of  the  amount  of  time  between  two  consecutive 
events  is  given  by  the  exponential  law 

PT(t)  = X e‘Xt  05t5»  1(5-2)) 

Cornel  (1974),  Esteva  (1968),  and  Milne  and  Davenport  (1969)  have, 
at  least  for  main  shocks,  verified  the  Poisson  model  favorably;  however, 
Knopoff  (1964)  reports  after  testing  different  regions  in  the  world  the 
Poisson  model  is  inadequate  to  explain  the  time  distribution  of  low 
magnitude  shocks.  Small  magnitude  earthquakes  may  contribute  heavily  to 
the  general  seismic  risk,  rendering  the  Poisson  model  inadequate.  A 
more  refined  model  of  generation  should  then  be  used.  In  seismic  risk 
analysis  the  Poisson  model  is  generally  an  acceptable  approach  except 
when  working  with  return  periods  of  several  thousand  years  or  in  analy- 
zing one  structure  that  has  just  been  damaged  (Oliveira,  1975).  The 
Poisson  distribution  is  intended  to  be  used  with  National  Bureau  of 
Standards  Seismicity  Tables  (Culver  et  al.,  1975).  The  main  deficiency 
of  the  simple  Poisson  model  is  that  it  ignores  the  tendency  of  earth- 
quakes to  come  in  groups  which  are  often  triggered  by  a large  main 
shock;  magnitude  is  treated  as  a separate  independent  phenomenon. 
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Figure  5-1  gives  the  probability  of  not  exceeding  specified  values 
during  the  life  of  a facility. 

An  improvement  over  the  Poisson  distribution  is  the  Weinbull  dis- 
tribution suggested  by  Chou  and  Fisher  (1975) 


where  u and  y are  scale  and  shape  parameters,  respectively.  Several 
methods  are  available  to  estimate  the  parameters.  The  maximum- 
likelihood  method  is  recommended  because  it  utilizes  the  available 
information  in  the  most  appropriate  manner.  The  shape  factor  y is 
estimated  by  solving  the  general  equation: 


and  the  scale  parameter  p is  obtained  by  equating 


where  n is  the  size  of  the  sample  and  t^  is  the  time  interval  involved. 

A graphical  method  of  plotting  historic  earthquake  data  is  very 
useful  and  widely  used  in  practice.  A new  random  variable  is  introduced 
as  Z = ln(pt^),  and 


The  list  of  earthquake  occurrence  data  is  arranged  in  groups  of  in- 
tensity or  magnitude  ranges.  Within  each  group  of  earthquake  events  of 
the  same  magnitude  range  the  data  are  ordered  by  time  intervals  between 
occurrences,  the  most  frequent  first.  The  plotting  position  of  any  data 
point  within  a group  is  [t^,  F(i)|  where  t.  is  the  i1*1  longest  time 
interval  and 

F(U  = J-jr  (5-7) 

where  i = position  on  the  list;  i.e.,  first,  second,  etc. 
n = total  number  of  events  in  the  group 


The  parameters  p and  y are  determined  by  the  intercept  and  slope,  re- 
spectively, of  the  plotted  data  following  the  relationship: 


Z.  = In  p + y In  t. 
1 1 


(5-8) 


where  t. 

l 


is  the  i 


th 


longest  time  interval  and 


is  given  by 


lnln 


(• 


n + 1 \ 

- i + 1/ 


(5-9) 


(5-10) 


When  y = 1,  the  Weibull  distribution  is  equivalent  to  the  Poisson  dis- 
tribution. 

Magnitude  Distribution.  The  number  of  earthquakes  of  magnitude  £ m 
to  be  expected  in  a given  interval  of  time  can  be  translated  into  a 
random  variable  using  an  exponential  distribution.  This  distribution  is 
similar  to  a Poisson  distribution;  however,  it  is  continuous  rather  than 
discrete. 
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The  magnitude  distribution  is  exponentially  related  to  magnitude 
through  a parameter  p.  The  probability  density  function  is 


f(m)  = p e 


-p(m-m  ) 


(5-11) 


and  the  cumulative  distribution  function  is 


FM(m) 


= 1 - e 


-p(m-m  ) 


m 2 m 


where  p = 


b log^lO 


1(5-12)] 


n»o  = a constant  defined  by  the  lower  magnitude  limit 

b = the  coefficient  of  Richter's  law  of  magnitudes 

(logjQ  = a - bm).  (See  Site  Recurrence  below). 

Oliveira  (1975)  states  "Experimental  data  show  that  the  linear  form, 
logjO  = a * bm,  of  Richter's  law  does  not  hold  for  very  large 

magnitudes.  This  suggests  the  use  of  truncation  at  a certain  value  of 
magnitude,  m^ , or  the  introduction  of  a quadratic  term  (Merz  and 
Cornell,  1973)."  Table  5-1  summarizes  the  probability  law  of  magnitudes 
for  the  above  three  cases. 

Seismic  Risk 


The  occurrence  and  magnitude  mechanisms  defined  above  characterize 
the  generation  of  earthquakes  in  time  and  intensity.  Randomness  in 
space  is  another  very  important  consideration  in  assessing  seismic  risk. 
Earthquake  activity  might  be  connected  with  the  presence  of  faults. 
Only  geological  and  geotectonic  considerations  associated  with  the  past 
history  of  occurrences  of  earthquakes  can  characterize  the  distribution 
laws  of  epicentral  locations.  Spatial  correlations  among  earthquake 
source  zones  have  been  detected  in  several  seismic  areas  (Ferraes, 
1973).  Sophisticated  models  have  been  developed  to  take  into  account, 
simultaneously,  space  and  time  generation  (Lomnitz,  1974). 


If  each  source  zone  is  considered  to  have  a constant  mean  rate  of 
occurrence  A^,  following  a Poisson  probability  law,  the  set  of  all 
independent  source  zones  will  generate  another  Poisson  process  with  mean 
rate  A = IA^  (the  sources  can  be  reduced  to  point  sources).  This  prin- 
ciple has  been  used  in  studies  of  risk  analysis  at  a site. 

From  the  Poisson  model  of  occurrences  and  the  exponential  type 
distribution  of  magnitudes  it  is  possible  to  show  that  the  extreme 
distribution  of  magnitude  follows  an  extreme  type  I distribution 
(Gumbel,  1958)  of  the  form 

-p(m-m  ) > 

GM(m)  = exp  - A e 0 m ~ “o  (5-13) 

Site  Recurrence 

Historical  events  within  a distance  around  the  site  can  be  tabu- 
lated to  evaluate  the  site  regional  seismicity.  This  gives  an  average 
measure  of  the  earthquake  proneness  of  a geographic  area  and  is  usually 
discussed  in  terms  of  an  earthquake  recurrence  equation  of  the  form 

log^  N = a - b m (5-14) 

where  N = number  of  earthquakes  per  year  exceeding  a magnitude  m 
a = seismicity  constant  for  an  area 
b = constant  for  an  area 

The  tabulated  historical  data  may  be  arranged  according  to  magni- 
tude, normalized  to  1 year  by  dividing  by  the  timespan  covered,  and  the 
number  of  events  greater  than  a magnitude  m plotted  against  m.  This 
graph,  plotted  on  semi log  graph  paper,  can  be  used  to  evaluate  the 
coefficients  a and  b. 

Some  specific  recurrence  plots  for  major  faults  having  historic 
activity  are  available  in  the  literature  (e.g.,  Environmental  Develop- 
ment Agency,  1975).  With  this  data  it  is  possible  to  develop  probabil- 
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ity  of  occurrence  data  for  selection  of  a design  and  collapse-level 
earthquake  magnitude.  The  distance  of  the  site  to  the  fault  is  impor- 
tant in  determining  the  site  motion  as  will  be  shown  later. 

Nuttli  (1974)  has  compiled  recurrence  information  for  the  central 
United  States.  The  recurrence  equation  for  the  region  of  his  study 
covering  the  states  of  Illinois,  Indiana,  Kentucky,  Tennessee, 
Mississippi,  Arkansas,  and  Missouri  is 

log  N = 3.55  - 0.87(±0. ll)m  (5-15) 

where  N is  the  number  of  earthquakes  per  year  occurring  in  the  magnitude 
range  m ± 0.2.  This  indicates  a return  period  of  6.3  years  for  a magni- 
tude 5 earthquake,  47  years  for  a magnitude  6 earthquake,  and  347  years 
for  a magnitude  7 earthquake. 

The  National  Bureau  of  Standards  (Culver  et  al.,  1975)  has  pub- 
lished a guide  for  the  evaluation  of  earthquake  hazards.  This  guide 
contains  tables  which  contain  the  coefficients  a and  b and  an  effective 
site  distance  for  1/2-degree  squares  of  the  entire  United  States.  This 
data  was  developed  considering  the  seismicity  of  each  sector  and  aver- 
ages the  results  of  the  sector.  When  specific  recurrence  data  is  avail- 
able for  a specific  fault  within  200  miles  of  the  site,  the  actual  fault 
data  and  distance  to  the  site  should  be  used. 

A probability  distribution  may  be  applied  to  recurrence  data  to 
determine  the  probability  of  occurrence  of  a specified  magnitude  event 
in  a number  of  years  (the  Poisson  model  is  the  most  common  as  noted 
earlier).  The  cumulative  probability  of  occurrence  of  an  event  £ m in 
time  t is 


P(t)  = 1 - e-A(m)t  (5-16) 

where  A(m)  is  the  rate  of  occurrence  of  an  event  £ m determined  from  the 
recurrence  plots. 
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Table  5-1.  Magnitude  Distribution  Law  (from  Oliveira,  1975) 


01  3oi 


Probability  distribution  of  maximum  for  the  time  life  period 


TIME  LIFE  (years) 


Figure  5-1.  Conversion  factors  as  function  of  probability 
of  not  exceeding  y-values  and  of  life-time 
(from  Oliveira,  1975). 
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Return  period  corresponding  to  the  maximum  yearly  distribution 
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CHAPTER  6 - GROUND  MOTION 


Earthquake-Induced  Ground  Motion  Levels 


The  rupture  of  a fault  generates  body  waves  propagating  within  the 
earth  and  surface  waves  propagating  on  the  surface.  The  body  waves  are 
composed  of  dilational,  longitudinal,  compressional  waves  and  distor- 
tional,  transve-  e,  shear  waves,  both  vertical  and  horizontal.  The 
surface  waves  are  waves  which  involve  a surface  layer  in  horizontal 
transverse  vibration  and  Rayleigh  waves  comprising  vibration  in  retro- 
grade orbit.  Each  of  these  types  of  waves  propagates  at  its  own  velo- 
city and  arrives  at  specific  locations  at  different  times.  Body  waves 
travel  fastest  through  the  higher  velocity  mediums  at  depth  (bedrock). 
The  waves  tend  to  be  refracted  toward  the  vertical  as  they  propagate 
upward  through  increasingly  softer  materials  near  the  surface.  (This 
follows  from  Snell's  law  because  shallower  layers  generally  have  lower 
wave  velocities.) 

Thus,  shear  waves  tend  to  approach  the  surface  traveling  in  the 
vertical  direction  and  vibrating  in  a horizontal  plane.  Most  structures 
are  well-designed  for  vertical  loads  by  use  of  standard  safety  factors. 
However,  horizontal  shear  waves  induce  motions  and  load  in  the  structure 
in  the  horizontal  plane.  Strong  motion  accelerograms  show  that  vertical 
accelerations  are  often  about  two-thirds  as  large  as  horizontal  vibra- 
tions. For  this  reason,  vertically  propagating,  horizontal  motions  are 
considered  of  greatest  importance  and  most  work  has  centered  on  under- 
standing them. 

Rock  motions  beneath  a particular  site  will  depend  upon  the  energy 
released  along  the  fault  (characterized  by  the  magnitude  of  the  earth- 
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quake  ami  type  of  fault  and  the  distance  of  the  site  from  the  zone  of 
energy  release).  Figure  6-1  Nhows  the  factors  affecting  seismic  wave 
transmission. 

Having  the  estimated  magnitude  of  a seismic  event  and  its  distance 
from  the  site  of  interest,  it  is  then  possible  to  predict  ground 
motions . 

In  order  to  estimate  the  response  of  any  point  on  the  surface,  it 
is  necessary  to  determine  the  surface  topography  and  the  underlying  rock 
configurations.  The  types  of  soils  and  their  characteristics  determine 
the  soil  response  under  dynamic  loading.  It  is  important  to  consider 
the  mechanism  of  underlying  soil  deposits  and  the  transfer  mechanisms  to 
the  surface.  (For  earthquake  engineering  purposes,  bedrock  may  be 
loosely  defined  as  a material  exhibiting  a shear  wave  velocity  of  2,500 
ft/sec  or  greater.)  However,  ground  motion  prediction  is  an  inexact 
procedure;  this  will  be  discussed  below. 

Western  United  States 

Several  correlations  are  available  to  relate  earthquake  magnitude, 
distance  from  causative  fault,  and  peak  acceleration.  An  early  correla- 
tion by  Housner  (1965)  estimated  peak  acceleration  for  seismic  motions 
experienced  at  the  ground  surface.  At  that  time  the  influence  of  over- 
lying  soils  on  the  peak  accelerations  could  not  be  evaluated.  Peak 
accelerations  represented  by  these  curves  are  not,  therefore,  directly 
applicable  to  the  bedrock  level. 

A study  by  Schnabel  and  Seed  (1972)  evaluated  the  peak  rock  outcrop 
motion*.  The  curves  given  in  Figure  6-2  resulted  from  this  study  and 
relate  peak  rock  acceleration  with  earthquake  magnitude  and  distance 
from  causative  fault. 

Seed  et  al.  (1975)  show  that  for  the  range  from  20  to  250  km  all 
measured  accelerogram  peak  data  from  6.5  magnitude  earthquakes  tends  to 
lie  within  a band  of  two  standard  deviations  (Figure  6-3).  Prediction 
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at  closer  distances  is  based  on  extrapolation.  Figures  6-4  and  6-5  show 
similar  data  for  stiff  soil  conditions  and  deep  cohesionless  sites. 
Figure  6-6  compares  acceleration  in  rock  to  acceleration  in  various  soil 
prof i les . 

Trifunac  and  Brady  (1975a)  suggest  that  peaks  of  ground  motion  may 

2 

be  scaled  by  the  following,  where  A is  in. /sec 

max 

1 " * 1°K|0|Ao(K)'  * ,°«lo'Ao(M)l 

They  characterize  site  conditions  s as 

0 = all uv i urn 

1 = intermediate 

2 = basement  rock 

Tables  6-1  and  6-2  give  values  of  log^lA  ^R))  and  log^l Aq(M) ) . They 
extend  the  applicability  of  the  equation  to  include  confidence  level  p, 
site  condition  s,  and  component  direction  v (v  = 0 for  horizontal,  v = 1 
for  vertical)  as  follows 

lo*10|Ao(M’»,’s’v)l  = *P  + l»M  + c + ds 

♦ e v ♦ f M2 

The  following  coefficients  in  the  expression  are  results  of  a least - 
squares-fit  to  available  strong  motion  data. 
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= ap  + bM  + c + ds 


+ ev  + fM 


- f(M  - M r M 5 M 

max  max 


1o810*Ao(M,P,S,v)' 


ap  + bM  ♦ c + ds 


+ ev  + fM 


M < M 5 M . 
max  nun 


= ap  + bM  . + c + ds 

r min 


+ ev  + fM  . 

min 


M £ M . 
min 


where 


= -0,898 

= -1.789 

= 6.217 

= 0.060 
= 0.331 


= 0.186 


Total  N Data  = 


The  Trifunac  and  Brady  (1975b)  equation  results  in  ground  motions 
(Tables  6-3  and  6-4)  near  the  fault  much  greater  than  Seed's  estimates 
(Figures  6-2  and  6-6). 

Page  et  al.  (1972)  reviewed  ground  motion  records  where  the  dis- 
tance to  fault  was  accurately  known.  They  found  that  peak  acceleration 
attenuates  with  distance  r at  a rate  in  the  range  r to  r at 
distances  beyond  10  km  for  magnitude  5,  about  20  km  for  magnitude  6,  and 
40  km  for  magnitude  7 (Figure  6-7).  For  distances  less  than  10  km  there 
are  no  strong  motion  data  for  shocks  larger  than  magnitude  6.  They 
estimate  the  near  fault,  free-field,  horizontal  ground  motion  in  Table 
6-5.  The  free-field  peak  acceleration  values  close  to  the  fault  are 
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higher  than  those  estimated  by  Schnabel  and  Seed  (1972).  This  table 
reflects  the  analysis  of  the  6.6-magnitude  1971  San  Fernando  earth- 
quake's Pacoima  Dam  record.  Motions  for  magnitudes  greater  than  this 
were  extrapolated  considering  that  peak  acceleration  increased  with 
magnitude  and  that  near-fault  peak  acceleration  is  proportional  to  the 
effective  stress  causing  slippage. 

Table  6-6  (after  McGuire,  1976)  summarizes  different  attenuation 
equations  developed  for  ground  motions.  Figure  6-8  (after  Donovan, 
1974)  compares  some  of  the  attenuation  equations  for  a magnitude  6.5 
earthquake  with  measured  data  for  the  1971  San  Fernando  earthquake. 
Although  site  conditions  are  not  defined,  other  than  as  soil  or  rock,  it 
is  clear  that  there  is  a range  in  predictions  of  motion  at  a given 
distance.  The  estimation  of  ground  motion  is  still  an  imperfect  proce- 
dure with  upper  and  lower  bounds  differing  as  much  as  100%  from  mean 
predictions  demonstrating  the  statistical  spread  in  the  data.  Since  no 
single  relationship  has  been  shown  to  give  more  reliable  results  than 
any  other,  it  is  recommended  that  the  data  presented  in  this  section  be 
reviewed  to  select  the  methodology  most  compatible  with  the  intended 
usage,  site  characteristics,  and  importance  of  the  structure. 

Central  and  Eastern  United  States 

Insofar  as  its  effects  on  land  form  are  concerned,  the  most  destruc- 
tive earthquakes  in  North  America,  since  its  settlement  by  Europeans, 
occurred  in  the  New  Madrid  fault  zone  between  December  1811  and 
February  1812  (Nuttli  and  Zellweg,  1974;  Nuttli,  1973a,  b,  c,  and  1974). 
This  area  includes  southeast  Missouri,  northeast  Arkansas,  western 
Tennessee,  western  Kentucky  and  southern  Illinois.  However,  in  this 
region  the  occurrence  of  earthquakes  is  very  infrequent.  The  historical 
data  is  limited  to  less  than  200  years.  Thus,  it  is  possible  that  there 
are  regions  in  the  United  States  that  have  not  experienced  a destructive 
earthquake  in  the  past  200  years  but  might  suffer  one  in  the  next  100 


years.  There  is  a significant  lack  in  strong  motion  records  in  the 
eastern  and  central  United  States.  Since  most  of  the  recorded  data  for 
the  United  States  have  come  from  California,  the  characteristics  of  the 
accelerograms  reflect  seismologic  and  geologic  conditions  different  from 
those  of  the  east  and  central  parts  of  the  country.  The  tectonic  forces 
in  California  are  dominated  by  the  San  Andreas  fault  system,  a strike- 
slip  fault  more  than  600  miles  long,  and  no  comparable  fault  system 
exists  elsewhere  in  the  United  States. 

Nuttli  has  compiled  a list  of  earthquakes  in  the  central  United 
States  since  1843.  He  concludes  that  the  attenuation  rate  of  motion 
with  distance  is  less  in  the  central  United  States,  and  consequently  a 
given  magnitude  earthquake  is  felt  over  a larger  area.  Nuttli  (1973b) 
defines  geographic  regions  (Figure  6-9)  and  associates  maximum  credible 
earthquakes  of  7.2  with  region  1,  6.2  with  region  2,  and  5.7  with  region 
3.  He  further  adds  however,  that  there  is  no  place  in  the  central 
United  States  that  can  be  considered  completely  aseismic. 

For  region  1 the  maximum  credible  earthquake  is  equal  to  that  of 
the  three  largest  earthquakes  of  the  1811-1812  New  Madrid  sequence. 
Since  the  attenuation  with  distance  is  less  in  the  central  United 
States,  the  earthquakes  can  be  expected  to  be  more  destructive.  Earth- 
quake surface  waves  are  dispersed  which  result  in  an  increase  in  the 
time  duration  of  the  wave  motion  with  an  increase  in  epicentral  dis- 
tance. This  is  not  a problem  in  California  because  strong  attenuation 
of  wave  energy  with  distance  reduces  the  motion  to  low  levels  at  dis- 
tances at  which  dispersion  becomes  important.  The  weak  attenuation  in 
the  central  United  States  produces  prolonged  shaking  of  as  much  as  1 or 
2 minutes  at  distances  of  a few  hundred  miles.  Further,  there  is  a 
danger  of  developing  resonant  conditions  for  structures  with  natural 
periods  of  oscillation  near  those  of  the  predominant  surface  motion. 
Nuttli  has  estimated  the  acceleration  with  distance  for  region  1 (7.2 
magnitude)  as  shown  in  Figure  6-10.  The  acceleration  given  in  Figure 
6-10  for  0.3,  1,  and  3 cps  waves  are  the  resultants  of  vertical  and 


horizontal  components  of  the  sustained  maximum  surface  wave  motion. 


(The  resultant  is  the  vector  sum  of  two  horizontal  and  one  vertical 
components.)  This  work  notes  the  attenuation  as  a function  of  wave 
frequency.  If  the  horizontal  motion  is  considered  as  two-thirds  of  the 
resultant,  Figure  6-2  can  be  plotted  as  shown  in  Figure  6-10  to  compare 
Western  attenuation  with  Central  attenuation.  As  can  be  seen  for  dis- 
tances of  less  than  60  miles  (100  km)  there  is  fair  agreement.  Within 
these  distances  attenuation  is  controlled  by  geometric  spreading  rather 
than  absorption  so  that  the  attenuation  does  not  vary  much  with  geo- 
graphic area.  However,  at  greater  distances,  lower  level,  long-duration 
shaking  may  be  a significant  problem,  causing  liquefaction. 

Table  6-7  from  Nuttli  (1973a)  gives  displacement  velocity  and 
acceleration  resultant  data  for  hard  rock  ground  motion.  It  is  impor- 
tant to  note  that  the  active  fault  zones  in  the  central  United  States 
are  not  well-delineated,  in  contrast  to  those  in  California.  There  is 
no  known  evidence  of  surficial  fault  breakage  for  any  central  United 
States  earthquake  except  possibly  those  of  1811  and  1812. 
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Earthquake  Characteristics 


In  developing  a ground  motion  prediction,  it  is  necessary  to  define 
the  distance  from  the  point  of  interest  to  the  fault  under  considera- 
tion. If  the  fault  break  is  short  in  length  and  the  site  is  located  a 
considerable  distance  from  the  fault,  the  significant  distance  from  the 
site  to  the  zone  of  energy  release  can  be  expressed  by  the  epicentral 
distance.  In  the  case  of  a long  fault  this  can  be  grossly  misleading. 
The  rupture  of  the  fault  propagates  along  its  length.  When  the  site  is 
close  to  a fault  the  distance  should  consider  the  release  of  energy  as 
the  rupture  propagates  along  the  fault  length.  In  such  a case,  the 
distance  to  the  site  from  the  zone  of  energy  release  is  more  appropri- 
ately characterized  by  the  shortest  distance  to  the  causative  fault 
rather  than  the  distance  to  the  epicenter.  To  give  perspective  to  this, 
consider  Table  6-8  which  gives  an  indication  of  magnitude  and  length  of 
slipped  fault. 
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After  the  design  earthquake  magnitudes  and  site  ground  motion 
levels  have  been  defined  by  the  geological  and  seismological  studies, 
characteristics  of  the  earthquake  motions  must  be  defined  for  use  in 
engineering  studies.  Time-motion  records  may  be  used  with  complex 
computer  programs  as  input  to  analyses  which  predict  the  response  of 
overlying  soils. 

The  important  characteristics  of  the  acceleration  record  are  its 
duration,  predominant  period,  and  peak  acceleration.  A present  practice 
is  to  scale  acceleration  from  an  existing  suitable  record  to  provide 
earthquake  time  histories.  The  general  problem  of  scaling  earthquake 
records  is  treated  by  Seed,  Idriss  and  Kiefer  (1969).  The  acceleration 
is  linearly  scaled  to  provide  the  design  peak  accelerations. 

The  duration  of  strong  motion  has  been  related  in  a general  way  to 
fault  length  and  the  time  required  for  the  fault  to  shear,  which  can 
also  be  related  to  earthquake  magnitude.  Figure  6-11  provides  a rela- 
tionship between  duration  of  strong  motion  and  earthquake  magnitude 
developed  by  Lee  and  Chan  (1972).  Note  the  definition  of  duration  on 
the  graph.  Data  from  Page  et  al.  (1972)  for  near  fault  horizontal 
motion  is  given  in  Table  6-4.  Studies  by  Bolt  (1974)  and  Kobayashi 
(1974)  have  developed  similar  duration  versus  magnitude  correlations 
(Figures  6-12  and  6-13).  In  general  the  studies  presented  in  Figures 
6-11,  6-12,  and  6-13  show  good  agreement.  Duration  defined  as  motion 
greater  than  O.OSg  can  be  estimated  by  an  average  of  the  data  shown. 
The  rate  of  rupture  is  believed  to  be  on  the  order  of  magnitude  of  2 
miles  per  second  so  that,  for  example,  a magnitude  7 earthquake  result- 
ing from  a 25-mile  (40  km)  fault  break  would  have  at  least  12.5  seconds 
duration.  Nuttli  (1973c)  suggests  that  Central  United  States  earth- 
quakes do  not  have  the  same  duration  (of  felt  motion)  characteristics  as 
Western  earthquakes  as  a result  of  the  increased  dispersion  effects  and 
weak  attenuation.  However,  for  distances  less  than  100  miles  from  the 
fault  this  may  not  be  significant. 

Figure  6-14  provides  a relationship  between  predominant  periods  of 
maximum  acceleration,  earthquake  magnitude,  and  distance  from  the  causa- 
tive fault.  The  scarcity  of  close-in  data  results  in  a constant  period 


72 


for  a given  earthquake  Magnitude  out  to  a distance  of  25  Biles  (40  km) 
froa  the  causative  fault.  Figueroa  (1960),  however,  indicates  a very 
wide  scatter  in  predoainant  period  data  out  to  a distance  of  50  ailes 
(80  ka). 

An  acceleration  tiae  history  for  use  in  seisaic  design  studies  can 
also  be  generated  with  desired  characteristics  by  use  of  randoa  number 
prograas  shaped  to  the  desired  spectral  characteristics. 

Scaling  an  individual  earthquake  record  has  the  disadvantage  that 
each  record  is  representative  of  a specific  event  and  site.  Each  indi- 
vidual record  will  be  deficient  in  soae  response  frequencies.  An  ensem- 
ble of  scaled  earthquake  records  can  be  used  as  a better  average  of  the 
individual  records.  Scaled  artificial  earthquake  records  that  do  not 
exhibit  a specific  bias  can  also  be  used.  Procedures  for  selecting  an 
ensemble  of  real  or  artificial  records  are  provided  in  Werner  (1970)  and 
Guzaan  and  Jennings  (1975a  and  b). 

Seisaic  Studies 

As  an  example  of  available  seisaic  studies,  Greensfelder  (1974) 
shows  contours  of  possible  peak  acceleration  froa  active  faults  using 
the  attenuation  after  Schnabel  and  Seed  (1972).  The  probability  of 
occurrence  is  not  considered  other  than  the  fact  that  the  faults  consid- 
ered active  include  Quaternary  movement.  Algeraissen  and  Perkins  (1976) 
developed  a similar  aap  for  the  entire  United  States.  This  may  assist 
in  predicting  site  motion.  However,  the  relatively  large  area  coverage 
of  the  maps  may  not  be  sufficiently  detailed  for  specific  site  studies. 

Discussion 

An  in-depth  geological  and  seisaological  investigation  is  desirable 
at  proposed  sites  to  locate  faults  and  evaluate  the  site  soil  profile. 
When  historical  data  is  available  the  design  earthquake  magnitude  should 
be  determined  in  terms  of  the  structure  life  and  earthquake  recurrence. 
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For  engineering  analysis,  use  of  the  Poisson  model  should  be  adequate. 

Figure  6-10  shows  that  although  the  attenuation  in  Central  and  Eastern 
United  States  is  less  than  in  the  West  within  the  range  less  than  60 
■iles  (100  km),  Figure  6-2  is  probably  adequate  since  this  is  the  range 
of  most  engineering  interest.  The  earthquake  duration  can  be  approxi- 
mated from  Figures  6-11,  6-12,  and  6-13. 

1 
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(km)  -log10AQ(R)  R (km)  -log^A^R)  | R (km)  -log^A^R)  I R (km)  -log 


Table  6-2.  Means  and  Standard  Deviations  of  the  Logaritha  of  the  Magnitude-Dependent 
Scaling  Function  A (M) 


Table  6-3.  Mean  Value  (Intermediate  Site)  (from  Ferritto  and  Forrest,  1977) 
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Table  6-3.  Continued 


Table  6-3.  Continued 
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Table  6-7.  Continued 


The  hard-ruck  ground  motion*  (displacement,  particle  velocity,  and  acceleration)  are  the  vector  resultants  of  the  vertical  and  horizontal  components  of  the 
sustained  maximum  surface  wave  motion.  At  distances  of  75  miles,  the  duration  of  this  motion  will  be  as  much  as  30  s.  At  distances  of  100  miles  and 
peater,  the  duration  can  be  as  much  as  1 to  2 minutes,  generally  increasing  with  increasing  epicentral  distance. 


(a)  Generation  ami  tranamtaatou  of  aeiaaiti'  wave*. 


Figure  6*1.  Factor*  affecting  ground  shaking  at  a aitr  (froai 
Shannon-Wi laon  and  Aghahian  Associates,  1975). 


(h)  l.iat  of  factor*  affecting  aeiaaiic  wave  trananiaaion. 
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Figure1  6-5 . Relationship  between  maximum  ground  acceleration  and 
distance  from  xone  of  energy  release  for  deep 
cohesionlesa  soil  conditions  (from  Seed  et  al.,  1975), 
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Figure  6-6.  Approximate  relationships  between  maximum 
accelerations  on  rock  and  other  local  site 
conditions  (from  Seed  et  al.,  1975). 


I arth«|uakc  Magnitude 

5 0 5 9 4 1970  l ytic  ( reck  m * 5 4 
1966  Park  field  m • 5.5 
1957  Daly  City  m • 5.3 

6.0- 6  9 % 1971  San  I crnando  m 

196K  Borrego  Mountain  m 

7.0-  7 9 4 1952  Kern  County  m * 7.7 

1949  Puget  Sound  m « 7. 1 


Distance  (km) 


Peak  horizontal  acceleration  versus  distance  to 
slipped  fault  as  a function  of  nagnitude. 
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Attenuation  equations  for  magnitude  0.5  compared 
to  data  from  strong  motion  stations,  recording 
the  February  9,  1971  San  Fernando  earthquake 
(from  Ferritto  and  Forrest,  1977). 
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Figure  6-10.  Ground  accelerations  on  hard  rock  for  the  maximum 

credible  earthquake  of  seisaiic  region  1 (New  Madrid 
faulted  zone)  (froai  Ferritto  and  Forrest,  1977). 
Accelerations  for  wave  frequencies  of  0.3,  1,  and 
3 Hz  are  the  resultants  of  the  vertical  and  horizontal 
components  of  the  sustained  maximum  surface-wave  motion 
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Figure  6-13.  Kobayashi  (1974)  results  for  duration 


Distance  from ( aiisattvv  I atilt  inn) 
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CHAPTER  7 - DEVELOPMENT  OF  RECURRENCE  ANALYSIS  PROCEDURE 


Requirements  for  Design 


The  current  NAVFAC  design  procedure  prescribes  the  selection  of  a 
design  level  earthquake  with  10%  chance  of  being  exceeded  in  25  years. 
As  such,  a probabilistic  approach  is  specified  which  must  relate  site 
activity  to  structural  design  levels.  Chapter  U demonstrated  the  limi- 
tations of  historic  data  when  used  as  the  sole  data  source  for  recur- 
rence trends.  Chapter  3 pointed  out  the  use  of  slip  rates  in  the  compu- 
tation of  recurrence  trends  and  the  associated  limitations.  The  best 
engineering  estimate  of  earthquake  recurrence  can  be  made  by  use  of 
historic  data  supplemented  by  fault  slip  data. 


Computer  Analysis 


An  automated  procedure  has  been  developed  to  perform  the  seismic 
analysis  using  the  NOAA  data  base  as  input  (Appendix  A).  The  specific 
region  of  the  study  is  specified  in  terms  of  bounding  latitude  and 
longitude.  The  general  data  file  of  earthquakes  is  sorted,  and  the 
subset  of  earthquakes  in  the  region  becomes  the  data  base  for  the  spe- 
cific study.  The  location  of  the  site  is  specified.  The  methodology 
proposed  has  been  automated  as  a computer  program.  The  program 
(Appendix  B)  determines  the  distance  between  the  site  and  each  epicenter 
in  the  study  data  base  and,  using  the  earthquake  average  magnitude  and 
the  Trifunac  and  Brady  attenuation  relationship,  computes  the  site 
acceleration  for  all  historic  earthquakes;  a histogram  is  created.  The 
earthquake  magnitudes  are  then  plotted  and  a least-squares  analysis  used 
to  obtain  recurrence  coefficients.  A map  of  epicenters  is  plotted. 
These  help  in  defining  the  regional  seismicity. 
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Further  study  ran  be  made  by  specifying  the  location  of  a fault  in 
terns  of  coordinatec  of  several  points  defining  line  segments.  All 
earthquakes  within  a specified  distance  from  the  fault  line  under  study 
are  made  a subset  and  the  recurrence  of  the  fault  calculated  as  above. 
A probability  analysis  is  then  perfomed  to  calculate  expected  site 
acceleration  and  causative  earthquake  sugnitude  and  epicentral  distance. 
The  progran  randomly  selects  the  epicenter  of  an  earthquake.  Using  the 
fault  recurrence  data  in  terms  of  Richter  coefficients  and  maximum 
earthquake  magnitude  associated  with  the  fault,  the  program  determines 
an  earthquake  magnitude  and  the  length  of  fault  break  which  would  occur. 
Using  the  random  epicenter  location  and  the  fault  break  (assumed 
centered  on  the  epicenter),  the  closest  distance  of  the  site  to  fault 
break  is  calculated.  This  distance,  along  with  the  acceleration- 
magnitude  attenuation  relationship,  gives  the  site  acceleration.  The 
process  is  repeated  to  produce  a list  of  site  accelerations  and  related 
causative  magnitude  and  epicentral  distance.  A Poisson  distribution  is 
used  to  compute  a probability  distribution.  Random  location,  magnitude 
and  acceleration  level  are  all  considered  in  the  determination  of  site 
acceleration.  Provisions  are  included  to  use  recurrence  data  from  slip 
analysis  or  regional  seismicity  in  lieu  of  fault  specific  data.  Further 
provision  is  included  for  floating  earthquakes. 

Experts  disagree  on  whether  main  shocks  alone  or  main  shocks  and 
aftershocks  may  be  used  to  estimate  recurrence;  decision  on  which  to  use 
is  independent  of  the  computer  analysis.  It  is  felt  that  the  use  of 
aftershocks  is  a valid  procedure  since  these  events  do  represent  the 
release  of  strain  energy.  Further,  the  inclusion  of  aftershocks  will 
yield  conservative  results. 

As  stated  above  a single  recent  large  event  may  release  strain 
built  up  over  hundreds  or  thousands  of  years.  As  such  it  might  indicate 
a period  of  less  activity  in  the  immediate  future.  However,  since  the 
data  base  is  relatively  short,  the  return  for  this  event  might  be  erro- 
neously indicated  as  much  greater.  For  example,  if  this  were  a 500-year 
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event  and  occurred  during  a 50-year  data  base  its  return  night  be  esti- 
mated at  0.02  rather  than  0.002.  The  plotted  data  points  and  line  of 
best  fit  are  determined  by  the  computer  analysis  using  regression  analy- 
sis techniques.  These  should  be  reviewed  and  judgment  used  to  adjust 
this  type  of  data  point  which  will  clearly  plot  significantly  higher 
than  the  linear  portion  of  the  recurrence  data. 

Comparison  With  Fault  Slip  Data 

Chapter  3 presents  procedures  for  estimating  recurrence  data  based 
on  fault  slip.  Table  3-1  gives  estimates  of  recurrence  data  for  several 
faults  based  on  slip  rate  and  fault  length.  Figures  7-1,  7-2,  and  7-3 
show  maps  locating  the  faults  and  epicenters  associated  with  the 
regions.  Figure  7-4  gives  plots  of  the  recurrence  data  showing  the  data 
points  based  on  historic  earthquakes  and  estimates  based  on  slip  rate. 
For  the  most  part  agreement  seems  reasonable.  Although  this  may  not  be 
the  case  for  all  faults,  agreement  between  slip  data  and  historic  data 
gives  confirmation  and  added  assurance  of  the  recurrence  coefficients. 

The  San  Andreas  fault  slip  data  were  adjusted  for  the  length  of 
fault  actually  covered  in  the  region  where  historic  data  were  taken.  It 
should  be  noted  that  the  activity  over  its  great  length  is  quite  vari- 
able, covering  the  whole  range  from  locked  to  very  active. 

The  White  Wolf  fault  is  predominantly  a reverse  fault.  The  fault 
was  most  active  during  the  Pleistocene  and  Holocene  period.  The  slip 
rate  was  estimated  based  on  10,000  feet  (3,000  meters)  of  vertical 
displacement  of  Miocene  rock  over  a period  of  7 million  years.  This 
averaging  is  probably  insufficient  to  account  for  more  geologically 
recent  activity  exhibited.  The  historic  data  is  biased  by  the  1952  Kern 
County  earthquake  magnitude  7.7  and  its  aftershocks  which  resulted  in 
vertical  displacement  of  up  to  1.2  m (Lamar  et  al.,  1973).  The  Garlock 
fault  has  had  relatively  little  historic  activity.  The  slip  rate  indi- 
cates recurrence  intervals  greater  than  historic  data. 
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Combination  of  Fault  Slip  Data  With  Historical  Data 

The  coabination  of  the  recurrence  data  based  on  fault  slip  rates 
and  the  recurrence  data  based  on  the  historical  data  is  an  area  involv- 
ing a great  deal  of  judgment.  Judgment  may  be  used  directly  to  adjust 
the  recurrence  coefficients. 

An  alternative  technique  developed  by  Campbell  (1977)  utilises  a 
Bayesian  technique  for  estimation  of  seismic  risk.  This  technique 
combines  the  two  data  bases,  which  are  founded  on  the  uncertainty  asso- 
ciated with  each.  Campbell  (1977)  has  derived  a procedure  for  estimat- 
ing slip  based  on  seismic  moment,  shear  modulus,  fault  area,  and  the 
historical  magnitude- frequency  coefficients.  Although  based  on  empiri- 
cal data  in  part,  this  model  has  not  as  yet  been  extensively  validated 
against  known  fault  behavior.  The  calculated  mean  rate  of  occurrence 
and  coefficient  of  variation  based  on  fault  slip  is  used  as  the  prior 
estimate.  This  is  "updated"  by  the  historical  data.  Although  the 
approach  has  been  automated,  considerable  judgment  is  still  required. 
Great  care  is  required  to  avoid  arbitrarily  biasing  the  results  in  favor 
of  either  the  slip  rate  estimates  or  historic  data.  The  size  of  the 
historic  data  base  can  totally  dominate  the  procedure.  Only  main  shocks 
are  used  in  this  procedure.  This  will  limit  the  number  of  events  con- 
siderably for  less  active  faults. 

The  procedure  developed  by  Campbell  (1977)  does  provide  a rigorous 
means  of  combining  seismotectonic  data  with  local  historic  seismicity. 
The  procedure  may  be  an  improvement  if  applied  with  good  judgment  where 
data  are  available. 

Case  Study 

The  Naval  Construction  Battalion  Center  (NCBC),  Port  Hueneme  is 
located  in  southern  California  at  longitude  11. 2W,  latitude  34.15N.  The 
time  period  for  the  study  was  25  years  to  comply  with  the  NAVFAC  design 
instruction.  During  the  60  years  of  historic  data,  3,394  earthquakes 


occurred  in  the  region  around  NCBC;  most  were  very  minor.  The  largest 
event  was  the  1952  7.7  magnitude  earthquake  in  the  White  Wolf  fault. 
This  event  produced  a site  acceleration  of  about  0.04g.  The  highest 
site  acceleration  was  from  the  1973  magnitude  5.6  Point  Mugu  earthquake 
and  is  estimated  to  have  had  a ground  acceleration  of  0.15g- 

As  an  example  of  the  methodology  developed  the  following  is  pre- 
sented. Figure  7-3  shows  the  California  coast  (longitude  118  to  120N, 
latitude  33.5  to  35N),  NCBC,  and  all  epicenters  and  faults.  Figure  7-5 
gives  an  outline  of  the  computational  procedures.  Figures  7-4  and  7-6 
give  the  regional  recurrence  data;  Figure  7-7  the  regional  probability 
data.  The  regional  recurrence  data  although  interesting  are  limited  in 
its  applicability  to  engineering  analysis.  Further,  the  data  are  very 
dependent  on  the  area  (or  radius  from  the  site)  of  coverage.  Figure  7-8 
shows  the  epicenters  associated  with  the  Big  Pine  fault.  Figure  7-9 
gives  the  recurrence  data  for  the  Big  Pine  fault.  Also  shown  is  the 
recurrence  data  from  fault  slip  studies  (Chapter  3).  As  can  be  seen, 
agreement  is  very  good  in  this  case.  Figure  7-10  gives  probability  of 
occurrence  as  a function  of  magnitude;  Figure  7-11  gives  probability  of 
acceleration  from  an  event  on  the  Big  Pine  fault  at  the  NCBC  site  not 
being  exceeded.  Associated  with  the  probability  of  acceleration  is  the 
related  causative  earthquake  magnitude  and  epicentral  distance.  In  a 
similar  manner  Table  7-1  was  prepared  for  all  the  faults  affecting  the 
site. 

With  the  probability-acceleration  data  for  each  fault  (such  as 
Figure  7-11)  it  is  possible  to  compute  the  total  risk  P of  not  exceeding 
an  acceleration  level  from  all  the  faults  by 

n 

P = n (P.) 
i*l 

where  i - 1 to  n faults 

P.  = probability  of  not  exceeding  acceleration 

In  this  example,  Table  7-2  summarizes  the  data. 
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This  exaaple  shows  a procedure  for  determining  site  acceleration. 
Response  spectra  can  be  produced  using  the  site  acceleration  for  each 
fault,  which  should  be  exaained  to  determine  local  characteristics  for 
structural  design.  This  will  be  discussed  in  the  next  section. 

It  is  iaportant  to  note  that  the  procedure  expresses  a cumulative 
total  aeisaic  risk  and  results  in  values  of  acceleration  at  a confidence 
level  in  exceas  of  those  for  any  individual  fault.  The  aore  traditional 
design  procedure  in  use  would  siaply  utilise  Table  7-1  directly  in 
constructing  a design  spectrua  assuaing  the  acceleration  produced  on  one 
or  several  faults  would  govern  structural  design.  The  selection  of 
individual  fault  risk  or  total  risk  is  independent  of  the  procedures 
specified  herein.  It  is  not  the  intent  of  this  report  to  prescribe  a 
design  criteria. 


Table  7-1.  Site  Acceleration 


Maximum  Peak  Site 
Acceleration  in  25  yr 


Fault 


Confidence  Level 


Big  Pine 
Malibu  Canyon 
Oakridge 
Santa  Ynez 
Arroyo  Parida 
San  Andreas 
San  Gabriel 
Santa  Cruz 
Santa  Susana 


Limited  by  maximum  credible  event  as  shown  in  parentheses 
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Figure  7-1.  Map  of  the  San  Andreas,  Elsinore,  San  Jacinto,  and  Sierra  Madre  fau] 
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(b)  San  Andreas  South. 


Figure  7-4.  Continued 


(d)  Sierra  Madre. 


Figure  7-4.  Continued 
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Figure  7-7.  Probability  of  occurrence,  Region  118-120  W 
33.5-35.0  N,  25  yeart. 


Figure  7-8.  Portion  of  sup  showing  epicenters  associated 
with  Big  Pine  fault. 
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Figure  7-10.  Probability  of  occurrence,  Big  Pine  fault,  25  years 
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Figure  7-11.  Probability  of  aite  acceleration  not  being 
exceeded  in  25  years,  Big  Pine  fault. 


CHAPTER  8 - RESPONSE  SPECTRA  AND  ANALYTICAL  TECHNIQUES 


Engineers  have  found  it  useful  to  examine  the  frequency  content  of 
earthquake- induced  ground  motion.  Strong  motion  accelerograms  have  been 
analyzed  as  a means  of  obtaining  further  insight  into  ground  motion. 

The  first  technique  available  is  Fourier  analysis.  The  Fourier 
spectrum  of  an  acceleration  shows  the  significant  frequency  characteris- 
tics of  the  recorded  motion.  The  Fourier  spectrum  is  defined  as 


over  the  interval  0<x<T.  The  acceleration  is  zero  outside  the  limits  0 

to  T.  The  Fourier  amplitude  spectrum  is  given  by  the  square  root  of  the 

sum  of  the  squares  of  the  real  and  imaginary  parts  of  F(u>).  Associated 

with  the  magnitude  F(u>)  is  the  phase  angle  $(u»),  defined  as  the  arc 

tangent  of  the  imaginary  part  divided  by  the  real  part.  Units  of  the 

spectrum  reflect  the  time  integration;  thus,  for  an  acceleration  in 
2 

ft/sec  , the  spectrum  magnitude  will  be  in  feet  per  second  and  the  phase 
angle  will  be  in  radians.  The  Fourier  spectrum  magnitude  and  phase 
angle  represent  only  the  input  motion.  Through  convolution  this  may  be 
combined  with  transfer  functions  for  other  elements  (such  as  soil  struc- 
ture interaction)  assuming  elastic  behavior.  The  Fourier  spectrum 
magnitude  and  phase  are  a complete  record  which  is  unique  and  maintains 
the  total  time  history  record.  The  time  history  record  may  be  recreated 
by  reverse  transformation.  Generally,  only  the  magnitude  of  the  Fourier 
spectrum  is  shorn  in  reference  illustrations. 
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In  earthquake  engineering  it  is  important  to  be  able  to  determine 
the  magnitude  of  maximum  response  of  a structure.  This  has  given  rise 
to  the  response  spectra  technique.  A single-degree-of-freedom,  spring- 
mass-damper  system  can  be  analyzed  and  its  time  history  of  displacement 
calculated  to  determine  relative  displacement  between  the  mass  (the 
structure)  and  the  excited  base  (the  ground).  Relative  velocity  and 
relative  acceleration  may  also  be  calculated.  However,  of  primary 
interest  for  engineering  applications  is  the  maximum  absolute  values  of 
structure  relative-displacement,  structure  relative-velocity  and  abso- 
lute structure  acceleration.  These  values  SD,  SV,  SA  are  functions  of 
the  critical  damping.  Plots  of  SD,  SV,  SA  versus  the  undamped  natural 
period  of  vibration  and  for  various  fractions  of  critical  damping  are 
called  response  spectra.  In  typical  engineering  structures  the  damping 
is  small  and  for  harmonic  excitation  the  following  holds 

SD  = (-L)SV 

SA  = (j^SV 

For  earthquake-like  excitations  which  are  not  strictly  harmonic 
excitation  an  engineering  assumption  is  made  that  the  above  is  still 
accurate.  The  following  definitions  are  used 

PSV  = (f^SD 
2 

PSA  = (j1)  SD 

where  PSV  is  the  pseudo-relative  velocity  and  PSA  is  the  pseudo-absolute 
acceleration.  The  term  "pseudo"  is  used  to  recognize  the  assumptions 
made  concerning  small  damping  and  harmonic  motion.  Thus  SD,  PSV,  PSA 
and  T make  up  a set  of  data;  knowing  any  two  makes  it  possible  to  deter- 
mine the  other  two.  This  unique  relationship  makes  it  possible  to  plot 
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response  spectra  in  tripartile  form.  The  response  spectra  shows  the 
response  of  a single-degree-of-freedom  system  (structure)  as  a function 
of  damping  of  the  system  and  period  for  the  given  input  acceleration. 

A comparison  of  Fourier  and  pseudo-velocity  response  spectrum 
reveals  that  both  are  of  the  same  units.  For  an  undamped  oscillator  a 
similar  mathematical  relationship  exists  between  the  Fourier  amplitude 
spectrum  and  exact  relative  velocity  response  spectrum.  (The  Fourier 
spectrum  may  be  viewed  as  the  maximum  velocity  of  the  undamped  oscilla- 
tor in  the  free  vibration  following  the  earthquake;  however,  the  exact 
velocity  response  spectrum  is  the  maximum  velocity  during  both  the 
earthquake  and  the  subsequent  free  vibration.) 

The  Fourier  amplitude  spectrum  is  the  quantity  most  frequently  used 
by  seismologists  in  their  investigations  of  the  earthquake  ground 
motion.  The  response  spectrum  is  generally  used  by  structural  engineers 
in  the  design  of  structures.  Generally  the  pseudo-velocity  response 
spectrum  is  the  upper  envelope  of  the  Fourier  spectrum. 

In  summary  the  most  representative  measure  of  the  driving  ground 
motion  is  the  Fourier  spectrum,  which  may  be  plotted  in  tripartile  form. 
The  response  spectrum  gives  the  response  of  a single-degree-of-freedom 
structure  to  the  ground  motion.  Amplification  of  motion  occurs  where 
resonant  components  of  motion  interact  with  the  structure.  The  response 
spectrum  is  most  widely  used  in  structural  engineering  with  the  modal 
analysis  technique.  This  approach  determines  the  eigenvalues  and  mode 
shapes  for  a number  of  the  modes  of  structure,  and  using  the  modal 
period  determines  the  acceleration  to  be  used  in  computing  an  equivalent 
maximum  static  force  to  be  applied  to  the  structure.  A response  spec- 
trum cannot  be  used  directly  as  input  to  a dynamic  structural  analysis 
to  generate  a time  history  response  because  the  response  spectrum  does 
not  contain  a measure  of  the  duration  of  excitation.  The  Fourier  spec- 
trum does  contain  the  full  earthquake  representation  and  may  be  used  in 
conjunction  with  transfer  functions  to  compute  the  structure  response. 
This  technique  is  not  used  in  general  structural  engineering  because  of 
the  difficulty  in  determining  the  transfer  function.  The  more  common 
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structural  analysis  techniques  which  produce  dynamic  time  histories  of 
response  require  input  acceleration  histories.  Random  analysis  programs 
exist  to  generate  time  histories  from  a given  spectral  envelope;  how- 
ever, the  duration  must  be  specified.  The  resulting  randomly  generated 
signal  is  not  unique,  and  any  number  may  be  generated  from  a single 
spectral  envelope. 

Time  history  records  of  actual  earthquakes  are  available  from  the 
California  Institute  of  Technology  National  Earthquake  Information 
Center.  These  may  be  used  when  required  for  time  history  input  of 
ground  motion.  Care  is  needed  in  selection  since  most  records  are 
recorded  in  structures  which  may  influence  the  recording.  Also  the 
region  around  the  site  may  influence  the  site  response.  The  propagation 
of  seismic  waves  is  influenced  by  the  local  geology  and  soil  conditions. 
The  greater  the  extent  of  softer  soils  over  bedrock,  the  less  the  bound- 
ary effects  of  the  bedrock  will  have  on  the  site.  The  depth  of  soil 
overlaying  bedrock  affects  the  period  of  vibration  of  the  ground.  This 
establishes  a fundamental  soil  frequency  of  particular  importance  on 
soil-structure  interaction.  Further,  this  is  a factor  in  determining 
the  frequencies  of  waves  filtered  out  by  the  soil,  thus  directly  affect- 
ing the  time  history  record. 


Newmark-Hall  Procedures  for  Elastic  Spectra 

Newmark  (1967)  based  on  studies  of  response  spectra  noted  that 
response  spectra  could  be  related  to  peak  ground  acceleration,  velocity, 
and  displacement.  From  this  study  it  was  possible  to  develop  standard 
shapes  for  use  in  structural  analysis.  The  standard  ground  motion 
spectra  (lg,  48  in. /sec,  36  in.)  could  be  scaled,  based  on  peak  hori- 
zontal ground  acceleration  expected  at  the  site.  Figure  8-1  gives 
amplification  factors  which  could  then  be  applied  to  estimate  structural 
response . 

McGuire  (1977)  also  demonstrates  ratios  of  damped  spr  ra . He 
found  that  the  5%  damped  pseudo-velocity  spectra  have  values  approxi- 
mately 70%  to  80%  of  the  2%  damped  spectra  and  that  the  10%  damped 
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spectra  have  values  of  approximately  55%  to  65%  of  the  2%  damped  spec- 
tra. These  ratios  are  constant  throughout  frequency  range  independent 
of  the  type  of  earthquake  distance  from  site  and  magnitude.  This  con- 
firms the  Newmark-Hall  approach. 

At  a frequency  of  about  6 cps,  the  amplified  acceleration  region 
line  intersects  a line  sloping  down  toward  the  maximum  ground  accelera- 
tion value  and  intersecting  that  line  at  various  frequencies,  depending 
on  the  damping.  The  intersection  is  at  a frequency  of  about  30  cps  for 
2%  damping.  These  lines  are  designated  as  the  acceleration  transition 
region  of  the  spectra.  Finally,  beyond  the  intersection  with  the  maxi- 
mum ground  acceleration  line,  the  response  spectrum  continues  with  the 
maximum  ground  acceleration  value  for  higher  frequencies. 

The  spectra  so  determined  can  be  used  as  design  spectra  for  elastic 
responses.  The  spectra  are  completely  described  when  the  maximum  ground 
motion  values  are  given  for  the  three  components  of  ground  motion,  and 
the  damping  is  known.  When  only  the  maximum  ground  acceleration  is 
given,  the  values  used  for  maximum  ground  velocity  and  displacement  are 
taken  as  proportional  to  those  in  the  figure,  or  as  scaled  by  the  same 
scale  factor  relative  to  the  maximum  ground  acceleration  compared  with 
lg- 

An  assumption  is  made  that  acceleration,  velocity,  and  displacement 
values  are  proportional  to  one  another,  independent  of  magnitude  of 
motion.  The  shape  is  thought  correct  for  sites  on  firm  ground,  soft 
rock,  or  competent  sediments.  However  for  soft  sediments,  the  veloci- 
ties and  displacements  must  be  increased.  Garcia  and  Roesset  (1970) 
performed  studies  comparing  actual  spectra  with  those  estimated  by  the 
Newmark-Hall  procedure.  Results  show  favorable  comparison  indicating 
the  utility  of  the  Newmark-Hall  procedure. 

Vertical  spectra  may  be  estimated  by  taking  two-thirds  of  the 
horizontal  spectra  when  fault  movements  are  horizontal  and  by  taking 
horizontal  spectra  when  large  vertical  motions  are  expected. 
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Newmark-Hal 1 Procedures  for  Inelastic  Response  Spectra 


The  elastic  spectra  discussed  in  the  previous  section  may  be  ad- 
justed to  approximate  inelastic  behavior  of  a structure.  The  displace- 
ment region  and  the  velocity  region  are  divided  by  the  structure  ductil- 
ity p to  obtain  yield  displacement  D'  and  velocity  V'  (Figure  8-2).  The 
acceleration  region  (right  side)  is  relocated  by  choosing  it  at  a level 
which  corresponds  to  the  same  energy  absorption  for  elasto-plastic 
behavior  as  for  elastic  for  the  same  period  of  vibration. 

The  extreme  right-hand  portion  of  the  spectrum,  where  the  response 
is  governed  by  the  maximum  ground  acceleration,  remains  at  the  same 
acceleration  level  as  for  the  elastic  case,  and  therefore  at  a cor- 
responding increased  total  displacement  level.  The  frequencies  at  the 
corners  are  kept  at  the  same  values  as  in  the  elastic  spectrum.  The 
acceleration  transition  region  of  the  response  spectrum  is  now  drawn 
also  as  a straight  line  transition  from  the  newly  located  amplified 
acceleration  line  and  the  ground  acceleration  line,  using  the  same 
frequency  points  of  intersection  as  in  the  elastic  response  spectrum. 

In  all  cases  the  "inelastic  maximum  acceleration"  spectrum  and  the 
"inelastic  maximum  displacement"  spectrum  differ  by  the  factor  p at  the 
same  frequencies.  The  design  spectrum  so  obtained  is  shown  in  Figure 
8-2.  Both  the  maximum  displacement  and  maximum  acceleration  bounds  are 
shown  for  comparison  with  the  elastic  response  spectrum. 

The  solid  line  DVAAq  shows  the  elastic  response  spectrum.  The 
heavy  circles  at  the  intersections  of  the  various  branches  show  the 
frequencies  which  remain  constant  in  the  construction  of  the  inelastic 
design  spectrum. 

The  dashed  line  D'V'A'A  shows  the  inelastic  acceleration,  and  the 

o 

line  DVA"A^  shows  the  inelastic  displacement.  These  two  differ  by  a 
constant  factor  p for  the  construction  shown,  but  A and  A’  differ  by  the 
factor  V 2p-l , since  this  is  the  factor  that  corresponds  to  constant 
energy. 
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Applied  Technology  Council  Spectra  Example 
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Using  the  Newmark-Hall  procedure  the  Applied  Technology  Council 
(1974)  developed  basic  response  spectra  for  the  "damage  threshold"  and 
"collapse  threshold"  earthquakes.  Using  a structure  life  of  70  years 
ground  acceleration  was  selected  for  that  study  (representative  for  a 
number  of  sites  in  southern  California)  such  that  probabilities  of 
occurrence  of  50%  and  10%  for  damage  and  collapse  thresholds  respective- 
ly (Figures  8-3  and  8-4).  Ground  acceleration  values  of  0.24g  for  the 
damage  threshold  and  of  0.28g  for  the  collapse  threshold  were  used  for 
their  specific  site. 

Table  8-1  shows  three  categories  of  structure  with  ductility  p'  and 
damping  0'  specified.  Figures  8-5,  8-6  and  8-7  show  inelastic  design 
spectrum  for  the  three  types  of  structures. 

Regulatory  Guide  1.60  Procedures 

Newmark  (1973),  Blume  (1973),  and  Newmark,  Blume  and  Kapur  (1973) 
presented  studies  of  the  response  spectra.  These  form  the  basis  for  the 
Nuclear  Regulatory  Commission  (AEC)  Regulatory  Guide  1.60  (1973).  The 
approach  is  similar  to  the  Newmark-Hall  procedure.  The  response  spectra 
are  summarized  in  Shannon  and  Wilson  (1975)  quoted  as  follows: 


The  response  spectra  bounds  are  established  by  five 
straight  lines,  similar  to  the  Newmark-Hall  method,  but  the 
control  points  are  designated  by  letters  A,  B,  C,  and  D,  and 
have  definite  assigned  frequencies  for  all  horizontal  spectra. 
Amplification  factors  for  constructing  the  response  spectra 
are  again  a function  of  the  percentage  of  critical  damping, 
and  are  specified  in  ***  [Figure  8-8]  for  the  four  control 
frequencies.  Response  spectra  for  O.lg  maximum  horizontal 
ground  acceleration  have  been  constructed  according  to 
Regulatory  Guide  1.60,  and  are  provided  in  ***  [Figure  8-8). 
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Vertical  response  spectra  are  constructed  in  the  same 
manner  as  horizontal  response  spectra;  however,  the  frequency 
for  control  point  C is  at  3.5  Hz  rather  than  at  2.5  Hz,  and 
the  amplification  factors  are  different  ***  (Figure  8-9). 

Also,  the  maximum  horizontal  components  of  displacement  and 
acceleration  are  used  as  scaling  factors  in  constructing  the 
vertical  response  spectra  ***.  The  amplification  factors  in 
Figure  ***  (8-9)  have  been  established  so  that  the  resulting 
vertical  response  spectra  will  be  two-thirds  of  the  values 
computed  for  the  horizontal  spectra,  at  frequencies  below  0.25 
Hz.  In  the  range  of  frequencies  from  0.25  Hz  to  3.5  Hz,  the 
ratio  of  vertical  to  horizontal  spectral  amplitudes  varies 
from  two-thirds  to  one;  whereas,  for  frequencies  between  3.5 
Hz  and  33  Hz,  the  horizontal  and  vertical  response  spectra  are 
identical.  Beyond  33  Hz,  the  vertical  spectral  accelerations 
decrease  from  a value  equal  to  the  peak  horizontal  ground 
acceleration  at  33  Hz  to  two-thirds  of  the  horizontal  acceler- 
ation at  50  Hz. 

This  specification  of  the  vertical  spectral  acceleration  amplitudes 
above  33  Hz  represents  a correction  to  that  currently  indicated  in 
Revision  1 of  Regulatory  Guide  1.60.  As  pointed  out  by  Valera  (1974), 
this  correction  results  in  a vertical  spectrum  shape  that  is  consistent 
with  the  original  Newmark,  Blume  and  Kapur  (1973)  paper.  The  resulting 
family  of  vertical  response  spectra  are  shown  in  Figure  8-9. 

Approximately  one-quarter  of  the  records  considered  in  the  develop- 
ment of  the  Regulatory  Guide  1.60  spectra  were  recorded  on  rock  sites, 
several  were  recorded  on  deep  cohesionless  soil  sites,  and  the  remaining 
records  were  obtained  on  stiff  soil  sites.  It  should  be  noted  that 
Newmark  Services  (1973),  and  Blume  and  Associates  (1973),  in  the  course 
of  their  independent  investigations  that  led  to  Regulatory  Guide  1.60, 
considered  potential  effects  of  local  soil  conditions;  in  fact,  the 
Newmark  study  provided  spectra  that  indicated  important  differences 
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between  motions  recorded  on  sediments  and  on  competent  rock,  particular- 
ly for  frequencies  below  about  3 Hz.  However,  the  recommended  design 
spectra  contained  in  the  subsequent  paper  by  Newmark,  Blume,  and  Kapur 
(1973)  and  in  Regulatory  Guide  1.60  do  not  differentiate  between  various 
local  soil  profiles  and,  therefore,  represent  a statistical  sample  in 
which  the  subsurface  soil  conditions  are  an  uncontrolled  variable. 


The  data  base  of  strong  motion  records  can  be  a useful  source  of 
seismic  data.  Records  may  be  selected  to  represent  seismologic,  geolog- 
ic, and  local  site  conditions.  Selection  is  complicated  by  a number  of 
factors.  Ideally,  the  records  should  be  selected  to  match  source-site 
transmission  path  source  mechanism  and  local  site  conditions.  These  are 
not  readily  quantifiable.  Thus,  reliance  is  made  on  earthquake  magni- 
tude site  acceleration  level,  site  classification,  and  duration  of 
motion.  Judgment  is  an  important  factor  in  selecting  and  scaling 
records . 

Appendix  C presents  a CF.1,  program  which  uses  the  available  data 
optimizing  the  selection  considering  magnitude  distance,  acceleration, 
and  site  category.  This  should  assist  in  reducing  the  problem  of 
selecting  the  closest  records  available.  Appendix  1)  describes  a CEL 
plotting  program.  Appendix  K describes  a CEL  program  to  plot  time 
histories. 

Guzman  and  Jennings  (1975a  and  b)  describe  a procedure  for  develop- 
ing criteria  spectra  based  on  site-matched  records.  Geological  and 
seismologi ca  1 factors  should  be  evaluated.  A group  of  strong  motion 
records  is  selected  representing  comparable  magnitude,  distance  geology, 
and  soil  conditions.  Accelerograms  are  added  to  the  group  until  condi- 
tions under  which  they  were  recorded  are  no  longer  representative  of  the 
earthquake  or  until  no  additional  reliability  is  obtained  by  further 
additi  ons . Attenuation  curves  are  used  to  adjust  the  peak  acceleration 
magnitude  and  distance  more  closely  to  the  desired  earthquake. 
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Seed,  Ugas  and  Lysmer  (1974)  studied  104  records  broadly  classified 
into  rock  (28  records),  stiff  soil  (31  records),  deep  cohesionless  soil 
(30  records),  and  soft  to  medium  soil  (IS  records).  The  spectra  were 
normalized  to  5%  damping  ratio  and  O.lg  zero-period  acceleration. 
Figures  8-10  and  8-11  show  the  effect  of  site  conditions.  The  soft-to- 
medium  soil  sites  and  deep  cohesionless  soil  sites  show  larger  ampli- 
tudes than  the  others  for  frequencies  less  than  3 Hz.  The  Applied 
Technology  Council  (1976)  have  normalized  these  spectra  as  shown  in 
Figures  8-12  and  8-13.  The  spectra  have  been  compared  to  Regulatory 
Guide  1.60  by  Shannon  and  Wilson  (1975)  as  follows: 

These  comparisons  indicate  that  differences  between  the 
site-dependent  spectra  and  Regulatory  Guide  1.60  are  most 
pronounced  for  frequencies  below  2 to  3 Hz;  less  significant 
but  still  noticeable  differences  are  noted  at  higher  frequen- 
cies. For  rock,  stiff  soil,  or  deep  cohesionless  soil  sites, 
higher  frequency  components  of  composite  spectra  differ  from 
Regulatory  Guide  1.60  by  as  much  as  25%.  Higher  frequency 
components  of  spectra  from  soft-to-medium  soil  sites  exhibit 
larger  differences  relative  to  Regulatory  Guide  1.60.  Figure 
8-14  shows  the  effect  of  shaking  intensity. 

Variability  in  Site-Hatched  Records 

Practical  limitations  in  the  number  of  records  available  limit  the 
selection  process.  Figures  8-15a  and  b show  the  extent  of  variation 
possible.  The  differences  are  probably  caused  by  variations  in  the 
source  mechanism  and  transmission  path.  It  is  only  through  obtaining  a 
number  of  records  that  it  is  possible  to  achieve  sufficient  reliability. 

The  scaling  of  spectral  shapes  implies  that  the  shapes  are  indepen- 
dent of  the  intensity  of  motion.  McGuire  (1977)  has  performed  a regres- 
sion analysis  giving  interesting  information.  It  is  known  that  larger 
magnitude  events  tend  to  produce  relatively  more  or  stronger  long  period 
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motions  than  do  small  magnitude  events.  Further,  the  high  frequency 
components  of  ground  motion  attenuates  with  distance  at  a faster  rate 
than  the  low  frequency  component. 

Figure  8-16  shows  the  mean  psuedo-velocity  spectra  computed  for 
earthquakes  of  several  magnitudes  and  distances,  expressed  as  a percent- 
age of  the  maximum  pseudo-velocity  observed  at  any  frequency  for  that 
value  of  M and  R.  McGuire  (1977)  notes  the  following: 

In  general,  it  is  observed  that  longer  hypocentral  dis- 
tances imply  decreased  relative  spectra  in  the  high  frequency 
region,  and  a lower  frequency  at  which  the  spectral  velocity 
is  a maximum.  Small  magnitude  earthquakes  have  higher  rela- 
tive spectra  in  the  high  frequency  region,  and  lower  relative 
spectra  at  low  frequencies,  than  do  large  earthquakes  observed 
at  the  same  distance. 

Influence  of  Site  Conditions  on  Response  Spectra 

Seed  and  Idriss  (1969)  give  response  spectra  for  sites  whose  soil 
conditions  are  shown  during  the  same  or  similar  magnitude  earthquakes, 
Figure  8-17.  All  of  these  represent  motions  recorded  at  very  consider- 
able distances  from  the  epicentral  regions.  The  spectra  are  arranged  in 
sequence  from  A to  F,  corresponding  to  increasing  degrees  of  "softness" 
of  the  soil  conditions  underlying  the  recording  stations.  To  eliminate 
the  influence  of  different  amplitudes  of  surface  accelerations,  the 
ordinates  of  the  spectra  have  been  normalized  by  dividing  the  spectral 
accelerations  by  the  maximum  ground  acceleration  at  each  site.  Thus  the 
different  forms  of  the  normalized  spectra  reflect  primarily  the  differ- 
ent frequency  characteristics  of  the  motions  from  which  they  were  ob- 
tained. 

It  may  be  seen  that  for  the  recording  made  at  Site  A,  the  peak 
ordinate  of  the  response  spectrum  is  developed  at  a period  of  about  0.3 
second,  indicating  a predominantly  high  frequency  in  the  ground  motion. 
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However,  for  the  still  stiff  (note  high  values  of  standard  penetration 
resistance)  but  slightly  softer  soil  deposit  at  Site  B,  the  peak  ordi- 
nate of  the  spectrum  occurs  at  0.5  second;  and,  as  the  ground  conditions 
become  progressively  softer,  as  evidenced  by  the  presence  of  increasing- 
ly greater  depths  of  soft  and  medium  stiff  clays  and  silts,  the  periods 
at  which  the  peak  spectral  accelerations  are  developed  change  as 
follows : 


Sites 


Approximate  Period 


A 

B 

C 

D 

E 

F 


0.3 


0.5 

0.6 


0.8 


1.3 


2.5 


The  sites  are  arranged  in  increasing  order  of  softness  of  soil  condi- 
tions. The  period  is  that  at  which  maximum  spectral  acceleration  is 
developed.  It  is  thus  apparent  that  the  frequency  components  of  the 
motions  at  the  different  sites  and  the  form  of  the  response  spectra 
change  in  a reasonably  consistent  fashion,  depending  on  the  softness  or 
hardness  of  the  soil  conditions.  The  natural  period  of  structures  in 
seconds  may  be  crudely  estimated  by  dividing  the  number  of  stories  by 
10.0.  For  the  sites  underlain  by  deposits  of  stiff  soils,  the  peak 
ordinates  of  the  acceleration  response  spectra  tend  to  occur  at  a low 
value  of  the  fundamental  period,  say  0.4  or  0.5  second  (see  Figure 
8-17),  indicating  that  at  these  locations  the  maximum  accelerations 
would  be  induced  in  relatively  stiff  structures  five  or  six  stories  in 
height.  On  the  other  hand,  for  the  sites  underlain  by  deep  deposits  of 
softer  soils,  the  peak  ordinates  of  the  acceleration  response  spectra 
tend  to  occur  at  a rather  high  value  of  the  fundamental  period,  say  1.5 
to  2.5  seconds  (see  Figure  8-17),  indicating  that  at  these  sites  the 


stories  in  height.  Thus,  lateral  forces  on  structures  and  related 
building  damage  in  the  same  general  area  may  develop  selectively;  multi- 
story structures  may  be  severely  affected  where  they  rest  on  relatively 
soft  soil  deposits  but  adjacent  stiffer  structures  on  the  same  deposits 
may  be  hardly  affected  at  all;  conversely,  multistory  buildings  on 

shallow,  stiff  soil  deposits  may  be  only  slightly  affected  while  adja- 
cent stiff  structures  are  subjected  to  large  lateral  forces  (Seed  and 
Idriss,  1969). 

Seed  and  Idriss  (1969)  also  demonstrate  the  effect  of  local  and 
distant  earthquakes  of  magnitudes  to  produce  the  same  level  of  site 
ground  motion.  One  example  they  cite  is  as  follows.  A strong  local 

earthquake  may  produce  high  frequency  rock  motions  having  a maximum 
acceleration  of  about  0.13g  and  a predominant  period  of  about  0.25 

second.  On  the  other  hand,  a great  earthquake  located  some  75  miles 

away  might  also  induce  a maximum  acceleration  in  the  rock  of  about 
0.13g,  but  the  motions  would  have  a lower  frequency  and  a somewhat 
longer  predominant  period,  perhaps  of  the  order  of  0.7  second.  The 
computed  ground  surface  accelerations  at  the  site  during  two  such  earth- 
quakes are  listed: 


Maximum 

Acceleration  in 
Underlying  Rock 


Maximum 
Ground  Surface 
Acceleration 


Amplification 

Factor 


e of  Earthquake 


Strong  Local  Shock 
Great  Distant  Earthquake 


The  acceleration  response  spectra  for  these  ground  surface  motions  are 
shown  in  Figure  8-18.  The  maximum  ground  surface  accelerations  differ 
by  some  70%  but  more  significantly,  the  spectral  accelerations  corre- 
sponding to  the  fundamental  period  of  the  building  of  the  site,  about 
1.3  seconds,  differ  by  several  hundred  percent.  Clearly  other  modal 
effects  need  to  be  considered,  but  these  figures  indicate  the  potential 
significance  of  the  characteristics  of  the  basic  rock  motions  in  deter- 
mining ground  surface  and  structural  response  (Seed  and  Idriss,  1969). 


The  influence  of  the  characteristics  of  the  base  rock  motion  on 
ground  response,  as  indicated  by  the  spectra  in  Figure  8-18  is  somewhat 
typical  of  that  found  for  other  site  analyses.  In  some  cases  the  accel- 
eration spectrum  shows  its  major  peak  at  a period  corresponding  approxi- 
mately to  the  fundamental  period  of  the  soil  deposit,  while  in  other 
cases  the  peak  seems  to  develop  at  a period  corresponding  approximately 
with  the  predominant  period  of  the  base  rock  motions.  Which  of  these 
situations  will  develop  in  any  given  case  will  probably  depend  in  large 
measure  on:  (1)  the  difference  between  the  predominant  period  of  the 
base  motion  and  the  natural  period  of  the  deposit  and  (2)  the  duration 
of  the  ground  shaking.  Because  there  will  always  be  some  uncertainty 
concerning  the  form  of  the  base  motions  likely  to  develop  at  any  given 
site,  it  is  clearly  desirable  that  analyses  should  be  made  for  a range 
of  possible,  but  realistic,  base  excitations  to  bracket  the  range  of 
surface  accelerations  to  be  expected.  It  is  clearly  shown  that  because 
the  predominant  period  of  the  rock  motions  for  the  strong  distant  earth- 
quake would  be  larger  than  that  for  the  local  earthquake,  a soft  soil 
deposit  would  tend  to  emphasize  these  long  period  effects,  producing  a 
much  higher  damage  potential  for  multistory  buildings  than  would  the 
shorter  period  motions  of  the  local  earthquake.  Such  a result  would  be 
expected  since  the  large  amplitude  motions  produce  greater  strains  in 
the  transmitting  soil  deposit,  resulting  in  correspondingly  lower  moduli 
and  a longer  fundamental  period  of  the  soil  deposit  (Seed  and  Idriss, 
1969).  Short  period  structures  develop  their  highest  damage  potential 
when  they  are  situated  on  shallow  or  stiff  deposits.  Long  period  struc- 
tures develop  their  highest  damage  potential  when  they  are  situated  on 
deep  or  soft  soil  deposits.  Both  reflect  the  undesirable  effects  of 
similarities  in  the  fundamental  period  of  a structure  and  the  soil 
deposit  on  which  it  is  constructed.  Seed  et  al.  (1975),  Figure  8-19, 
gives  approximate  relationships  between  maximum  accelerations  for  dif- 
ferent soil  types  for  the  same  depth  to  bedrock.  For  low  intensity 
motion  virtually  all  soil  deposits  tend  to  amplify  underlying  rock 
accelerations.  For  extremely  high  intensity  motions  many  soil  deposits 


■ay  attenuate  underlying  rock  acceleration.  Soft  soil  deposits  are  more 
effective  in  attenuating  peak  rock  accelerations  than  shallow  stiff 
soils. 

Gates  (1976)  reports  the  influence  of  soil  depth  (Figure  8-20) 
based  on  analytical  procedures.  A Soil  Amplification  Factor  was  deter- 
mined as  a function  of  acceleration  level  and  site  description.  The 
Soil  Amplification  Factor  is  the  spectral  ratio  between  the  computed 
surface  motions  and  the  input  rock  motion  for  various  frequencies  for  5% 
damping.  The  ratio  is  primarily  dependent  on  the  soil  properties  and  is 
not  influenced  to  a great  extent  by  the  frequency  content  of  the  input 
motion. 


Computer  Analysis  Techniques,  One-Dimensional  Models 


A soil  profile  may  be  analyzed  as  a one-diuiensional  shear  wave 
problem,  assuming  the  stress  wave  to  be  only  a vertically  propagating 
shear  wave  (generally  this  is  true  for  deep  alluvium  sites  and  not  true 
for  shallow  sites  near  the  fault).  The  differential  equations  of  motion 
can  be  solved  in  closed  form  for  linear  elastic  soil  properties.  This 
has  been  done  by  Seed  and  Idriss  (1969)  and  Kanai  (1961)  to  provide  a 
one-dimensional  analysis  of  sites  of  simple  geometry.  However,  the 
stress-strain  characteristics  of  a site  are  highly  nonlinear,  hystere- 
tic,  and  strain-dependent. 

Streeter  et  al.  (1974)  developed  a computer  program  using  the 
method  of  characteristics  for  calculating  one-dimensional  dynamic  behav- 
ior of  soils.  A soil  profile  is  divided  into  layers  down  to  bedrock. 
Dynamic  excitation  of  the  soil  is  introduced  at  the  rock-soil  interface. 
The  response  of  the  soil  can  he  evaluated  on  the  basis  of  elastic, 
viscoelastic,  or  nonlinear  (Ramberg-Osgood)  soil  behavior.  The  program 
determines  shear,  velocity,  and  displacement  information. 

An  analytical  technique  for  analyzing  the  response  of  horizontal 
soil  profiles  to  earthquake  motion  is  described  by  Seed  and  Idriss 
(1969,  1970a,  1970b)  and  Idriss  and  Seed  (1968,  1970).  The  soil  profile 


is  idealized  by  a series  of  discrete  masses  and  springs  with  linear 
viscous  dampers.  The  nonlinear  and  hysteretic  stress-strain  character- 
istics of  the  soil  are  introduced  by  using  an  equivalent  shear  modulus 
and  an  equivalent  viscous  damping  factor  which  can  vary  with  each  layer 
of  soil  profile  and  with  the  strain  level  within  the  layer.  The  equiva- 
lent shear  modulus  for  a given  strain  level  is  taken  as  the  slope  of  the 
diagonal  line  (average  slope)  drawn  through  the  hysteresis  loop,  which 
is  shown  in  Figure  8-21  for  a cyclically  loaded  laboratory  specimen. 
The  average  equivalent  viscous  damping  coefficient  is  proportional  to 
the  ratio  of  area  of  the  hysteretic  loop,  as  shown  in  the  figure,  to  the 
maximum  stored  energy  during  the  cycle. 

An  iterative  procedure  is  used  to  obtain  strain  compatible  values 
for  shear  modulus  and  damping.  The  response  of  the  soil  profile  modeled 
as  discrete  masses  is  computed,  and  strains  are  determined. 

Another  automated-analysis  technique,  more  widely  used  today  for 
treating  horizontal  soil  layers,  has  been  developed  by  Schnabel,  Lysmer, 
and  Seed  (1972),  based  on  the  one-dimensional  wave  propagation  method. 
This  program,  SHAKE,  can  compute  the  responses  for  a given  horizontal 
earthquake  acceleration  specified  anywhere  in  the  system.  The  analysis 
incorporates  nonlinear  soil  behavior,  the  effect  of  the  elasticity  of 
the  base  rock,  and  variable  damping.  It  computes  the  responses  in  a 
system  of  homogeneous  viscoelastic  layers  of  infinite  horizontal  extent, 
subject  to  vertically  traveling  shear  waves.  The  program  is  based  on 
the  continuous  solution  of  the  wave-equation  adapted  for  use  with  tran- 
sient motions  through  the  Fast  Fourier  Transform  algorithm.  Equivalent 
linear  soil  properties  are  obtained  by  an  iterative  procedure  for  values 
of  modulus  and  damping  compatible  with  the  effective  strains  in  each 
layer.  The  following  assumptions  are  made: 

1.  The  soil  layers  extend  infinitely  in  the  horizontal  direction. 

2.  The  layers  are  completely  defined  by  shear  modulus,  critical- 
damping  ratio,  density,  and  thickness. 


3.  The  soil  values  are  independent  of  frequency. 

4.  Only  vertically  propagating,  horizontal  shear  waves  are  con- 
sidered. 

The  soil  model  is  similar  to  that  developed  by  Seed  and  Idriss 
(1970c),  using  data  similar  to  Hardin  and  Drnevich  (1970).  The  absolute 
• range  of  soil  parameter  variation  may  be  stipulated  by  merely  in-putting 

factors  whose  numerical  values  may  be  derived  from  simple  soil  strength 
properties.  These  strength  properties  may  be  the  undrained  shear 
strength  of  a clay  or  the  relative  density  for  sands.  The  program 
requires  the  definition  of  the  soil  profile  down  to  bedrock  (assumed  as 
seismic  velocity  2,500  ft/sec)  as  well  as  an  earthquake  time  history 
record  in  digital  form. 

The  motion  used  as  a basis  for  the  analysis  can  be  given  in  any 
layer  in  the  system,  and  new  motions  can  be  computed  in  any  other  layer. 
Maximum  stresses  and  strains,  as  well  as  time  histories,  may  be  obtained 
in  the  middle  of  each  layer.  Response  spectra  may  be  obtained  and 
amplification  spectra  determined. 

Effects  of  Soil  and  Site  Parameters 

Frequently,  the  parameters  needed  in  the  response  studies  are 
poorly  defined  at  a given  location.  Often,  the  values  of  these  parame- 
ters must  be  assumed  in  order  to  perform  the  ground  response  analyses. 
Experience  has  shown  that  variations  in  the  value  of  any  one  of  the 
parameters  may  affect  the  solution  differently  from  site  to  site,  and  no 
general  rules  may  be  formulated  at  this  time  to  establish  the  influence 
of  the  variables. 

Earthquake  motions  are  produced  by  a stress  wave,  which  is  trans- 
mitted more  rapidly  and  with  less  energy  loss  through  the  bedrock  than 
through  the  overlying  soils.  When  the  bedrock  has  a horizontal  surface 
of  great  extent  and  the  overlying  soil  layers  are  also  horizontal,  it  is 
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frequently  assumed  that  the  earthquake  motion  within  the  soil  is  pro- 
duced essentially  by  horizontal  shear  waves  which  propagate  upward 
through  the  soil  from  the  bedrock  surface.  This  assumption  greatly 
simplifies  the  analysis  since  the  problem  can  be  reduced  to  a one- 
dimensional shear  wave  problem.  This  is  a simplification,  since  verti- 
cal components  of  the  earthquake  motion  are  always  present  and  the  wave 
transmission  problem  may  be  more  complex  than  can  be  simulated  in  a 
one-dimensional  model. 


When  the  bedrock  or  soil  layers  are  inclined,  a one-dimensional 
shear  wave  assumption  is  questionable,  and  a two-dimensional  model  may 
be  required  to  account  for  the  more  complex  geometry  and  wave  motion. 

Lysmer,  Seed,  and  Schnabel  (1970)  have  shown  that  under  identical 
boundary  conditions,  the  lumped  mass  solution  and  the  wave  propagation 
solution  are  basically  the  same.  Arango  and  Dietrich  (1972)  have  inves- 
tigated the  variation  of  parameters  for  the  two  methods.  They  note 
close  agreement  in  peak  levels  of  motion  with  some  differences  in  com- 
puted time  histories. 

Depth  to  Bedrock 

In  many  cases  the  depth  to  bedrock  is  not  well-defined.  A prelim- 
inary analysis  may  be  required  to  assess  the  influence  of  depth  to 
bedrock  on  the  ground  response.  Dezfulian  and  Seed  (1969)  have  shown 
that  an  increase  in  thickness  of  the  deposit  may  or  may  not  cause  a 
substantial  change  in  surface  motion.  Their  studies  show  that  for 
shallow  deposits,  an  increase  in  thickness  of  medium  sand  from  38  to  50 
feet  (12  to  15.6  meters)  reduced  the  response  significantly.  Increasing 
the  thickness  to  80  feet  (25  meters)  reduced  the  response  still  more, 
but  a further  increase  from  80  to  100  feet  (25  to  31.2  meters)  did  not 
produce  any  additional  reduction  in  the  response. 

For  much  deeper  deposits,  1,000  feet  (330  meters),  Kiefer  et  al. 
(1970)  analyzing  the  conditions  at  Osaka,  found  that  the  response  was 
not  very  sensitive  to  the  range  of  depths  investigated. 
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Arango  and  Dietrich  (1972)  studied  the  variation  for  depths  to 
bedrock  equal  to  600,  800,  and  900  feet  (180,  250,  and  270  meters).  The 
values  of  the  maximum  acceleration  and  the  velocity  and  acceleration 
spectra  are  shown  in  Figure  8-22.  The  acceleration  spectra  for  two 
depths  to  bedrock  at  Study  Site  A are  also  shown  in  Figure  8-22. 

The  above  examples  show  that  for  shallow  soil  deposits,  the  depth 
to  bedrock  may  or  may  not  significantly  affect  the  response;  deep  soil 
deposits  are,  in  general,  less  sensitive.  Preliminary  studies  using  a 
reasonable  range  of  depth  to  bedrock  should  precede  any  ground  response 
calculation  when  uncertainties  regarding  the  actual  depth  are  present. 

Influence  of  Soil  Profile 

The  frequency  characteristics  of  the  ground  motions  and  the  form  of 
the  ground  response  spectra  may  be  influenced  by  the  nature  of  the  soil 
conditions  underlying  the  sites.  This  is  illustrated  by  the  studies  by 
Arango  and  Dietrich  (1972).  Different  soil  profiles  were  used  in  the 
response  analysis  as  shown  in  Figure  8-23.  The  values  of  the  maximum 
acceleration  and  displacement  obtained  are  shown  in  Figure  8-24.  The 
corresponding  response  spectra  are  shown  in  Figure  8-25.  Significant 
changes  in  response  can  result  from  variation  in  soil  profile,  and  great 
care  must  be  placed  on  the  correct  site  stratigraphic  representations. 
The  importance  of  the  time  history  of  the  ground  motion  on  the  response 
values  is  also  apparent  by  comparing  the  spectra  from  Figure  8-25a  to 
that  of  Figure  8-25b. 

Soil  Rigidity 

Since  the  stiffness  of  the  soil  deposits  can  only  be  approximated, 
it  is  often  desirable  to  run  preliminary  response  analyses  using  the 
most  reasonable  values  of  the  shear  moduli  for  the  various  soils  and 
values  (say  50%  to  100%)  greater  than  those  judged  to  be  the  most  rea- 
sonable. Arango  and  Dietrich  (1972)  calculated  the  maximum  ground 


surface  acceleration,  ground  displacement,  the  fundamental  period  of  the 
soil  coluarn,  and  the  response  spectra  by  using  the  average  values  of  the 
shear  moduli  and  values  50%  higher.  The  results  of  the  calculations  are 
shown  in  Figure  6-26.  In  some  cases,  errors  in  the  estimated  shear 
moduli  cause  minor  differences  in  the  calculated  ground  response  which 
have  no  practical  significance  for  engineering  purposes.  In  other 
cases,  however,  it  has  been  found  that  great  differences  may  occur  as  a 
result  of  varying  the  values  of  the  shear  moduli. 

Amplitude  of  Rock  Acceleration 

Schnabel  and  Seed  (1972)  have  indicated  that  spectral  acceleration 
values  are  often  not  significantly  influenced  by  substantial  reductions 
in  maximum  acceleration  levels  in  rock.  It  was  found  that  generally  a 
reduction  of  15%  to  25%  in  maximum  rock  acceleration  values  will  affect 
the  spectral  acceleration  by  less  than  10%. 

Frequency  Content  of  the  Rock  Motions 

The  form  and  frequency  characteristics  of  the  base  input  may  have  a 
very  significant  influence  on  the  response  of  soils.  Arango  and 
Dietrich  (1972)  studied  a site  under  two  different  earthquakes  (Figure 
8-27).  As  shown  at  the  right  side  of  Figure  8-27,  the  two  acceleration 
histories  applied  to  the  outcrop  rock  had  the  same  peak  acceleration  and 
the  spectra  were  similar.  However,  the  small  differences  in  frequency 
caused  the  motion  to  be  amplified  differently  in  the  three  soil  columns. 
Therefore,  two  or  more  histories  of  acceleration  should  be  considered  in 
any  given  response  study  in  order  to  define  the  relative  magnitude  of 
the  ground  response  at  any  given  location. 


Computer  Analysis  Techniques,  Two-Dimensional  Models 

As  pointed  out  earlier,  when  the  ground  surface  or  the  soil  layers 
are  inclined,  one-dimensional  wave  assumptions  may  not  be  valid  and  a 
two-dimensional  model  may  be  required  to  represent  the  more  complex 
geometry.  Although  two-dimensional  liquefaction  analyses  are  not  in 
routine  soil  practice,  the  same  procedures  for  evaluation  of  a stress 
history  can  be  utilized.  Finite  element  representations  have  been  used 
to  study  dams  and  embankments. 

Idriss  et  al.  (1973)  have  developed  a two-dimensional  finite  ele- 
ment program  - QUAD-4  - for  the  evaluation  of  seismic  response  of  soil 
deposits.  This  program  allows  for  variable  damping  in  each  element 
using  a Rayleigh  damping  expression  for  that  element.  The  damping 
matrix  for  the  entire  assemblage  of  elements  is  obtained  by  appropriate 
addition  of  the  damping  submatrices  of  all  the  elements. 

The  response  is  evaluated  by  the  solution  of  the  equations  of 
motion  using  direct  numerical  integration  methods  with  a time  increment 
small  enough  to  provide  stability.  The  program  uses  plane  strain  quad- 
rilateral and  triangular  elements.  An  iteration  procedure  is  used  to 
determine  the  strain -dependent  modulus  and  damping  for  each  element, 
based  on  the  average  strain  developed  in  that  element.  The  relation  of 
modulus  and  damping  is  based  on  Seed  and  Idriss  (1970c).  The  solution 
is  obtained  using  the  modulus  and  damping  for  each  element  which  is 
compatible  with  the  average  strain.  The  developers  of  the  program 
report  that  comparison  with  one-dimensional  methods  shows  that  the 
finite  element  solution  values  of  shear  stress  are  about  10%  greater. 
The  response  spectra  of  one-  and  two-dimensional  methods  are  of  similar 
shape.  Major  differences  on  response  spectra  occur  only  when  the  input 
motion  has  large  amounts  of  high  frequency  components  or  when  the  finite 
element  model  is  very  coarse.  The  addition  of  variable  damping  makes 
the  response  calculation  results  in  better  agreement  with  recorded  data. 
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Lysmer  et  al.  (1974)  have  developed  a two-dimensional  finite  ele- 
ment program,  LUSH  (revised  version  called  FLUSH),  which  solves  the 
transient  response  problem  in  soil  sites  by  complex  frequency  response. 
It  can  calculate  the  response  of  sloping  soil  layers  and  can  include  the 
soil-structure  interaction  effect.  The  program  accounts  for  the  non- 
linear effects  which  occur  in  soil  masses  by  a combination  of  the  equiv- 
alent linear  method  described  in  the  section  on  one-dimensional  analyses 
(Seed  and  Idriss,  1969)  and  the  method  of  complex  response  with  complex 
moduli  allowing  for  different  damping  properties  in  all  elements. 

The  model  consists  of  plane  quadrilateral  or  triangular  elements. 
Three  different  material  types  are  provided  for:  nonlinear  clays  and 
sands,  elastic  solids,  and  rigid  solids.  Typical  relationships  between 
stiffness,  damping,  and  effective  shear  strains  for  sand  and  clay  are 
provided  within  the  program.  These  are  similar  to  the  curves  used  in 
SHAKE.  Viscous  damping  is  introduced  by  using  complex  moduli  in  the 
formation  of  the  stiffness  matrix  which  leads  to  the  same  amplitude 

response  as  nodal  analysis  with  a uniform  fraction  of  critical  damping. 

-4 

The  initial  soil  properties  are  specified  at  low  strain  level  (y  = 10  % 
strain)  and  the  program  iterates  to  find  material  properties  at  strain 
levels  compatible  with  the  specified  motion. 

The  mesh  size  of  elements  in  the  model  should  be  small  compared 
with  the  wavelength  of  shear  waves  propagating  through  the  model.  A 
suggested  maximum  height  element  is 
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Case  Study 


Chapter  7 presented  a portion  of  a case  study  as  an  example  of  the 
computational  procedure.  This  section  will  demonstrate  the  preparation 
of  response  spectra.  Assume  it  is  of  interest  to  generate  the  response 
spectra  associated  with  motion  from  the  Santa  Cruz  fault.  Table  7-1 
shows  a 90%  confidence  level  of  0.283g.  Associated  with  the  accelera- 
tion is  a causitive  magnitude  7.2  and  epicentral  distance  25  miles. 
This  information  is  used  as  input  to  Program  0PT1REC.  The  program  using 
the  weighting  specified  by  the  user  (in  this  case  equal  weights), 
searches  the  list  of  records,  and  produces  an  ordered  list  of  records 
representing  the  site  conditions.  Table  8-2  gives  the  first  eight 
entries.  Figure  8-28  shows  time  histories  of  the  records  obtained  by 
Program  T1MH1S.  Figure  8-29a  gives  the  response  spectra  for  5%  damping 
for  each  record.  Figure  8-29b  gives  the  mean  value  and  the  mean  plus 
one  standard  deviation.  Figure  8-29c  gives  the  maximum  envelope  of 
points.  These  were  obtained  using  the  Program  RESPLOT. 

As  can  be  seen  by  Figure  8-29a  there  is  some  variability  in  the 
group  of  records.  Previously,  the  randomness  of  the  site  acceleration 
in  terms  of  attenuation  and  dispersion,  source  description,  epicentral 
location  and  size  of  event  have  been  considered.  This  section  must  only 
consider  site  characteristics  which  may  cause  amplification  of  motion. 
The  best  choice  of  spectra  from  a design  viewpoint  would  probably  be  a 
mean  plus  one  standard  deviation.  Care  must  be  taken  that  the  inclusion 
of  acceleration  dispersion  not  be  taken  twice;  however,  there  is  no 
practical  way  of  separation  of  these  effects.  Superimposed  on  Figure 
8-29b  is  the  response  spectra  which  would  be  indicated  by  the  Newmark- 
Hall  procedure.  Judgment  would  indicate  that  the  mean  plus  one  standard 
deviation  should  be  appropriate  for  design  usage. 

The  results  in  Figure  8-29  were  dominated  by  the  San  Fernando 
earthquake.  Additional  data  were  included  from  Santa  Barbara, 
Long  Beach,  Eureka,  and  Taft  earthquakes.  Results  remained  essentially 


Tabic  8-1.  Structure  Categories  and  Modification  factors 
(from  Applied  Technology  Council,  1974) 


Structural  Categoric*. 

CATEGORY  A 
Space  Frame 

Ductile  Moment  Resisting  Concrete 
(Code  K * 0.67) 


Space  Frame 

Ductile  Moment  Resisting  Steel 
(Code  K • 0.67) 


Ductile  Moment  Resisting  Concrete 
Frame  with  Concrete  Shear  Walls 
(Code  K < 0.80) 


CATEGORY B 

Ductile  Moment  Resisting  Steel 
Frame  with  Vertical  Bracing 
(Code  K > 0.80) 

Concrete  Shear  Walls  with 
Vertical  l.oad-Bcaring  Frame 
(Code  K * 1.00) 


CATEGORY  C 

Vertical  Load  Steel  Frame 
With  Vertical  Bracing 
(Code  K - 1.00) 


Concrete  Sheat  Bearing  Walls 
With  Rigid  Diaphragm 
(Code  K ■ 1.33) 


Flexible  Diaphragms  on 
Masonry  or  Tilt-up  Walls 
(Code  K • 1.33) 


Spectral 

Modification  Factors 


Damage 


Collapse 


4.00  10% 


Damage  Collapse 

**D  Pd  m'c  Pc. 

1.50  7%  2.00  7% 


Damage  Collapse 

r'd  0'd  **’c  Pc 

1.50  5%  2.00  5% 


Project  Building* 

NO.  LATERAL  FORCE  SYSTEM 

1 5 stories 

ductile  reinforced  concrete  frame 

2 19  stories 

ductile  steel  frame  E-W 
(See  also  Category  B) 

3 10  stories 

ductile  reinforced  concrete 
frame  N-S,  ductile  reinforced 
concrete  frame  with  shear  walls  E-W 

10  2 stories 

ductile  steel  frame 


NO.  LATERAL  FORCE  SYSTEM 

2 19  stories 

braced  steel  frame  N-S 
(See  also  Category  A) 

8 9 stories 

reinforced  concrete  frame 
with  shear  walls 


NO.  LATERAL  FORCE  SYSTEM 

4 14  stories 
concrete  shear  walls 

5 7 stories 
concrete  shear  walls 

6 2 stories 

steel  frame  with  vertical  bracing 

7 3 levels 

masonry  shear  walls 

9 1 story 

concrete  tilt-up  shear  walls 

1 1 1 story 

masonry  shear  walls 
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Table  8-2.  Earthquake  Records  From  San  Fernando  Earthquake 


(Magnitude  was  6.6.) 


CEL  Record 

Epicentral 

''max 

Identification 

Distance  (mi) 

(8) 

Percent  of 
Critical  Damping 


Displacement  Velocity  Acceleration 


(a)  Relative  values  of  spectrum  amplification  factors 

(from  "Seismic  design  spectra  for  nuclear  power  plants 
by  N.  M.  Newmark,  J.  A.  Blume,  and  K.  K.  Kapur,  in 
Journal  of  the  Power  Division,  ASCE,  vol  99,  no.  P02, 
Nov  1973,  table  5). 


♦.v 


Ground  Motion  Man  mo 


FrtquMicy , cps 

(b)  Velocity  versus  frequency. 

Figure  8-1.  Basic  design  spectra  normalized  to  l.Og  (from 
Applied  Technology  Council,  1974). 
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Figure  8-2.  Design  spectra  for  inelastic  analysis 

(from  Applied  Technology  Council,  1974). 
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Figure  8-3.  Elastic  response  spectra,  damage  threshold  level 
(from  Applied  Technology  Council,  1974). 
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Figure  8-4.  Elastic  response  spectra,  collapse  threshold  level 
(from  Applied  Technology  Council,  1974). 
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CATEGORY  A 

DTSS- Damage  Threshold  Spectrum  for  Strength  p’j)  = 1.5  = 1046 

CTSD— Collapse  Threshold  Spectrum  for  Deformation  (j'q  = 4.0  = 10% 


Figure  8-5.  Modified  design  spectra  for  Category  A buildings 
(from  Applied  Technology  Council,  1974). 
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Pseudo  Velocity  (in. /sec  ) 


200 


100 

80 

60 


Figure  8-6.  Modified  design  spectra  for  Category  B buildings 
(from  Applied  Technology  Council,  1974). 
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Amplification  factors  for  Control  Points 
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acce lerat ion , and  is  36  in.  for  ground  acceleration  of  1.0  gravity. 
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(b)  Spectrum  shapes 


Regulatory  Guide  1.60  response  spectra — horizontal 
ground  motions  (from  Shannon-Wilson  and  Agbabian 
Associates.  1975). 


Figure  8-8 
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Amplification  Factors  for  Control  Points 
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Total  number  of  records  onolysed  04 


Spectra  tor  5%  dompmg 


■ Soft  to  medium  cloy  and  sand  - 15  records 
/ — Deep  cohesionless  soils  l >250  ftl  - 30  records 


/ / — Stiff  soil  conditions  < < 1 50  ft ) -3i  records 


r- Rock  - 28  records 


Period  - seconds 


(a)  Average  acceleration  spectra  for  different  site  conditions 


Period  - seconds 

(b)  84th  percentile  acceleration  spectra  for  different 
site  conditions. 

Figure  8-10.  Seed-Ugas-Lysmer  spectrum  shapes  (arithmetic 
plot)  (from  Seed,  Ugas,  and  Lysmer,  1974). 
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(b)  Ensemble  MSD 


Seed-Ugas-Lysmer  spectrum  shapes  based  on  local 
soil  conditions  (from  Shannon-Wilson  and  Agbabian 
Associates,  1975). 


Figure  8-11 
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Normalized  spectral  curves  recommended  for  use  in 
building  code  (from  Applied  Technology  Council,  l1 
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Figure  8-13.  Spectral  curves  for  use  in  building  code 

normalized  to  EPA  for  stiff  site  conditions 
(from  Applied  Technology  Council,  1976). 
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Figure  8-14.  Effect  of  shaking  intensity  on  normalized  acceleration  spectra  for 
deep  cohesionless  soil  sites  (from  Seed,  Ugas,  Lysmer,  1974). 
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Comparisons  of  individual  spectra  with  similar 
site  conditions  (from  Shannon-Wilson  and 
Agbabian  Associates,  1975). 
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Figure  8-16.  Predicted  pseudo-velocities  at  different  frequencies 
for  earthquakes  of  several  magnitudes  and  distances, 
normalized  to  maximum  pseudo-velocity  at  any  frequency 
for  that  magnitude  and  distance  (reproduced  with 
permission  from  "Seismic  design  spectra  and  mapping 
procedures  using  hazard  analysis  based  directly  on 
oscillator  response,"  by  R.  K.  McGuire,  Earthquake 
Engineering  and  Structural  Dynamics,  5 211-234,  1977. 
Copyright  © 1977,  John  Wiley  & Sons  Limited). 


Effect  of  soil  conditions  on  form  of  response 
spectra  (from  "Influence  of  soil  conditions  on 
ground  motion  during  earthquakes,"  by  H.  B.  Seed 
and  I.  M.  Idriss,  in  Journal  of  the  Soil  Mechanics 
and  Foundations  Division,  ASCE,  vol  95,  no.  SMI, 
Jan  1969,  figure  7). 


Figure  8-17 
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Acceleration  spectra  for  computed  ground  surface 
accelerations  at  a site  during  two  earthquakes  with 
equal  maximum  base  rock  accelerations  (from 
"Influence  of  soil  conditions  on  ground  motions 
during  earthquakes,"  by  H.  B.  Seed  and  I.  M.  Idriss 
in  Journal  of  the  Soil  Mechanics  and  Foundations 
Division,  ASCE,  vol  95,  SMI,  Jan  1969,  figure  25). 


Figure  8-18 
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(a)  Shallow  stiff 


(b)  Medium 


Soil  amplification  spectral  ratio  (from 
"California  seismic  design  criteria  for 
bridges,"  by  J.  H.  Gates,  in  Journal  of 
the  Structural  Division,  vol  102,  no.  ST12 
Dec  1976,  figure  4). 


Figure  8-20 
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Figure  8-21.  Equivalent  linear  shear  moduli  and  damping 

used  in  discrete  mass  model  (from  "Influence 
of  soil  conditions  on  ground  motions  during 
earthquakes,"  by  H.  B.  Seed  and  I.  M.  Idriss 
in  Journal  of  the  Soil  Mechanics  and 
Foundations  Division,  ASCE,  vol  95,  no.  SMI, 
Jan  1969,  figure  13). 
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Figure  8-23.  Four  soil  profile  representations  of  same  site  (from  "Soil  and  earthquake 
uncertainties  on  site  response  studies,"  by  I.  Arango  and  R.  J.  Dietrich 
[formerly  of  Shannon  & Wilson,  Inc.],  in  Proceedings  of  the  International 
Conference  on  Microzonation  for  Safer  Construction  Research  and  Application, 
30  Oct-3  Nov  1972). 


Figure  8-23.  Continued 
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(a)  El  Centro  earthquake. 
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(b)  Site  intensity  = 8+  earthquake. 


Figure  8-25.  Response  spectra  (from  "Soil  and  earthquake 
uncertainties  on  site  response  studies,"  by 
I.  Arango  and  R.  J.  Dietrich  [formerly  of 
Shannon  & Wilson,  Inc.],  in  Proceedings  of 
the  International  Conference  on  Microzonation 
for  Safer  Construction  Research  and  Application, 
30  Oct-3  Nov  1972) . 
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(a)  Generalized  soil  profile. 


Figure  8-26.  Effect  of  variation  of  material  properties  (from 

"Soil  and  earthquake  uncertainties  on  site  response 
studies,"  by  I.  Arango  and  R.  J.  Dietrich  [formerly 
of  Shannon  & Wilson,  Inc.],  in  Proceedings  of  the 
International  Conference  on  Microzonation  for 
Safer  Construction  Research  and  Application, 

30  Oct-3  Nov  1972). 


180 


Fundamental  Period  (s>  ^ Maximum  Ground  Surface  Acceleration  (g) 


Base  motion  with  a predominant 
period  of  t).  3 5 \econd 

computed  using  shear 
moduli  50%  higher  than 
typical  average  values 


Base  motion  with  a predominant 
period  of  0.35  second 


computed  for  typical 
average  values  of 
shear  moduli 


computed  for  typical  average 
values  of  shear  moduli 


computed  using  shear 
moduli  50%  higher 
than  typical  average 
values 


Maximum  Acceleration  in  Base  Rock  (g) 


) Relationship  between  max- 
imum ground  surface  dis- 
placement and  maxi- 
mum base  rock 
acceleration. 


Maximum  Acceleration  in  Base  Rock  (g) 


(c)  Relationship  between  max- 
imum ground  surface  dis- 
placement and  maximum 
base  rock  accel- 
eration. 


computed  for  typical  average 
values  of  shear  moduli 


Acceleration  Response  Spectra 


M * 6.5  A * 50  kms  (31  miles) 
M » 7.0  d » BO  kms  (50  miles) 
M » 7.5  a » 100  kms  (62  miles) 
base  motion  with  a predominant 
period  of  0.35  second  and  a maxi- 
mum acceleration  of  0.04  g. 


computed  using  shear  moduli  50% 
higher  than  typical  average  values 


base  motion  with  fundamental 
period  >0.35  second 


computed  for  shear  moduli  50% 
higher  than  typical  values 


0.04  0.08  0.12  0.16  0.20  0.24 

Maximum  Acceleration  in  Base  Rock  (g) 


for  average  values 


(d)  Fundamental  period  of 
soil  deposit. 


Period  (s) 


(e)  Effect  of  soil  moduli 
on  the  acceleration 
response  spectra. 


Figure  8-26.  Continued 


(a)  Response  spectra 


Shallow  Section 


(b)  Maximum  surface  accelerations. 
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Figure  8-27.  Response  spectra  and  maximum  surface  accelerations 
(from  "Evaluation  of  soil  liquefaction  effects  on 
level  ground  during  earthquake,"  by  H.  B.  Seed  in 
ASCE  Preprint  2752,  Liquefaction  Problems  in 
Geotechnical  Engineering,  ASCE  Annual  Convention, 
Philadelphia,  Pa.,  27  Sep-1  Oct  1976). 
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Figure  8-28.  Continued 





Figure  8-28.  Continued 
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Figure  8-29.  Response  spectra  with  5%  damping 
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Figure  8-29.  Continued. 
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CHAPTER  9 - DESIGN  LEVELS  AND  ALLOWABLE  RISK 


Risk 

It  is  generally  impossible  to  design  and  construct  a building  with 
a 1001  guarantee  that  the  building  will  not  fail  in  some  way  and  en- 
danger people  as  a result  of  an  earthquake;  the  resources  available  to 
construct  buildings  are  limited.  The  anticipated  benefits  of  various 
life  protecting  programs  must  be  weighed  against  the  cost  of  implement- 
ing such  programs.  Further,  it  is  impossible  to  specify  with  exact 
certainty  the  strongest  earthquake  ground  shaking  that  might  possibly 
occur  at  any  specified  location.  As  long  as  uncertainties  exist  in  the 
loading  it  is  impossible  to  design  for  exactly  zero  risk.  A decision  to 
design  a building  for  a specified  capacity  has  associated  with  it  an 
implicit  risk.  This  implied  risk  may  be  quite  small  (e.g.,  1 chance  in 
10,000  that  a building  will  fail  during  an  earthquake),  but  it  is 
greater  than  zero.  In  general,  losses  may  be  in  the  form  of  damage  and 
repair  costs,  injuries  and  fatalities,  and  the  indirect  adverse  effects 
upon  the  area.  The  emphasis  of  seismic  design  should  be  the  risk  of 
failure  of  buildings , where  such  a failure  would  imply  a threat  to  life 
and  the  risk  of  failure  of  critical  facilities  to  function  at  the  re- 
quired level  to  insure  operational  requirements. 

Both  the  time  of  occurrence  and  the  magnitude  to  affect  a particu- 
lar region  is  unknown.  Future  losses  are  uncertain,  and  some  measure  of 
probability  must  be  used  in  evaluation  of  the  losses. 

One  way  is  the  use  of  average  annual  losses,  which  can  be  expressed 
as  the  average  dollar  losses  per  year,  the  average  major  casualties  per 
year,  the  average  number  of  building  failures  per  year,  etc.  Losses 
expressed  in  this  way  are  annual  risks.  However,  large  earthquakes  are 
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very  rare  events,  and  losses  averaged  over  such  infrequent  events  may 
not  give  a meaningful  portrayal  of  the  large  loss  that  might  occur  over 
one  such  event. 


- 


tJ 


1 
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The  second  way  is  to  define  a threshold  of  loss  and  to  estimate  the 
probability  that  the  threshold  will  be  equalled  or  exceeded  during  some 
earthquake.  For  example,  we  might  speak  of  the  probability  that  the 
dollar  cost  of  damages  and  repairs  will  exceed  $1  billion  during  at 
least  one  earthquake  during  the  next  50  years.  The  threshold  might 
alternatively  be  failure  of  a number  of  critical  facilities. 

Chapter  5 discusses  the  use  of  time  as  a random  variable.  Using  an 
exponential  distribution  the  probability  of  exceedence  may  be  expressed 
as: 


PIE1 


1 


where  A 


average  rate  of  occurrence 


Tf  = average  recurrence  interval  or  return  period 
T^  = exposure  period 


If  AT  » T /T  « 1,  then: 

e r 

P ( E ) 5 A Te/Tr 


The  following  results  are  of  particular  interest: 
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0.5 

39% 

0.2 

18% 

0.1 

10% 

0.05 

5% 
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Note  that  T » T aieana  that  it  ia  quite  (63%)  certain  that  the  event 

will  occur.  If  the  probability  of  the  event  ia  to  be  leaa  than  10%,  Tr 

amat  be  More  than  10T  . 

e 

Acceptable  Kiak 

No  law  in  the  United  Statea  puta  an  aaK>unt  on  the  acceptable  nuaiber 
of  fatalitiea  per  peraon  expoaed  per  year,  or  on  any  other  propoaed 
definition  of  acceptable  riak. 

The  Applied  Technology  Council  (1976)  auawariaea  data  for  the 

probability  of  Man-atade  and  natural  diaaatera  causing  a greater-than- 

apecified  nuaiber  of  fatalitiea  (Figures  9-1  and  9-2).  It  the  total 

"Man-caused"  and  total  "natural"  curvea  are  reduced  by  1,000  (ao  as  to 

give  a level  of  riak  that  would  not  contribute  significantly  to  total 

overall  riak),  the  "acceptable"  rate  of  eventa  cauaing  100  or  wore 

-4 

fatalitiea  would  be  about  5 x 10  With  a 50-year  exposure  period, 
there  would  be  a 2.5%  probability  of  one  or  wore  auch  eventa. 

Considering  all  of  the  relevant  factors,  judgMent  would  indicate 
that  the  following  riak  levels  should  not  tie  exceeded;  that  ia,  riak 
levela  above  those  indicated  should  be  considered  unacceptable. 

1.  The  highest  level  of  safety  attainable  with  present  technology 
ia  required  for  critical  aiiasion  essential  facilities  and  those  contain- 
ing  nuclear  Material  which  could  be  released.  Further,  explosives 
Manufacturing  facilities  with  large  quantities  of  sensitive  Material 
that  could  detonate  should  be  included.  Critical  sea  walls,  failure  of 
which  could  inundate  a base,  Must  also  be  included.  Hedundant  facili- 
ties Must  be  considered  when  a high  en  <gh  level  of  survivability  cannot 
be  assured. 

2.  The  second  category  of  structures  - those  needed  after  a 
disaster  should  also  eMhody  a high  level  of  safety,  although  less  than 
in  the  previous  category.  The  principal  safety  goal  should  be  struc- 
tures and  facilities  that  will  continue  to  function  after  experiencing 
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an  earthquake  disaster.  Exaaples  would  be  sea  walls,  important  utility 
centers,  hospitals,  fire,  police,  and  emergency  communicat ion  facili- 
ties, and  critical  transportation  elements  such  as  bridges  and  over- 
passes. These  should  he  designed  and  built  as  safely  as  Modern  technol- 
ogy permits.  An  extra  project  cost  of  5%  to  25%  to  provide  that  level 
of  safety  would  not  be  excessive.  In  sosw  cases,  where  present  technol- 
ogy cannot  assure  aaxiaua  attainable  safety,  redundance  May  be  neces- 
sary. 


3.  In  a third  category  are  other  structures  that  have  occasional 
high  occupancy  or  are  needed  for  use  after  a disaster.  These  should 
have  as  little  risk  of  collapse  as  is  economically  possible  and  also 
protect  the  occupants  from  serious  hazard  to  the  highest  degree  possi- 
ble. Under  these  criteria,  such  structures  could  be  damaged  to  a point 
where  their  functions  were  impaired;  but  their  continued  functioning, 
though  under  difficult  circumstances,  would  still  be  possible.  Some 
injury  would  possibly  occur  from  falling  articles  and  other  hazards,  but 
loss  of  life  would  be  rare.  Facilities  in  this  risk  category  would 
include  schools,  churches,  theaters,  other  places  normally  attracting 
large  concentrations  of  people,  secondary  utility  structures,  extremely 
large  structures,  roads,  alternate  or  noncritical  bridges  and  over- 
passes. The  degree  of  safety  required  in  this  class  of  structures  would 
be  attainable  by  a 5%  to  15%  increase  in  construction  cost. 


A.  The  vast  majority  of  structures  consist  of  industrial  build- 
ings, quarters,  or  family  housing,  for  which  an  "ordinary"  degree  of 
i'isk  should  be  acceptable.  These  criteria  require  that  buildings  should 
be  able  to: 


a.  Resist  minor  earthquakes  without  damage 


b.  Resist  moderate  earthquakes  without  structural  damage,  but 
with  some  nonstructural  damage 
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II 


c.  Resist  major  earthquakes , of  the  intensity  or  severity  of 
the  strongest  experienced  without  collapse,  but  with  some 
structural  as  well  as  nonst ructural  damage 


In  most  structures,  it  is  expected  that  structural  daaiage,  even  in  a 
major  earthquake,  could  be  limited  to  repairable  damage. 


Building  Failure  and  Total  Risk 


Exceeding  the  design  earthquake  by  itself  does  not  determine  the 
risk  of  failure  of  a structure;  total  risk  is  also  a function  of  the 
inherent  building  strength.  The  seismic  hazard  may  be  represented  by  a 
probability  distribution  of  site  acceleration.  This  must  be  related  to 
the  damage  done  during  each  acceleration  occurrence.  The  probability  of 
facility  damage  may  be  estimated  from  structural  analysis.  Having 
defined  both  the  seismic  hazard  and  the  probability  of  facility  failure, 
the  total  risk  could  be  obtained  by  summing  for  all  accelerations  the 
product  of  probability  of  facility  failure  at  a specific  acceleration 
times  the  number  of  occurrences  of  that  acceleration. 

Whitman  et  al.  (1975)  illustrate  one  approach  for  computing  total 
risk.  They  construct  a damage  matrix  (a  variation  of  which  is  shown  in 
Figure  9-3 a)  in  which  a damage  state  is  shown  as  a function  of  ground 
motion.  Damage  probabilities  are  estimated,  based  on  uncertainties  of 
facility  construction.  For  each  ground  motion  level  individual  damage 
probabilities  are  multiplied  with  their  corresponding  damage  ratio  and 
summed  to  give  the  total  expected  damage  for  that  acceleration.  This  is 
plotted  in  Figure  9-3b.  This  data  can  be  combined  with  the  probability 
that  a site  will  experience  various  ground  acceleration  levels  to  deter- 
mine total  risk  and  damage  for  all  ground  motion  possible.  Thus,  site 
acceleration  may  be  included  by  use  of  the  computational  technique 
described  in  Appendix  B for  determining  site  microzonat ion  analysis. 
One  of  the  output  items  in  this  analysis  technique  is  data  on  accelera- 
tion and  probability  of  not  being  exceeded  in  a number  of  years 
(Figure  9-4). 
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With  Figures  9-3b,  9-4,  and  9-5  it  is  possible  to  examine  minimum 
loss  conditions  for  various  levels  of  earthquake  ground  motion.  The 
methodology  proposed  is  as  follows: 


1.  Select  design  level  acceleration 

2.  Using  Figure  9-5,  estimate  cost  increase 

3.  Estimate  loss  for  design  level 


n 

Loss  = 

i=l 


Probability  site  acceleration  a.  (Figure  4) 
Damage  a.  (Figure  9-3b)  summed  over  n incre- 
ments in  site  acceleration  level 


4.  Compute  total:  Cost  = Loss  + Construction  Cost 

5.  Repeat  for  various  design  levels 

The  result  for  a typical  example  is  shown  in  Figure  9-6.  They 
indicate  the  most  economical  approximation  is  about  95%  confidence  of 
not  exceeding  a site  acceleration  for  the  acceleration  probability  shown 
in  Figure  9-4.  This  is  also  expressed  as  approximately  an  event  with  a 
500-year  return  period.  These  results  are  extremely  crude  using  only 
approximations  for  cost  and  damage.  Further  effort  should  be  spent  in 
this  area  to  obtain  the  maximum  return  and  minimum  risk.  The  optimum 
design  level  will  be  a function  of  site  acceleration.  The  problem  was 
repeated  using  different  site  acceleration  levels.  Figure  9-7  shows 
least  total  cost  (expected  damage  plus  increased  construction  cost) 
return  periods  for  90%  confidence  level  site  acceleration.  Figure  9-7 
clearly  shows  with  site  acceleration  it  becomes  more  economical  to 
accept  greater  risk.  For  example,  if  a 90%  probability  site  accelera- 
tion were  0.2g,  the  least  cost  event  return  period  would  be  about  210 
years.  This  event  could  be  determined  from  the  computer  program  output 
listing.  Economics  is  not  suggested  as  the  criterion  for  design  of 
structures.  However,  any  criteria  proposed  should  be  reviewed  on  an 
economic  basis  as  well  as  a structural  basis. 
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Figure  9-2.  Natural  events  (from  Applied  Technology 
Council,  1970). 
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Figure  9-3.  Determination  of  damage  ratio 
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Figure  9-5.  Acceleration  versus  construction  cost  increase 
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Appendix  A 


NOAA 

PROGRAM  CH42C 


Program  CH42C  in  a CDC  6600  COBOL- FORTRAN  program  used  to  search 
the  National  and  Geophysical  and  Solar-Terrestrial  Data  Center's  INGSDC) 
hypocenter  data  file.  This  data  file  contains  hypocenter  data  from  163R 
to  the  present.  The  tape  records  are  90  characters  long  (punched  card 
images)  blocked  SO. 

There  are  five  types  of  input  cards  (search  set)  needed  to  search 
the  file.  As  many  as  SO  different  searches  may  he  made  during  a single 
pass  through  the  tape.  A separate  printout,  card  deck,  or  tape  file  is 
generated  for  each  of  the  searches.  The  input  cards  must  be  sequenced 
as  follows: 


-This  card  is  used  to  specify  the  types 
of  output  desired  for  each  set. 


use  a 1 if  tape  output  is  desired 
use  a 1 if  a listing  is  desired 
use  a 1 if  punchout  is  desired 
use  a 1 if  plotter  tape  is  desired 


leader  care 


-This  card  is  used  to  print  a header 
on  the  printout  for  each  set. 


CC  1 must  be  blank 

2-80  free  form  (area,  customer  name,  etc.) 


Type  3 (Parameter  cards) These  cards  are  used  to  search  the 

file  for  specific  types  of  events. 

At  least  one  of  the  15  different  types 
of  parameter  cards  (A  thru  P)  must  be 
present  for  each  search  set.  A 
description  of  the  parameter  card  formats 
is  attached. 

Tyj»e  4 i(Special_ Parameter  Combination/"or"  card)---- This  card  is 

used  to  combine  parameter  card  sets  in  an 
"or"  condition  (e.g.,  magnitude  4.0  or 
intensity  TV).  This  card  is  optional. 


CC  1-2 
3-80 


or 

blank 


CC  - Card  Column 
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This  card  is  used  to  end  each  search  set. 


Type  5 (End  card)- 
CC  1-3  ZZZ 


INPUT  TAPE  UNIT:  TAPE  10. 

OUTPUT  TAPE  UNIT:  TAPE  11.  for  block  50 

Parameter  Card  Formats  for  Program  CH42C 

The  parameter  cards  A thru  P described  below  may  be  in  any  se- 
quence. Except  for  card  K (of  which  tnere  can  be  several)  there  can  be 
only  one  of  each  for  each  search  set.  All  parameter  cards  used  in  a 
search  must  be  satisfied  before  a record  is  considered  a hit.  However, 
if  more  than  one  magnitude  card  (C,F,G,H)  is  used,  only  one  of  the 
magnitude  cards  need  be  satisfied  for  a record  to  be  considered  a hit. 
For  those  parameter  cards  with  several  items  (cards  B,I ,K,L,M,N) , a hit 
on  only  one  of  the  items  is  needed  to  satisfy  the  parameter  card. 

Card  A (Interval) 


This  card  is  used  to  search  on  dates,  only  the  low  and  high  year 
must  be  present  on  the  card. 


card 


CC  1 

A 

2 

blank 

3-6 

year  (1964) 

7-8 

month  or  blank 

9-10 

day  or  blank 

11-12 

hour  or  blank 

13-16 

year  (1973) 

17-18 

month  or  blank 

19-20 

day  or  blank 

21-22 

hour  or  blank 

B (Area) 

This  card 

is  used  to  search  for  events  in 

reas.  From  one  to  four  different  areas 

CC  1 

B 

2 

blank 

3-5 

low  latitude,  example  900  = 90 

6-8 

high  latitude 

9 

direction  of  latitude  (N  or  S) 

10 

blank 

11-14 

low  longitude,  example  1500  = 

15-18 

high  longitude 
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19 

direction  of  longitude  (E  or  W) 

20 

blank 

21-23 

low  latitude 

24-26 

high  latitude 

27 

direction  of  latitude 

28 

blank 

29-32 

low  longitude 

33-36 

high  longitude 

37 

direction  of  longitude 

38 

blank 

39-41 

low  latitude 

42-44 

high  latitude 

45 

direction  of  latitude 

46 

blank 

47-50 

low  longitude 

51-54 

high  longitude 

55 

direction  of  longitude 

56 

blank 

57-59 

low  latitude 

60-62 

high  latitude 

63 

direction  of  latitude 

64 

blank 

65-68 

low  longitude 

69-72 

high  longitude 

73 

direction  of  longitude 

Card  C (MB  Magnitude) 

This  card  is  used  to  search  for  events  with  MB  magnitudes  within 
specified  range. 

CC  1 C 

2 blank 

3-4  low  MB  magnitude,  example  55  = 5.5 

5- 10  blank 

11-12  high  MB  magnitude 
Card  D (Depth) 

This  card  is  used  to  search  for  events  with  depth  within  specified 
range.  If  low  and  high  depths  are  blank,  all  tape  records  with  blank 
for  a depth  will  be  considered  hits. 

CC  1 D 

2 blank 

3-5  low  depth  (in  km) 

6- 10  blank 

11-13  high  depth  (in  km) 
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This  card  is  used  to  search  for  events  derived  from  the  specified 
minimum  number  of  stations.  Any  tape  record  whose  number  of  stations  is 
equal  to  or  greater  than  the  number  indicated  on  the  card  is  considered 
a hit. 


2-10  blank 

11-13  number  of  stations 


Card  F (MS  Magnitudes 


This  card  is  used  to  search  for  events  with  MS  magnitudes  within 
the  specified  range. 


blank 

low  MS  magnitude,  example  55 
blank 

high  MS  magnitude 


Card  G (Magnitudes  From  Other  Sources 


This  card  is  used  to  search  for  events  with  magnitudes  from  other 
sources  within  the  specified  range. 


blank 

low  magnitude,  example  500 
blank 

high  magnitude 


Card  H (Local  Magnitudes 


This  card  is  used  to  search  for  events  with  local  magnitudes  within 
the  specified  range. 


blank 

low  local  magnitude,  example  500 
blank 

high  local  magnitude 


Card  I (Associated  Phenomena 


This  card  is  used  to  search  for  events  with  various  associated 
phenomena.  Whatever  is  in  CC  2 indicates  the  field  to  be  searched. 


I 


"D"  diastrophic  phenomena 

"T"  tsunamis 

"S"  seiches 

"V"  volcanism 

"N"  non-teetonic  events 

"0"  infrasonic  waves 


Whatever  is  in  CC  11-17  indicates  what  that  field  is  being  checked 
against.  Therefore,  this  card  allows  search  of  one  of  six  different 
fields  and  allows  search  of  that  field  for  from  one  to  seven  different 
i terns . 


CC  1 
2 

3-10 

11-17 


Diastrophic 


I 

field  to  be  searched  either  D,  T,  S,  V,  N,  or  0 
blank 

item  to  be  searched  for.  If  search  is  made  on  only  one 


item  put  symbol  in  CC  11. 
are  as  follows: 


F = surface  faulting 

U = uplift/subsidence 

D = faulting  and  uplift/subsidence 

Non-Tectonic 

R = rockburst 
C = coal  bump  or  rockburst 
in  coal  mine 
M = meteoritic  source 
E = explosion  or  suspected  explosion 
I = collapse  . 

L = lights  or  other  visual  phenomena 


Items  which  can  be  selected 


Tsunami 

T = tsunami  generated 
7 = possible  tsunami 


Seiche 

S = 

? = 


seiche  generated 
possible  seiche 


Volcanism 

V = earthquake  associated 
with  volcanism 


Infrasonic 

T = T-wsve 
A = acoustic  wave 
G = gravity  wave 
B — both  A & G 

Card  J (Intensities) 

This  card  is  used  to  search  for  events  with  intensities  within  the 
specified  range.  Any  tape  record  whose  intensity  is  equal  to  or  greater 
than  the  intensity  indicated  on  the  card  is  considered  a hit. 
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CC  1 
2 
3 


J 

blank 

intensity  (1  thru  9 
T = XII) 


= I thru  IX,  X = X,  E = XI, 


Card  K (Flinn-Engdahl  Regions) 


This  card  is  to  search  for  events  falling  within  the  specified 
Flinn-Engdahl  geographic  region.  Several  K cards  can  be  used;  however, 
not  more  than  120  different  Flinn-Engdahl  geographic  numbers  may  be 
used. 


CC  1 
2 

3-5 

6-80 


K 

blank 

first  Flinn-Engdahl  geographic  number 

next  25  Flinn-Engdahl  geographic  numbers  (acceptable  to 
have  many  K cards  with  less  than  the  26  F-E  numbers) 


Card  L (Sources) 


This  card  is  used  to  search  for  events  derived  by  specific  sources. 
From  1 to  13  comparisons  can  be  made. 


CC  1 L 

2 blank 

3-5  first  source 

9-41  next  12  sources 


Card  M (Cultural  Effects) 


This  card  is  used  to  search  for  events  with  cultural  effects, 
one  to  four  different  cultural  effects  may  be  searched. 


From 


CC  1 
2 


3 

4 

5 


M 

blank 

first  cultural  effect 
second  cultural  effect 
third  cultural  effect;  C 
F = felt 

fourth  cultural  effect 


= casualties,  D = damage, 


Card  N (Depth  Control) 


This  card  is  used  to  search  for  events  with  specific  Depth  Control 
Designations.  From  one  to  four  different  depth  control  designations  may 
be  searched. 


CC  1 
2 
3 


N 

blank 

first  depth  control  designation 
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4 second  depth  control  designation 

5 third  depth  control  designation 

6 fourth  depth  control  designation 

N = held  at  33  km 

G = depth  restrained  by  geophysicist 

D = restrained  depth  based  on  2 or  more  reported  pP's 

A = assigned 

Card  P (Radius  Card) 


This  card  is  used  to  search  for  events  within  the  specified  radius 
from  the  specified  geographic  point. 


CC  1 
2 

3-6 

7 

8 

9-13 

14 

15 

16-20 

21 


P 

blank 

latitude,  example  8900  = 89.00° 

latitude  direction  (N  or  S) 
blank 

longitude,  example  10,000  = 100.00° 

longitude  direction  (E  or  W) 
filler 

radius  in  km 

blank  = radius  in  km 

M = radius  in  miles 


Type  5 (End  Card) 


This  card  is  used  to  end  each  search  set. 


Appendix  B 
CEL 

PROGRAM  RECAL 


ii 

i 


This  program  uses  a subset  of  data  obtained  from  the  NOAA  histori- 
cal data  (see  Appendix  A)  as  input.  The  program  computes  and  plots  the 
annual  recurrence  interval  as  a function  of  magnitude  for  the  region. 
It  then,  using  this  data  and  assuming  a Poisson  distribution,  computes 
the  cumulative  probability  of  not  exceeding  a magnitude  over  a period  of 
time.  Faults  may  be  defined  as  line  segments,  and  all  events  within  a 
specified  distance  are  associated  with  that  fault;  recurrence  data  and 
cumulative  probability  are  computed  as  above  for  each  fault  input.  The 
program  has  options  to  omit  or  include  regional  analysis  and  to  calcu- 
late an  input  fault  recurrence.  Using  the  fault  recurrence,  the  program 
(by  use  of  a Monte  Carlo  technique)  assigns  2,000  years  of  earthquakes 
random  in  magnitude  and  epicentrally  located  to  the  fault.  The  shortest 
distance  from  the  fault  break  to  the  site  is  calculated  and  an  accelera- 
tion (mean  and  standard  deviation)  is  determined.  A Poisson  distribu- 
tion is  used  to  account  for  acceleration  uncertainty.  Where  faults  are 
not  clearly  defined  a floating  earthquake  may  be  assigned  for  the  re- 
gion. 

In  assigning  line  segments  to  a fault  one  should  take  care  to 
assign  all  events  in  the  area  to  specific  faults;  thus,  events  will  not 
be  overlooked  or  recurrence  interval  reduced.  Uncertainties  exist  in 
the  actual  location  of  many  epicenters.  Further  discussion  of  this  is 
in  the  text. 

It  is  intended  that  the  historic  recurrence  data  calculated  be 
adjusted  based  on  judgment,  taking  into  account  recurrence  intervals 
estimated  by  other  techniques  such  as  slip  rate,  etc.  Provisions  are 
included  to  allow  for  input  of  recurrence  coefficients. 

Since  the  program  uses  random  number  generation  and  the  Monte  Carlo 
technique,  there  will  be  some  variation  in  solution,  particularly  for 
values  of  confidence  above  70%.  However,  below  this  number,  sufficient 
samples  have  been  taken  to  insure  good  repeatability. 

User's  Guide 


Card  1 

Col 

Format 

Col 

(110  7F10.0) 

1 

10 

NPT 

number  points  in  data  set  if  not  specified 
seek  blank  card  after  earthquake  data 

11 

20 

SEALEX 

scale  for  plotting  inches  per  degree 
longitude 

21 

30 

SCALEY 

scale  for  plotting  inches  per  degree 
latitude 

31 

40 

XMAX 

left  side  maximum  longitude 

41 

50 

XMIN 

right  side  minimum  longitude  first  plot 

223 

1 1 

J L 


51 

60 

1/MAX 

top  maximum  latitude 

61 

70 

1/MIN 

bottom  minimum  latitude 

71 

80 

XMINN 

last  longitude  covered  for  multiple  plots 

Card 

2 

Format 

(8  F10.0) 

1 

10 

AMIN 

minimum  magnitude  default  = 0 

11 

20 

AMAX 

maximum  magnitude  default  = 8.6 

21 

30 

TIME 

time  interval  for  probability 
default  = 25.0  years 

31 

40 

SITEX 

site  longitude 

41 

50 

SITEY 

site  latitude 

51 

60 

FAC 

plot  scale  miles  per  inch 

61 

70 

R 

distance  from  fault  for  inclusion 
default  = 10  miles 

71 

75 

0PT1 

option  to  omit  regional  analysis 

76 

80 

0PT3 

option  to  perform  regional  analysis 
(floating  earthquake) 

Card 

_3_ 

Format 

(7  F8.0,  4X, 

2A10) 

1 

8 

west mo st  fault  longitude 

9 

16 

westmost  fault  latitude 

17 

24 

center  fault  longitude 

25 

32 

center  fault  latitude 

33 

40 

eastmost  fault  longitude 

41 

48 

eastmost  fault  latitude 

49 

56 

maximum  fault  magnitude 

57 

60 

option  to  read  A and  B from  next  card 

61 

80 

fault  name 

(Repeat  Card  3 as  required;  end  group  with  blank  card.) 

Data  cards  conpatible  with  output  from  Program  CH42. 

Data 

Cards  (All  fixed  point) 

Col 

Col 

5 

6 

first  two  digits  of  year 

11 

12 

second  two  digits  of  year 

20 

25 

latitude 

26 

32 

longitude 

36 

38 

body  wave  magnitude 

54 

55 

surface  wave  magnitude 

61 

63 

other  magnitude 

73 

75 

local  magnitude 

NOTE 

l 

If  NPT 

not  specified 

, supply  one  blank  card  after  data. 

Appendix  C 
CEL 

PROGRAM  OPTREC 


This  program  uses  the  data  in  Tables  C-l  through  C-4  in  selecting 
the  closest  records  to  a prescribed  magnitude,  epjcentral  distance, 
acceleration,  and  site  category.  Weighting  factors  are  used  to  combine 
the  effects  of  magnitude,  distance,  and  acceleration.  By  assigning  zero 
weight  to  a variable,  that  variable  may  be  excluded,  thus  it  is  possible 
to  search  by  magnitude  or  distance  or  acceleration  or  any  combination  of 


these. 

Default 

parameters 

consider  equal  weight. 

User' s 

Guide 

Group 

Data  Set 

Col 

Col 

1 

5 

ID 

record  1 

6 

12 

DiST 

epicentral  distance,  miles 

13 

19 

DUR 

record  duration,  sec 

20 

24 

NUM 

CALTEC  record  number 

25 

29 

DIR 

direction,  degrees 

30 

37 

ACC 

acceleration  peak,  cra/sq  sec 

38 

45 

VEL 

velocity  peak,  cm/sec 

46 

53 

D1SP 

displacement  peak,  cm 

54 

63 

11,12 

reference 

64 

68 

AMAG 

magnitude 

69 

70 

ISITE 

site  category 

0 unclassified 

1 alluvium 

2 intermediate 

3 hard  rock 

(Repeat  for 

all.  End 

with  blank  card. ) 

Group  2 

Col 

Col 

1 

10 

SM 

magnitude 

11 

20 

SD 

site  distance,  miles 

21 

30 

SA 

acceleration,  g's 

31 

35 

1 CODE ( 1 ) 

1 = include  unclassified  site 
0 = omit 

36 

40 

1 CODE ( 2 ) 

1 = include  alluvium  sites; 

Li 
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41 

45 

IC0DE(3) 

1 = include  intenaediate  sites; 
0 = omit 

46 

50 

IC0DE(4) 

1 = include  hard  rock  sites; 

0 = oaiit 

51 

60 

WH 

weighting  for  nagnitude 

61 

70 

WD 

weighting  for  distance 

71 

80 

WA 

weighting  for  acceleration 

Oaitting  all  weighting  defaults  to  equal  weighting,  set  any  weight  to 
0.0  to  oaiit  from  weighting. 
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Table  C-l.  Unclassified  Sites 
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1* 

nr$T 
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Nun  ni» 

»rc 

VII 

0 t PP 

»»r 

T 

• 5.PP 

TV. IP 

5 

IPO 

•00,50 

-0 , >5 

>.6P 

1 

1 

7.7 

• 

•5.55 

77, Bp 

5 

>70 

•5> .10 

5.07 

• >,PP 

1 

1 

».7 

IP 

>5.  PI 

5# . ap 

• 

?1 

155.0'' 

•15.70 

• 6 , 7 0 

1 

1 

7.7 

1 1 

?5.P1 

5u.an 

• 

1 1 1 

170.00 

•17,70 

• 5.15 

1 

l 

».7 

1* 

• T ,55 

• ?.#o 

0 

1 PO 

•5a  .10 

-0.1? 

•5.11 

1 

l 

7.7 

IT 

• T .55 

• >.50 

0 

• 0 

#5.50 

5.5P 

• 5 , •• 

1 

1 

7.7 

IP 

• 7.55 

TP. 00 

7 

( PO 

•5* . 1 0 

• 6.57 

• a ,•• 

1 

l 

7.7 

>0 

• T .55 

TP. OP 

7 

•0 

• 1 .?o 

• . 5 1 

6,05 

1 

1 

7.7 

>2 

15. OP 

77, *0 

P 

5<IP 

105.00 

•51,00 

1>.0" 

1 

1 

6.5 

M 

15.00 

75.00 

P 

7* 

•>55.00 

?*.«0 

1 • . 1 0 

1 

1 

6.5 

5 1 

75. ?6 

• 0,00 

1 1 

1*0 
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1 
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1 
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•1 1 .PO 
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1 ,65 

1 

1 

6.0 

55 
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1? 

>70 

•15, BO 

? . • 7 

• >.50 

1 

1 

6. a 

57 

in. 27 

5P.00 

15 

• 5 

•#5,*0 

->.•• 

• 1.10 

1 

1 

5.5 

O 

1 0.2T 

55,00 

15 
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•*«.•« 

a. 06 
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1 

1 

5.5 

a* 
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55.70 

1« 

551 
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• i,5l 

1 

1 

5.5 

« 1 

5.25 

5*  .oo 

1 a 

• 1 

•«5.#0 

• >,i> 

1 ,00 

1 

1 

5.5 

1 S.Pi 

• o.po 

1 p 

1*1 

• 05. «o 

• 7.77 

->.Pa 

1 

1 

5.6 

M 

1 5. 01 

00,50 

1* 

?7  1 

•170.00 

17.10 

5, Pa 

1 

1 

5.6 

55 

• 5.P« 

• T,«0 

IP 
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1 

1 

6.5 

•6 
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1 

6.5 

61 
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>1 
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-15,50 

1 

> 

6.5 

6? 

100. T? 

1 >0.00 

?1 
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.1 5>.op 
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-17.50 

1 

? 

6.5 

67 

>5.6* 

75.00 

>5 

0 

5> , 1 0 

-1  .56 

".76 

1 

? 

5.a 

6P 

>1.M 

C. 

c 

• 
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r- 

? 5 

>70 

•>o.«o 

->.?1 

.a« 

1 

? 

5.0 

70 

'll  .10 

• 0,50 

>tt 

1*0 

• 157, Oi> 

•?0 , 50 

• a,>r 

1 

2 

6.5 

»1 

<11,10 

00,50 

?0 

>70 

•170,00 

1 1 .60 

• 5,66 

1 

> 

6.5 

7* 

•?.!• 
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?5 

( PO 
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1 

? 

6.0 
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•>.lo 

51 ,00 

?5 
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? 
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« 
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>• 
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1 

? 

7.1 
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?• 
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•17.10 

10. ao 

1 

2 

7.1 

HP 

56.77 
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• 0 

55 
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1 

2 

6.5 
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2 
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Table  C-l.  Continued 
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710 

1*9 

•07.10 

-?.*  7 

• 1.10 

1 

8? 

s.o 

• l 9 

11.17 

7.19 

711 

1«« 

•0* , 00 

8.12 

.*7 

1 

89 

4.0 

917 

11.17 

7.19 

771 

90 

• 1*  . 9 0 

-1.11 

.«? 

1 

88 

4.0 

• 88 

8.87 

1*.T" 

710 

11* 

11*. "0 

• *,** 

•9.81 

1 

87 

S.o 

• 88 

8.87 

19.7ft 

710 

80* 

190,00 

•9.60 

1 . *0 

1 

97 

S.4 

97a 

18.87 

89, *0 

71* 

ft 

-110.00 

-4.7* 

• 1 ,7*. 

1 

87 

*.4 

918 

18.81 

89.80 

11* 

*0 

•S? . *0 

•1,10 

1 .** 

1 

91 

S . 4 

• 17 

19.81 

98,80 

119 

1*0 

•90,80 

-?.** 

-.9* 

1 

87 

S.4 

•18 

19. SI 

98.00 

719 

• 0 

•7*.9P 

-1  .77 

• , 70 

1 

91 

*.4 

9*8 

19.79 

87.1* 

198 

rt 

•t*.0" 

I.S« 

1.7* 

1 

93 

5.4 

• *1 

19.79 

87.90 

198 

9* 

1 *.70 

1 .9* 

1.19 

1 

81 

*.4 

• •8 

1*  • *9 

?1.ao 

loo 

9* 

• 10, *0 

!.*1 

-1  ."7 

1 

81 

S.4 

9*7 

18  . *9 

27.90 

190 

1** 

89.10 

-8.0* 

9,1* 

1 

81 

S.4 

•8? 

*0.98 

*9.90 

770 

1*0 

•81  .*" 

-7. SI 

«.2S 

7 

8* 

*.4 

•*1 

• 0.98 

*4,8* 

170 

*8 

?*,1" 

>.71 

•2.11 

1 

8* 

*.« 

•88 

1*7.7? 

*8.70 

771 

17* 

•11.8* 

-«.!• 

7.97 

1 

9* 

*.4 

98* 

1*7,7? 

• 8.1* 

171 

?** 

1 1 .70 

-9.8* 

-9.** 

1 

?S 

*.« 

•88 

1»7. 89 

*1  .*0 

IT? 

319 

9.79 

7.19 

•9.9* 

1 

8* 

9.4 

989 

187.8* 

SI  .** 

77? 

?09 

• . * 1 

-?.«T 

•8.18 

1 

?S 

*.4 

9*1 

117.00 

98.90 

777 

9* 

7.7  * 

-1  .1* 

• .si 

1 

2* 

*.4 

9*? 

117.00 

98.9(1 

171 

1M 

-7,01 

• 1.18 

-.7* 

1 

?S 

9.4 

9*7 

118.9* 

*9.1* 

77* 

0 

9 , *1 

-8.80 

1.71 

1 

2* 

9.4 

9*8 

118.9* 

SI  .*0 

77* 

90 

• 1 ft, 7* 

-8.8S 

• 1 .99 

1 

85 

9.4 

•78 

171 .98 

*«.io 

77* 

1*0 

• 7,0ft 

8.1* 

9,  *1 

1 

85 

9.4 

• 71 

111.98 

*8.10 

77* 

870 

1*,1* 

-2.9* 

•1  .*? 

1 

?S 

9.4 

979 

178.8* 

*8.1* 

17* 

877 

1 • .SO 

9.87 

•?.** 

1 

?S 

9.4 
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?«» 

>5.38 

58.20 

38 

Table  C-2 . 

12?  -216.00 

Cont inued 

-21.80  -11,20 

1 

8 

8.8 

?•! 

>5.33 

38,20 

3* 

217 

-132,00 

•18,80 

11.50 

1 

8 

8,8 

101 

83.33 

3.5? 

101 

1*0 

17.30 

-2.62 

1 .08 

1 

8 

8.8 

10? 

83.38 

1 1 .20 

101 

30 

•1".00 

2.21 

-1.28 

1 

8 

8.8 

107 

>8.03 

27.10 

101 

0 

-31 ,30 

8.81 

-2.51 

1 

8 

8.8 

10» 

>8.03 

>7,00 

101 

270 

121.00 

-3,00 

2.00 

1 

8 

8.8 

m 

21.21 

81,80 

10* 

1*0 

83,10 

-".11 

0.05 

1 

8 

8,8 

It  o 

>3.?1 

81.80 

1 05 

30 

77.80 

8,08 

• 0,95 

1 

8 

8,8 

31* 

>3.88 

2*  ,80 

107 

0 

31.30 

-7.9T 

2.98 

1 

7 

8,8 

J>0 

>1.83 

>3,8« 

107 

«0 

-1 07,00 

10.10 

-7.17 

1 

7 

8,8 

3?? 

>1.88 

33,20 

1 0* 

0 

-198,00 

• 9,80 

2.7? 

1 

7 

8,8 

3*3 

>3.38 

33,30 

1 0* 

30 

-1*2.00 

-18.00 

-8,90 

1 

7 

8.8 

3341 

>3.50 

32.00 

1 1 > 

1* 

102.00 

17.0(1 

11,00 

1 

7 

8.8 

335 

>5.50 

31  .30 

112 

V>3 

78,30 

-13.70 

•9.25 

1 

7 

8,8 

3410 

?o.?7 

37,70 

1 18 

120 

111.00 

10.20 

-1.80 

1 

7 

8,8 

344  1 

20  . ?7 

37.70 

118 

210 

118.00 

-3.15 

2.78 

1 

7 

8.8 

3413 

17.83 

80.10 

118 

1 1 

>21.00 

•?*,?0 

-11.80 

1 

8 

8.8 

34141 

17.81 

00.10 

113 

2*1 

-188,00 

-21.50 

-10,10 

1 

8,8 

33? 

10.1* 

38,10 

M« 

113 

13,70 

•1 1 .30 

8,81 

1 

8 

8,8 

333 

10.1* 

88,10 

11* 

2?3 

-12.70 

9.08 

7. *2 

1 

8 

8,8 

331 

>5.31 

03.70 

121 

270 

-1 1*,00 

-17.20 

• 8.88 

1 

8 

8.8 

S3? 

>5.31 

03,70 

121 

1*0 

11?.oo 

•10,80 

0.00 

1 

8 

8.6 

170 

83.33 

11.00 

1?o 

27o 

*10,30 

0,02 

>.11 

1 

3 

8,6 

331 

83.38 

11.10 

120 

1 80 

•10.50 

-3.80 

-2.71 

1 

6.6 

333 

>?.?! 

2*. 00 

1 ?• 

0 

-80,30 

11.20 

•7.28 

1 

9 

6,6 

330 

>?.?  1 

>7.80 

1 ?• 

270 

•31,80 

15.00 

*.10 

1 

9 

6.6 

3*3 

>2,87 

83.10 

1 H 

30 

-188.00 

17. 20 

9.2* 

1 

9 

6.8 

333 

>2.87 

O*,10 

Hi 

5»0 

-181.00 

• 10.10 

8.11 

1 

9 

8.6 

3*3 

>1.28 

83.30 

1 1« 

30 

-37,30 

18,70 

1 1 .00 

1 

9 

6.6 

333 

>3.23 

03.50 

1 10 

1 «0 

-82.10 

•10,70 

-8.21 

1 

9 

8,6 

303 

17.25 

38,80 

117 

33 

180,00 

-18.10 

-7.10 

1 

9 

8.6 

«0T 

17.23 

58,80 

117 

1*3 

•123,00 

22.10 

-*.05 

1 

9 

6,6 

3*3 

>0.*> 

33.80 

1 03 

1*0 

110,00 

11  .80 

-17.80 

1 

10 

6.6 

3*3 

?0.3> 

33.80 

1 03 

270 

101,00 

•28.80 

15.10 

1 

1 0 

6.8 

333 

>1.3> 

1 » , 80 

1 0* 

0 

-108,00 

18,20 

7.1? 

1 

10 

6.6 

337 

21.9? 

18,80 

1 8* 

270 

112.00 

17.50 

-11.10 

1 

1 0 

6.6 

837 

>3.71 

*3.>0 

178 

17 

•81.00 

20,90 

11.80 

1 

11 

6.6 

333 

>5.71 

*3.20 

178 

1?7 

•1 18.00 

-17,80 

•11,70 

1 

11 

6.8 

833 

31  .58 

103,00 

1 SO 

1 *0 

21.3" 

-8.71 

-1.35 

1 

1 1 

6.6 

500 

31  .54 

103,00 

1*0 

270 

23.90 

-*,08 

8,55 

1 

tl 

6.6 

547 

12.78 

33,70 

188 

181 

33.70 

-3.78 

• 0 ,9(| 

1 

1 0 

6,6 

513 

12.78 

33.70 

188 

231 

-38.70 

-9,70 

5.0? 

1 

1 0 

6.8 

5?0 

<18. ?7 

10.00 

1*7 

1 3 

-53.70 

• 1.08 

-.71 

1 

i a 

6,6 

5?1 

08.27 

10.00 

187 

2*5 

73,30 

-1.87 

-.78 

1 

1 4| 

6.6 

5*3 

*3.3> 

83,  >0 

1 88 

50 

-110.00 

17.00 

10.80 

1 

1 0 

6.6 

5*8 

>1.3? 

83,00 

188 

120 

127.00 

1».?0 

5.01 

1 

1 0 

6,6 

588 

75.37 

•3.80 

133 

11 

-00.90 

-1,38 

-2.19 

1 

1 0 

6.6 

585 

73.17 

33.80 

133 

101 

-11.00 

0.80 

-2.05 

1 

10 

6.6 

587 

83.37 

31.10 

138 

2*0 

13,00 

9,09 

-7.97 

1 

10 

6,6 

583 

03.37 

31.10 

138 

1 30 

-11.?o 

-9  , ?8 

8,70 

1 

1 0 

6.6 

550 

118.83 

82.30 

137 

05 

-23.80 

• 2.19 

-1.18 

1 

1 0 

6.6 

551 

118.83 

82.30 

137 

115 

•13.00 

-2.8? 

•1  .00 

1 

10 

6.6 

553 

>3.18 

83,80 

1 33 

?8 

118.00 

•17,80 

9.78 

1 

18 

6.6 

537 

>5.13 

03,80 

1 33 

238 

219.00 

21.10 

to, 10 

1 

13 

6.6 

533 

85.00 

8*. 70 

?oa 

0 

-28,00 

-8,17 

-5.8? 

1 

15 

8.6 

533 

83.00 

83,70 

?0« 

30 

20.80 

9.33 

-7.27 

1 

13 

6,6 

538 

88.83 

1 1?.oo 

>03 

113 

-2«.oo 

-7,1* 

8.(10 

1 

15 

6,6 

583 

00.83 

112,00 

>03 

?4IO 

-28.10 

10,00 

-8.71 

1 

13 

6,6 

571 

88 . 8 > 

3>. *0 

>08 

0 

17,00 

1,08 

•1.11 

1 

15 

6.6 

57? 

83.82 

3>,30 

208 

• 0 

-01,90 

2. *7 

-1,08 

1 

13 

6,6 
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3«0 

• 5.5* 

53.50 

>10 

135 

34. *0 

4.4* 

1 ,4* 

1 

13 

4.4 

1 

Ml 

• 5.55 

53. *0 

41  * 

445 

34.4" 

• 4. *4 

• 1.34 

1 

13 

*.4 

1 

M< 

45.** 

*1  .*0 

41* 

1 4 ft 

•104.0" 

•14,30 

• *,*a 

1 

1* 

4.4 

1 

Mk 

*10 

U:U 

*1  .*0 
54.10 

*0 

140 

•44.4" 

•?3.*o 

•w; 

*.<)• 

• 4.3* 

1 

1 

1*4 

t:t  1 

41 1 

•*,45 

54.10 

444 

4*0 

43.40 

5.31 

•4. *4 

1 

14 

4.4 

41  It 

51.51 

• 4,40 

431 

0 

• 1.3" 

10, *0 

■ 0.40 

1 

14 

4,4 

41  V 

31.  *1 

• 4,40 

431 

4*0 

3*. 40 

13,30 

10.30 

1 

14 

4,4 

41? 

1 * . 44 

34.50 

433 

1*4 

•443.00 

31  .40 

14.30 

1 

1* 

4,4 

4?) 

1 *.** 

54.50 

433 

4*4 

15*. 00 

• 1 *, 30 

*.«* 

1 

1* 

4,4 

411 

40.** 

• 4.3" 

434 

1*0 

•10*. 00 

•13,30 

•4,14 

1 

1* 

4,4 

41? 

40.T5 

• 4.30 

• 0 

144.00 

10,30 

3.4* 

1 

1* 

4.4 

440 

44. •• 

44, *0 

43* 

1«0 

1 t*.0" 

• 1 * « 40 

• * ,40 

1 

1* 

4.4 

44| 

44.** 

• • ,*0 

43* 

• 0 

1*4,00 

• 1*,10 

•1 1,40 

1 

1* 

4.4 

4*4 

45,04 

41  ,*0 

444 

30* 

14*. 00 

14.40 

5.00 

1 

10 

4.4 

*3* 

45.04 

41  ,*0 

44* 

41* 

•14*. 00 

•1*,*0 

*.*• 

1 

14 

4.4 

441 

41.40 

• 3.50 

444 

1*0 

•114.00 

•14. TO 

4.30 

1 

10 

4.4 

44? 

41.40 

43,50 

4*4 

• 0 

•10*. 00 

•14,30 

•10,30 

1 

10 

4.4 

44T 

41  .44 

• 5.00 

44* 

140 

144.0" 

• 1 • , *0 

*.*• 

1 

10 

4.4 

444 

>1.44 

• 5.00 

44* 

• 0 

•1*4.00 

■14,40 

•1 ",4fl 

1 

10 

4.4 

4*0 

>5.50 

• 1.40 

44* 

• 4 

•**,•" 

14.40 

11.30 

1 

10 

4.4 

4*1 

45.50 

41,40 

44* 

134 

•44.40 

10.10 

-*.33 

1 

to 

4.4 

4*4 

45.45 

34.40 

453 

330 

4*4.00 

1*,30 

ll.*0 

1 

10 

4.4 

4*5 

>5.45 

34.40 

453 

440 

•441 .00 

14.10 

•14.50 

1 

10 

4.4 

44* 

>4.10 

• 4.10 

4*4 

4* 

•34.40 

1 *.4o 

• 1 ",  40 

1 

1* 

4.4 

4*4 

44.10 

44,10 

45* 

11* 

•43,40 

•1*,50 

•10,30 

1 

14 

4.4 

*04 

>5.** 

34. *0 

441 

5* 

•**,40 

•14,30 

•14.40 

1 

1* 

4,4 

TO* 

>5. *« 

34.50 

441 

34* 

•104.0" 

11.40 

3. *5 

1 

1* 

4.4 

*41 

45. •• 

41  .*0 

444 

0 

134,0" 

•1  * . 30 

4.04 

1 

1* 

4,4 

*4? 

>5.** 

• 4,00 

444 

4*0 

•150,00 

41  .30 

•1 1 .40 

1 

1* 

4,4 

*>4 

51. *1 

44. *0 

44* 

0 

53.30 

•13,30 

0.44 

1 

1* 

4.4 

*45 

51  ,«1 

• 4»*0 

44* 

• 0 

41.30 

•13,00 

*.3* 

1 

1* 

4,4 

*1* 

*4.5* 

*4.40 

A** 

313 

3*.30 

4,04 

• ,*0 

1 

41 

3.* 

4io 

*4.54 

*4.40 

4*4 

445 

33. *0 

4. *4 

.*• 

1 

41 

5.* 

445 

• 0.** 

4*  « *0 

300 

31* 

•1 1*,00 

4.54 

• 4«*0 

1 

41 

4.4 

44* 

«0.*T 

4*  ,*0 

500 

445 

114,0" 

• 5. *3 

4.34 

1 

41 

4.4 

444 

4*  , 1 • 

41  ,30 

313 

4*1 

•13.1" 

4,40 

•4.4* 

1 

>1 

3.4 

445 

4*.  1* 

41,30 

313 

141 

14.3" 

1 , *4 

4.00 

1 

41 

3.4 

4** 

14.  *5 

44.00 

31* 

1*0 

•1*.*0 

1.33 

."3 

1 

44 

*.« 

• *5 

14.53 

44.4" 

31* 

440 

n.A" 

*.*3 

.*• 

1 

44 

3.4 

44* 

10, *0 

• 4,*0 

344 

•5 

4,3" 

• ,03 

• ,40 

1 

44 

*.4 

4*0 

10. *0 

• 4,*0 

315 

•44.4" 

4,41 

1.1* 

1 

44 

4,4 

*10 

3. •* 

4*. 00 

34* 

140 

■1*4,00 

-1*.*0 

4.03 

1 

44 

*.* 

*1 1 

5.** 

44, *0 

54* 

4*0 

44.40 

0.0* 

4.44 

1 

44 

«.* 

*1* 

•«.|4 

• 4,*0 

334 

1 «o 

14.30 

• 1.** 

• , *4 

1 

44 

4.3 

•40 

•*.14 

• 4.50 

334 

•0 

•14,30 

• 1 , *3 

•.*3 

1 

44 

4.3 

**S 

134.14 

• 3.40 

3"4 

*.** 

• 4.54 

1 .*0 

1 

43 

4.4 

**• 

134.14 

• 3,40 

5** 

414 

•1 1 .*0 

• 3,0* 

4.34 

1 

43 

4.4 

*** 

154.14 

41.50 

3*4 

144 

4.** 

4.4* 

• 1.0* 

1 

4.4 

*** 

134.14 

>1  .50 

3*4 

414 

•11,30 

3.00 

4.31 

1 

4.4 

231 
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1ft 

Ml? 

nua 

Wu“  ft|3 

arc 

vet 

ft  1 3> 

•22  •*» 

*3 

7.13 

33,40 

13 

10 

-41 ,40 

• 4,4? 

-2.23 

1 

1 

3.3  2 

aa 

?.I3 

33.40 

14 

100 

-103,00 

-a.41 

-.•1 

1 

1 

3.)  2 

as 

4.4a 

00,40 

14 

171 

-43,40 

• 9,03 

• 1,1a 

1 

1 

3.3  2 

at 

4.4a 

ao.?o 

14 

241 

39.10 

a.oa 

.•? 

1 

1 

3.)  2 

aa 

13.03 

ao.30 

1? 

24 

34,00 

1 .33 

• 1 ,40 

1 

1 

3.3  2 

30 

13.03 

ao.io 

1? 

114 

23.70 

1.13 

-1.1) 

1 

1 

3.1  2 

H» 

aa.3? 

30,00 

3* 

3?a 

•la.oo 

• 1 .04 

1 .2? 

1 

2 

3.4  > 

Ml 

aa  .3? 

?3 ,30 

3* 

23a 

1 1 .ao 

.40 

• ,34 

1 

? 

3.4  2 

M3 

31  .aa 

30,10 

3* 

144 

•?o,ao 

• ?.)? 

• ,44 

1 

3.4  2 

ll( 

31. aa 

30.10 

34 

74 

13,30 

-2.42 

1 ,)4 

1 

2 

140 

?3.13 

47.30 

3a 

304 

ia?,oo 

17.40 

1 1 .40 

1 

) 

4.4  2 

»3l 

>3.1 4 

3?,30 

30 

?ta 

•117,00 

- 1 7 • 30 

1 1 .*0 

1 

) 

4.4  2 

IK 

14.33 

41  .40 

34 

21 

-304,00 

•17.20 

«.?) 

1 

a 

4.4  ? 

14? 

14.3* 

41  .30 

34 

2*1 

-243 ,00 

•?7.40 

4,44 

1 

a 

4.4  2 

?U 

>3. 4a 

33.40 

?> 

243 

-42.20 

20.40 

ia,?o 

1 

9 

^ • 4 P 

219 

>3. aa 

3?.?0 

?? 

13 

114.00 

21  .30 

1 1 .30 

1 

3 

4.4  2 

?v 

>3.3? 

34.40 
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Appendix  D 


CEL 

PROGRAM  RESPLOT 


This  program  uses  the  CALTEC  spectra  data  and  plot  Fourier  and 
response  spectra.  The  program  can  produce  SCR  4020  plots  or  CALCOMP 
plots.  Options  exist  in  the  program  to  plot  different  records  on  the 
same  plot,  to  average  records  or  to  plot  an  envelope.  Data  auy  be 
displayed  in  acceleration-period  plots  or  tripartile  plots. 

User's  Guide 


Card  1 

Col  Col 

1 10 


Card  2 


0 CALCOMP  plot 

1 SCR  4020  plot 


Col 

1 

2 

6 


Col 


5 

10 


L = last  of  group  for  plot 
S = stop  after  plot 
ID  file  (see  Appendix  C) 
code 

1 Fourier  spectra 

2 0%  damping  response  spectra 

3 2%  damping  response  spectra 

4 5%  damping  response  spectra 

5 10%  damping  response  spectra 

6 20%  damping  response  spectra 


11 

20 

scale  value  (default  = 1.0) 

21 

30 

response  acceleration  (scaling),  g's 

31 

40 

ground  motion  acceleration  (scaling),  g's 

41 

45 

0 = plot  acceleration-period;  1 = omit 

46 

50 

0 = plot  tripartile;  1 = omit 

51 

55 

1 = average  group;  0 = omit 

56 

60 

1 = maximum  group;  0 = omit 

(Repeat  as  required.) 
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Appendix  E 


CEL 

PROGRAM  TIMHIS 


This  program  scales  plots  and  punches  time 
California  Institute  of  Technology  accelerograms. 


histories  of  the 


User's  Guide 


Card  1 


Col 

Col 

1 

5 

ID  number 

10 

20 

SCALE;  default  = 1.0 

21 

30 

AMAX,  g's 

31 

35 

flag  to  plot  CALC0MP 

35 

40 

flag  to  PUNCH 

41 

45 

flag  to  plot  4020 

46 

80 

format  for  punch  output 
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NAVFACENGCOM  - l.ANT  DIV.;  Fur.  BR  Deputy  Dir.  Naples  Italy;  RDTAELO 00P2.  Norfolk  VA 
N A V F ACENGCOM  - NORTH  DIV.  (Boretsky)  Philadelphia,  PA;  Code  WP(l  CDR  A.J.  Stewart);  Code  1028. 

RDTAELO.  Philadelphia  PA;  Design  Div.  (R.  Masino).  Philadelphia  PA;  ROICC.  Contracts,  Crane  IN 
NAVFACENGCOM  - PAC  DIV.  Code  402.  RDTAE,  Pearl  Harbor  HI:  Commander.  Pearl  Harbor.  HI 
NAVFACENGCOM  - SOUTH  DIV.  Code  SO,  RDTAELO.  Charleston  SC;  Dir.,  New  Orleans  l.A 
NAVFACENGCXIM  - WEST  DIV.  Code  04B;  RDTAELO  Code  201 1 San  Bruno.  CA 

NAVFACENGCOM  CONTRACT  AROICC,  Point  Mugu  CA;  AROICC.  Quantico.  VA;  Eng  Div  dir.  Southwest  Pac. 
Manila.  PI;  OICC.  Southwest  Pac.  Manila.  PI;  OICC/ROICC.  Balboa  Canal  Zone;  ROICC  (Ervin)  Puget  Sound 
Nuval  Shipyard.  Bremerton.  WA,  ROICC  AF  Guam:  ROICC  l.ANT  DIV.,  Norfolk  VA;  ROICC,  Diego  Garcia 
Island;  ROICC,  Kcfluvik.  Iceland.  ROICC',  Pacific.  San  Bruno  C'A 
NAVHOSP  l.T  R.  Elsbernd,  Puerto  Rico 
NAVMACi  SCE.  Guam 
NAVMIROOIC.  Philadelphia  PA 

NAVCX'EANSYSCEN  Code  4(N  (D.  G.  Moore).  San  Diego  CA;  Code  531  l(T)  (E.  Hamilton)  San  Diego  CA;  Code 
6565  (Tech.  Lib.).  San  Diego  CA 
NAVORDSTA  PWO.  Louisville  KY 
NAVPETOFF Code  30.  Alexandria  VA 
NAVPGSOOl.  Code  61  Wl.  (().  Wilson)  Monterey  CA 

NAVPHIBASE  CO.  ACB  2 Norfolk.  VA;  C'ode  S3T,  Norfolk  VA;  Harbor  Clearance  Unit  Two.  Little  Creek.  VA; 

OIC.  UCTONE  Norfolk.  Va 

NAVREGMEDCEN  Code  3041.  Memphis,  Millington TN;  SCE(D.  Kaye);  SCE(l.CDR  B.  E.  Thurston).  Sun  Diego 
CA;  SCE.  Camp  Pendleton  C'A 

NAVSC'CM.C’ECOFF C35  Port  Hueneme,  CA;  C44A  (R.  Chittenden).  Port  Hueneme  C'A;  C'O,  Code  C44A  Pori 
Hueneme.  CA 

NAVSFASYSCOM  Code  OCR’  (l.T  R MacDougal).  Washington  DC 
NAVSEC Code  6034  ( Library).  Washington  IX' 

NAVSF.CGRUACT  PWO.  Torn  Sta.  Okinawa 
NAVSHIPRFPFAC  Library  . Guam;  SCE  Subic  Bay 

NAVSHIPYOCode  202.4,  Long  Beach  C'A;  C'ode  202.5  (Library)  Puget  Sound.  Bremerton  WA;  C'ode  400.  Puget 
Sound.  Code  404  ( l.T  J.  Riccio).  Norfolk.  Portsmouth  VA;  C'ode  410.  Mare  Is..  Vallejo  CA;  Code  440  Portsmouth 
NH;  Code  440.  Norfolk;  Code  440.  Puget  Sound.  Bremerton  WA;  C'ode  440.4,  Charleston  SC';  L.D.  Vivian; 

I ibrary.  Portsmouth  NH;  PWO,  Mare  Is.;  Tech  Library.  Vallejo,  C'A 
NAVSTA  CX)  Naval  Stalmn.  May-port  Fl.;  CO  Roosevelt  Roads  P.R.  Puerto  Rico;  F.ngr.  Dir.,  Rota  Spuin;  Maint.  Div. 
Dir/Code  531.  Rodman  Canal  Zone;  PWD(l.TJG,P.M.  Motolenich),  Puerto  Rico;  PWO  Midway  Island;  PWO. 
Kellavik  Iceland;  PWO.  Mayport  FL;  ROICC.  Rotu  Spain;  SCE.  Guam;  SCE.  Subic  Bay.  R.P.:  Utilities  Engr  Off. 
(l.TJG  A.S.  Ritchie).  Rotu  Spuin 
NAVSUBASE  LTJG  D.W  Peck.  Groton,  CT 

NAVSUPPACT  C'O,  Scuttle  WA;  Code  413  Seattle  WA;  Engr  Div.  (F.  Mollica),  Nuples  Italy  ; l.TJG  McGamth. 
Vallejo  C'A 

NAVSURFWPNCEN  PWO.  White  Oak.  Silver  Spring.  MD 
NAVTF.C'HTRACEN  SCE.  Pensacola  FL 

NAVWPNCEN  Code  2636  (W  Bonner),  China  Lake  C'A;  PWO  (Code  26).  China  lake  CA;  ROICC  (Code  702),  (Tuna 
Lake  CA 
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NAVWPNSTA  EARl.E  ENS  G.A  Lowry,  Fallbrook  CA;  PW  Office  (Code  09CT)  Yorktown,  VA 
NAVWPNSUPPCEN  Code  09  (Bocnmghausen)  Crane  IN 
N<  blJ  405  OIC.  San  Diego.  CA 

NCBC  CELfCAPT  N.  W.  Petersen),  Port  Hueneme.  CA:  CEL  AOIC  Port  Huencme  CA;  Code  10  Davisville.  RL 
Code  155.  Port  Hueneme  CA;  Code  15ft.  Port  Hueneme.  CA;  PW  Engrg.  Gulfport  MS;  PWC)(Code  80)  Port 
Hueneme.  CA 

NCBU4IIOIC.  Norfolk  V A 
NCR  20,  Commander 

NMCB  1 33 (ENS  T.W  Nielsen);  5.  Operations  Dept.;  Eorty.  CO;  THREE.  Operations  Off. 

NORDA  Ctxle  440  (Ocean  Rsch  Off)  Bay  St.  Louts  MS 

NRL  Cixle  8400(J.  Walsh).  Washington  DC;  Code  8441  (R.A.  Skop).  Washington  IX';  Rosenthal.  Ctxle  8440,  Wash. 
DC 

NSD  SCE.  Subic  Bay.  R.P. 

NTC  Ctxle  54(ENS  P.  G.  Jackel).  Orlando  FL;  Commander  Orlando.  FI. 

NUSC  Ctxle  131  New  London,  CT;  Ctxle  EAI 23 ( R.S.  Munn).  New  London  CT;  Ctxle  TAI3I  (G.  De  la  Cruz).  New 
London  CT 

ONR  Ctxle  700F  Arlington  VA;  Dr.  A.  Laufer.  Pasadena  CA 
PHIBCB  I P&E,  Coronado.  CA 

PMTC  Ctxle  4253-3.  Point  Mugu.  CA;  Pat.  Counsel.  Point  Mugu  CA 

PW C ACE  Office  (LTJG  St.  Germain)  Norfolk  VA;  CO  Norfolk,  VA;  CO.  Great  Lakes  IL;  Ctxle  I Ift  (LTJG.  A. 
Eckhart)  Great  Lakes.  IL;  Ctxle  120.  Oakland  CA:  Ctxle  I20C  (Library)  San  Diego.  CA;  Ctxle  128.  Guam:  Ctxle 
200.  Great  l akes  IL;  Code  200.  Oakland  CA;  Ctxle  220  Oakland.  CA;  Code  220. 1 . Norfolk  VA;  Ctxle  30C 
(Boettcher)  San  Diego.  CA;  Ctxle  4(8).  Pearl  Harbor,  HI;  Ctxle  ft80.  San  Diego  CA;  Library.  Subic  Bay.  R.P.;  OIC 
CBU-405.  San  Diego  C'A;  XO  Oakland.  CA 
SPCC  Ctxle  I22B.  Mechanicsburg.  PA:  PWO  (Ctxle  120)  Mechanicsburg  PA 
U.S.  MERCHANT  MARINE  ACADEMY  Kings  Point.  NY  (Reprint  Custodian) 

US  DEPT  OF  INTERIOR  Bureau  of  Land  MNGMNT  - Ctxle  733  (T.E.  Sullivan)  Wash.  DC 

US  GEOlXXilCAL  SURVEY  Off.  Marine  Geology.  Piteleki.  Resttm  VA 

USAF  SCH(X)L  OF  AEROSPACE  MEDICINE  Hyperbaric  Medicine  Div.  Brooks  AFB.  TX 

USCG  (G-ECV)  Washington  IX;  (G-ECV/ftl)  (Burkhart)  Washington.  DC;  G-EOE-4/61  (T.  Dowd).  Washington  DC 

USCG  ACADEMY  LT  N.  Stramandi.  New  London  CT 

USCG  R&D  CENTER  D.  Motherway.  Groton  CT;  Tech.  Dir.  Groton.  CT 

USNA  Ch.  Mech.  Engr.  Dept  Annapolis  MD;  Ocean  Sys.  Eng  Dept  (Dr.  Monney)  Annapolis.  MD;  PWD  Engr.  Div. 

(C.  Bradford)  Annapolis  MD;  PWO  Annapolis  MD 
AMERICAN  CONCRETE  INSTITUTE  Detroit  Ml  (Library) 

CALIF.  DEPT  OF  NAVIGATION  A OCEAN  DEV.  Sacramento.  CA  (G.  Armstrong) 

CALIFORNIA  STATE  UNIVERSITY  LONG  BEACH.  CA  (CHELAPATI);  LONG  BEACH.  CA  (YEN) 
COLORADO  STATE  UNIV..  FOOTHILL  CAMPUS  Fort  Collins  (Nelson) 

CORNELL  UNIVERSITY  Ithaca  NY  (Serials  Dept.  Engr  Lib.) 

DAMES  A M(X)RE  LIBRARY  LOS  ANGELES.  CA 

DUKE  UNIV  MEDICAL  CENTER  B.  Muga.  Durham  NC;  DURHAM.  NC  (VESIC) 

FLORIDA  ATLANTIC  UNIVERSITY  BOCA  RATON.  FL  (MC  ALL1STER);  Boca  Raton  FL  (Ocean  Engr  Dept..  C. 
Lin) 

FLORIDA  ATLANTIC  UNIVERSITY  Boca  Raton  FL  (W.  Tessin) 

FLORIDA  TECHNOLOGICAL  UNIVERSITY  ORLANDO.  FL  (HARTMAN) 

GEORGIA  INSTITUTE  OF  TECHNOLOGY  Atlanta  GA  (School  of  Civil  Engr..  Kahn);  Atlanta  GA  (B.  Mazanti) 
IOWA  STATE  UNIVERSITY  Ames  IA  (CE  Dept.  Handy) 

VIRGINIA  INST.  OF  MARINE  SCI.  Gloucester  Point  VA  (Library) 

LEHIGH  UNIVERSITY  BETHLEHEM.  PA  (MARINE  GEOTECHNICAL  LAB..  RICHARDS);  Bethlehem  PA 
(Fritz  Engr.  Lab  No.  13,  Beedle);  Bethlehem  PA  (I.inderman  Lib.  No. 30.  Flecksteiner) 

LIBRARY  OF  CONGRESS  WASHINGTON,  DC  (SCIENCES  A TECH  DIV) 

MICHIGAN  TECHNOLOGICAL  UNIVERSITY  Houghton.  Ml  (Haas) 

MIT  Cambridge  MA:  Cambridge  MA(Rm  10-500,  Tech.  Reports.  Engr.  Lib.);  Cambridge  M A (Whitman) 

NEW  MEXICO  SOLAR  ENERGY  INST.  Dr.  Zwibel  Las  Cruces  NM 

NORTHWESTERN  UNIV  Z.P.  Bazant  Evanston  IL 

NY  CITY  COMMUNITY  COLLEGE  BROOKLYN.  NY  (LIBRARY) 

UNIV.  NOTRE  DAME  Katona.  Notre  Dame.  IN 

OREGON  STATE  UNIVERSITY  (CE  Dept  Grace)  Corvallis.  OR:  CORVALLIS.  OR(CE  DEPT.  BELL); 
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CORVALLIS,  OR  (CE  DEPT.  HICKS);  Corvalis  OR  (School  of  Oceanography) 

PENNSYLVANIA  STATE  UNIVERSITY  UNIVERSITY  PARK.  PA  (GOTOLSKI) 

PURDUE  UNIVERSITY  Lafayette  IN  (Leonards):  Lafayette.  IN  ( Allschueffl);  Lafayette.  IN  (CE  Engr.  Lib) 

SAN  DIEGO  STATE  UNIV.  I.  Nooruny  San  Diego.  CA;  Dr.  Krishnamoonhy,  San  Diego  CA 
SEATTLE  U Prof  Schwacglcr  Seattle  WA 

SOUTHWEST  RSCH  INST  King.  San  Antonio.  TX;  R DeHart.  San  Antonio  TX 
STANFORD  UNIVERSITY  Engr  Lib.  Stanford  CA;  Stanford  CA  (Gene) 

STATE  UNIV  OF  NEW  YORK  Buffalo.  NY 

IFXAS  A AM  UNIVERSITY  COLLEGE  STATION.  TX  (CE  DEPT);  College  Station  TX(CF.  Dept.  Herbich) 
UNIVERSITY  OF  CALIFORNIA  BERKELEY.  CA  (CE  DEPT.  GERWICK);  BERKELEY.  CA  (CE  DEFT. 
MITCHEl  l ).  Berkeley  CA  (B.  Bresler);  Berkeley  CA  (Dept  of  Naval  Arch.);  Berkeley  CA  (R.  Williamson); 
DAVIS.  CA  (CE  DEPT.  TAYLOR);  LIVERMORE.  CA  (LAWRENCE  LIVEI  MORE  LAB.  TOKARZ);  La  Jolla 
CA  ( Act|.  Dept.  I ib.  C-075A);  M Duncan.  Berkeley  CA;  SAN  DIEGO.  CA.  LA  JOLLA.  CA  (SEROCKI) 
UNIVERSITY  OF  DELAWARE  Newark.  DE  (Dept  of  Civil  Engineering.  Chesson) 

UNIVERSITY  OF  HAWAII  HONOLULU.  HI  (SCIENCE  ANDTECH.  DIV);  Honolulu  HI  (Dr  S/dard) 
UNIVERSITY  OF  ILLINOIS  Met/  Ref  Rm.  Urbana  IL;  URBANA.  IL  (DAVISSON);  URBANA,  II  (LIBRARY); 

URBANA.  IL  (NEWARK);  Urbana  11  (CF.  Dept.  W.  Gamble) 

UNIVERSITY  OF  MASSACHUSETTS  (Heronemus).  Amherst  MA  CE  Dept 
UNIVERSITY  OF  MICHIGAN  Ann  Arbor  Ml  (Richart) 

UNIVERSITY  OF  NEBRASKA-LINCOLN  Lincoln.  NE  (Ross  Ice  Shelf  Proj.) 

UNIVERSITY  OE  NEW  MEXICO  J Nielson-Engr  Mails  A Civil  Sys  Div.  Albuquerque  NM 
UTTVKRSITY  OF  TEXAS  Inst.  Marine  Sci  (Library).  Port  Arkansas  TX 

UNIVERSITY  OE  TEXAS  AT  AUSTIN  AUSTIN.  TX  (THOMPSON);  Austin  TX  (R  Olson);  Austin.  TX  (Breen) 
UNIVERSITY  OE  W ASHINGTON  Dept  of  Civil  F.ngr  (Dr.  Mattock).  Seattle  WA;  SEATTLE.  WA  (MERCHANT); 
SEATTLE.  WA  (OCEAN  ENG  RSCH  LAB.  GRAY);  Seattle  WA  (F..  l inger);  Seattle.  WA  Transportation. 
Construction  A Geom.  Div 
ALFRED  A.  YEE  A ASSOC.  Honolulu  HI 
AMETEK  Offshore  Res.  A Engr  Div 
APPLIED  TECH  COUNCIL  R.  Scholl.  Palo  Alto  CA 
ARVID  GRANT  OLYMPIA.  W A 
ATLANTIC  RICHFIELD  CO.  DALLAS.  TX  (SMITH) 

AUSTRALIA  Dept.  PW  (A.  Hicks).  Melbourne 
BECHTEL  CORP.  SAN  FRANCISCO,  CA  (PHELPS) 

BELGIUM  HAECON.  N.V..  Gent 
BETHLEHEM  STEEL  CO.  Dismuke.  Bethelehem,  PA 
BROWN  A CALDWELL  E M Saunders  Walnut  Creek,  CA 
BROWN  A R(X)T  Houston  TX  (D.  Ward) 

CANADA  Mem  Univ  Newfoundland  (Chari).  St  Johns;  Surveyor,  Nenninger  A Chenevert  Inc..  Montreal;  Warnock 
Herscy  Prof.  Srv  Ltd.  La  Sale.  Quebec 
CF  BRAUN  CO  Du  Bouchel.  Murray  Hill.  NJ 

CHEVRON  OIL  FIELD  RESEARCH  CO.  LA  HABRA.  CA  (BROOKS) 

CONCRETE  TECHNOLOGY  CORP.  TACOMA.  WA  (ANDERSON) 

CONRAD  ASSOC.  Van  Nuys  CA  ( A.  Luisoni) 

DRAVO  CORP  Pittsburgh  PA  (Giannino);  Pittsburgh  PA  (Wright) 

NORWAY  DETNORSKE  VERITAS  (Library).  Oslo 
EVALUATION  ASSOC.  INC  KING  OF  PRUSSIA.  PA  (FEDELE) 

FORD.  BACON  A DAVIS.  INC.  New  York  (l.ibraiy) 

FRANCE  Dr.  Dutertrc.  Boulogne;  L.  Pliskin,  Paris;  P.  Jensen.  Boulogne 
GEOTECHNICAL  ENGINEERS  INC.  Wtnchesler.  MA  (Paulding) 

GLIDDENCO.  STRONGSVILLE.  OH  (RSCH  LIB) 

GRUMMAN  AEROSPACE  CORP.  Belhpuge  NY  (Tech.  Info.  Ctr) 

HALEY  A ALDRICH.  INC.  Cambridge  MA  (Aldrich.  Jr.) 

HONEYWELL.  INC.  Minneapolis  MN  (Residential  Engr  Lib.) 

HUGHES  AIRCRAFT  Culver  City  CA  (Tech.  Doc.  Ctr) 

ITALY  M.  Caironi,  Milan:  Sergio  Tattoni  Milams;  Torino  (F.  Levi) 

MAKAI  OCEAN  ENGRNG  INC.  Kailua.  HI 
JAMES  CO.  R.  Girdley.  Orlando  FI. 

I.AMONT-DOHERTY  GEOLOGICAL  OBSERV.  Palisades  NY  (McCoy):  Palisades  NY  (Selwyn) 


I OCKHKED  MISSILES  * SPACE  CO.  INC.  Mgr  Naval  Arch  A Mar  Eng  Sunnyvale.  CA;  Sunnyvale  CA 
(Rvnewtcr).  Sunnyvale.  CA  (Phillip*) 

LOCKHEED  OCEAN  l ABORATORY  San  DiegoCA<F.  Simpson) 

MARATHON  Oil  CO  Houston  TX  (C  Seay) 

MC  O KI  LAND  ENGINEERS  INC  Houston  TX  (B  McClelland) 

MCDONNEl  AIRCRAFT  CO.  Dept  501  (R  H.  Kay  man).  St  Louts  MO 
ME  DAI  l A ASSOC.  INC  J.T.  GAFFEY  II  SANTA  ANA.  CA 
MEXICO  R Cardenas 

MOBII  PIPE  l INK  CO.  DALLAS.  TX  MGR  OK  ENGR  (NOACK) 

MUESKR.  RUTLEDGE.  WENTWORTH  AND  JOHNSTON  NEW  YORK  (RICHARIXS) 

NEW  / E ALAN D New  Zealand  Concrete  Research  Assoc  (Librarian).  Porirua 
NEW  PORT  NEWS  SHIPBLDG  A DRYDOCK  CO.  Newport  News  VAtTech.  Lib.) 

NORW  AY  I'L  l NORSK!  VERITAS  (Roren)  t)slo;  I.  Foss.  Osh>;  J.  Creed.  Ski;  Norwegian  Tech  Univ  (Brandt/aeg). 
Trondheim 

OFFSHORE  DEVELOPMENT  ENG.  INC.  BERKELEY.  CA 
PACIFIC  MARINE  TECHNOLOGY  Long  Beach.  CA  (Wagner) 

FOR  IT  AND  CEMENT  ASStX'  SKOKIE.  IL  (CORELY);  Skokie  IL(RschA  Dev  lab.  Lib.) 

PRESCON  CtWP  TOW  SON,  MD( KELLER) 

RANDCORP  Santa  Monica  CA  (A.  I aupa) 

RAYMOND  INTERNATIONAL  INC.  E Colle  Soil  Tech  Dept,  Pennsauken.  NJ 
RIVERSIDE  CEMENT  CO  Riverside  CA(W  Smith) 

SANDIA  l ABORAIOklES  Library  IXv.,  Livermore  CA 

sch u pack  assoc  so  norw  alk.  ctischupack) 

SEATFCHCORP  MIAMI.  FT  (PERONI) 

SHEl  I DEVELOPMENT  CO.  Houston  TX(E  Doyle) 

SHEl  I Oil  CO  HOUSTON.  TX  (MARSHALL) 

SOU  TH  AMERICA  N.  Nouel.  Valencia.  Venezuela 
SWEDEN  GewTevh  Inst;  VHH<I  rbrary).  Stivkholm 
TIDE  W ATER  CONSTR  CO  Norfolk  V A (Fowler) 

TRW  SYSTEMS  Cl  KVEl  AND.  OH  (ENG  LIB.);  REDONDO  BEACH.  CA(DAI) 

UNITED  KINGDOM  Cement  A Concrete  Assoc  Wcxham  Spnngs.  Slough  Bucks;  Cement  A Concrete  Assoc,  (l-it. 
Ex).  Bucks.  D.  New.  G.  Maunsell  A Partners.  I ondon.  I ibrary , Bristol;  Shaw  A Hatton  (E  Hansen).  London; 
Taylor.  Wnodrow  Constr  (0I4P).  Southall.  Middk  s«"x.  Unis  of  BrtstoltR.  Morgan).  Bnstol 
W ATT  BRIAN  ASStX'  INC.  Houston.  TX 

WESTINGHOUSK  El  KCTRIC  CORP.  Annapolis  MDtOceamc  IXv  l ib.  Bryan);  l ibrary.  Pittsburgh  PA 
WISS.  JANNEY.  El.STNER.  A ASStX'  Northbrook,  It  (D  W Pfeifer) 

WtX)DW  ARD-CLYDK  CONSULTANTS  (A.  Harngan)  San  Francisco;  PLYMOUTH  MEETING  PA  (CROSS,  III) 

AL  SMIXVTS  Los  Angeles.  l'A 

BARA.  JOHN  P I akew.wxl.  CO 

BROW  N.  ROBERT  University.  Al 

BULLOCK  la  Canada 

E HEUZE  Boulder  CO 

CAPT  MURPHY  Sunnyvale.  CA 

GREG  PAGE  EUGENE.  OR 

R E BESIER  OW  Savbrook  CT 

f W MERMEl  Washington  DC 


